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"Streuung an einem Kastenpotential 


Von D. Geißler 


‘ 


Inhaltsiibersicht 
Die S-, P- und D-Phasenverschiebungen ö,, 6,, 6. bei der Streuung zweier 
Teilchen mit gleicher Masse aneinander, deren Wechselwirkung durch ein 
Kastenpotential beschrieben wird, werden mit der gewöhnlichen und mit der 
relativistischen Schrédinger-Gleichung berechnet und graphisch dar- 
gestellt. Dabei dient 6 als Parameter und eine Kombination von (im wesent- 
lichen) Potentialtiefe und -reichweite erscheint als Funktion einer Kombi- 
nation von Reichweite und Energie der einfallenden Teilchen. Einige Mög- 
lichkeiten der Anwendung der Kurven 6 = const werden diskutiert. 


Einleitung 

Es ist bisher noch nieht möglich, die Form der Wechselwirkung zwischen 
zwei Nukleonen eindeutig festzulegen. In der phänomenologischen Theorie 
der Kernkräfte werden verschiedene Potentiale diskutiert, darunter das Kasten- 
potential, das sich dadurch auszeichnet, daß damit zahlreiche Probleme relativ 
einfach analytisch behandelt werden können. 

Eine wichtige und experimentell sehr häufig angewendete Methode zur 
Untersuchung der Eigenschaften der Kernkräfte ist die Streuung von Nu- 
kleonen an Nukleonen. Beobachtet wird der differentielle oder totale Wir- 
kungsquerschnitt, der theoretisch durch die Phasenverschiebungen, die die 
einfallenden Wellen am Kernpotential erfahren, bestimmt ist. 

Es erscheint daher nützlich, die durch ein Kastenpotential hervorgerufenen 
Phasenverschiebungen zu berechnen und in einer Form darzustellen, die die 
Abhängigkeit von der Energie der einfallenden Teilchen und von den Poten- 
tialparametern auf einfache Weise abzulesen gestattet. Das wird in der vor- 
liegenden Arbeit unternommen, indem Kurven angegeben werden, bei denen 
die S-, P- und D-Phasen jeweils als Parameter dienen und eine Kombination 
von (im wesentlichen) Potentialtiefe und -reichweite als Funktion einer Kom- 
bination von Reichweite und Energie der einfallenden Teilchen aufgetragen ist. 

Da in letzter Zeit Streuexperimente bei hohen Energien im Vordergrund 
stehen, erschien es angebracht, die Massenveränderlichkeit zu berücksichtigen. 
Das wird dadurch getan, daß neben der gewöhnlichen auch die relativistische 
Schrédinger-Gleichung zugrunde gelegt wird. 


$ 1. Anschlußbedingungen 


Die Schrödinger-Gleichung für die Streuung eines Teilchens an einem 
anderen mit gleicher Masse lautet nach Abseparation der Schwerpunkts- 
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bewegung 
Hierbei Bedienen A den Laplace-Operator bezüglich der Relativkoordinessh 
M die Masse, E’ die Gesamtenergie der Partikeln im Schwerpunktssystem 
und V(r) die als kugelsymmetrisch angenommene potentielle Energie. 


Für ein Kastenpotential 
r<a 


gewinnt man über die Substitution RE een 


y= (cos @) 


Vir) = 


und die Lösung der Differentialgleichung für wu, die zur Bestimmung der 
Phasen 6, dienenden mf bei r = a: 
ph (&) — sin - («) 1 
MB 
B=2E’ = Energie im Laborsystem, 


Die Funktionen j, (8) und n, (ß) sind definiert durch!) 


wobei unter J, die Besselfunktion n-ter Ordnung zu verstehen ist. Die 


Striche in Gl. (1) bedeuten Ableitungen nach dem Argument. Durch die 
Formel!) 


ji (B) = (B) — jx (B) 
läßt sich Gl. (1) für 2 > 0 schreiben als 


pie 1 (B) _ xo — sin 6; (x) 
(B) cos 0, 7, (&) — sin (&) 


Die relativistische Schrödinger-Gleichung des Problems heißt 
1 ” 
Ay + —- 


mit BE” =2Me (1 + B/2 M c?)t, woraus genau wie oben die Anschlub- 


bedingungen folgen, die dieselbe Form wie G1. (1) bzw. Gl. (2) haben nur a 


jetzt AR 


= = + (2 B" + V,) 

Dabei gilt übrigens für B<Me, V,< Me: 
> 


1) Siehe z.B. L.I. Schiff, Quantum Mechanics, New York, McGraw-Hill Book 
Comp., 1949. 
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Die Funktionen 7,(6) und n,(ß) lassen sich bekanntlich!) in elementarer 
Weise durch ß, sin ö und cos ß ausdrücken. Tut man dies, so nehmen die An- 
schlußbedingungen für / = 0, 1 und 2 folgende Gestalt an: 


I=0: & cotg (x + do) = Beotg ß, (3) 

2/1 —« cotg (x + = {1 —B cotg}, (4) 

l= 2: 3a-*—{1— (x + 6,)}-! = 38 * —{1—f cotg (5) 

Die Vieldeutigkeit von 6 um ganzzahlige Vielfache von z in den Gln. (3) 

bis (5) wird durch die Verabredung beseitigt, daß fir 6 = a, d.h. wenn kein 

Potential vorhanden ist, 6 = 0 sein soll und daß von diesem Werte aus 6 

für wachsendes 9? — x?, also bei Verstärkung des Potentials, stetig zunehmen 
möge. 

§ 2. Bereehnung der Kurven 6 = const 


Zur Darstellung der Kurven 6 = const werden als Abszisse die GréBe 


A 


(nicht-relativistisch) 


(relativistisch) 


Es ist zu bemerken, daß der relativistische Ausdruck für 6? — «? im Falle 
B, V,< Mc? in den nicht-relativistischen übergeht, während x? für beide 
Fälle genau übereinstimmt. 

Um bei gegebenen « und 6 das zugehörige zu finden, müssen die Gln.(1.3), 
(1.4) und (1.5) gelöst werden. Dazu wurden Zahlentafeln für die rechten 
Seiten dieser Gleichungen berechnet, und zwar für 6 = 0 bis 10 in Schritten 
von 0,01 und für 8 = 10 bis 17,3 in Schritten von 0,12). Der gesuchte 6-Wert 
ergibt sich durch Interpolation. Die Tafeln sind so eingerichtet, daß im all- 
gemeinen lineare Interpolation ausreicht, um ß auf 4 Stellen nach dem Komma 
genau zu erhalten (für 6 < 10; für größere 5 auf 3 Stellen). Eine Ausnahme 
bilden die Stellen, wo der Funktionswert sehr groß ist. Zur Berechnung dieser 
ß-Werte wurden spezielle Näherungsformeln benutzt. 

Die Ergebnisse finden sich für größere Werte von «? in den Abb. 1 bis 6, 
für kleinere (mit Näherungsformeln berechnet) in den Abb. 7 bis 13. 


§ 3. Einige Anwendungen der Kurven 6 = const 


1. Es mögen kinetische Energie B der einfallenden Teilchen, Reichweite a 
des Potentials und Phasenverschiebung 6 gegeben sein, während die Potential- 
tiefe V, gesucht ist. 

Die durch % = a? = M Ba?/2h? festgelegte Abszisse liefert auf der vor- 
geschriebenen Kurve 6 = const einen Ordinatenwert 7 = ?— x, aus dem 
V, zu errechnen ist: 


= (nichtrelativistisch) 


*) Diese Tafeln werden an anderer Stelle veröffentlicht. er ung 


— 
‘ 
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Vo = — + 2 B Me’ 
mit gegebenem E” = 2 Mc? (1 + B/2 Me2)} | (relativistisch) 


Min? 


a! 
6 8 


1 2 4 6 6 1 2 | 4 
Abb.1. Die Funktion f? — «? = f(x?) für ö, = const [definiert durch Gl. (1.3)]; 02 
V,(2E” + V,)a? (relativistisch); «= Reichweite, V, = Tiefe des Potentials, B= 
Energie im Laborsystem, Z” = 2 Mc? (1+ B/2M c?)t. Die gestrichelten Linien geben 

die Lage der Extrema an 


Fr V, a? (nichtrelativistisch), P?— « 


U6 


| maf 

2 4 6 8 mw 2 6 8 10% 2 4 6 

Abb. 2. Die Funktion — «? = f(a*) für 6, = const [definiert durch Gl. (13)]; 
3,2 < 6) < 6,4. Bezeichnungen wie bei Abb. 1 
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Abb. 3. Die Funktion — «? = f(a?) fiir 6, = const [definiert durch Gl. (1.4)]; 
0,05 < 6, < 0,5. Bezeichnungen wie bei Abb. 1 


T 


24 


\ 
\ 
\ 

\ 


2 
it ow 
1 2 6 6 8 10 2 4 6 8 10° 2 4 6 810 


Abb. 4. Die ; Funktion 6? — a? = f(a?) für 6, = const [definiert durch Gl. (1.4)]; 
0,4< 6, < 3,2. Bezeichnungen wie bei Abb. 1 
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Abb. 5. Die Funktion f? — «? = f(x?) für 6, = const [definiert durch Gl. 
3,2 <ö, < 6,4. Bezeichnungen wie bei Abb. 1 
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1 | 1 

2 6 8 10 2 1 6 
Abb. 6. Die Funktion — = f(a?) für const [definiert durch Gl. (1.5)) 
0,10 < 6, < 0,50. Bezeichnungen wie bei Abb. 1 


2. Gegeben seien B, V, und 6, gesucht die Reichweite a. 
Im nichtrelativistischen Fall ist zur Lösung die Gerade { 


y = (2 V,/B) x 
zu zeichnen, die die durch 6 = const bestimmte Kurve y = K (x) in einem 


Punkt mit der Abszisse schneidet. Dann ist 
(2h?/MB 
Im relativistischen Fall muß man die Abszisse % des Schnittpunkts der 
Geraden 
1 1 MB 


Yy == . the Vo (2 _ + V,) x, ]:2 
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Abb. 7. Die Funktion f? — a? = Abb. 8. Die Funktion B? — a? = f (a?) 
f(a?) fiird, = const und kleine fiir 6, = const und kleine Werte von 
Werte von a? [definiert durch x2 [definiert durch Gl. (1. 3)]; 
Gl. (1.3)]; 0,2 < 6) < 2,8. Bezeich- a < 09 < 6.0. Bezeichnungen wie 
nungen wie bei Abb. 1 bei Abb. 1 =) 
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mit der gegebenen Kurve y = f(x) bestimmen. Dann ist wieder Be. 


2 


3. Gegeben seien B, 6, und es soll V, in Abhängigkeit von a Auslesen 
werden. 

Die Abszissenskala läßt sich sofort in eine a?-Skaia umändern, und zu jedem 
Wert a ergibt sich wie bei 1. auf der durch das gegebene 6 festgelegten Kurve 
y = f(x) ein y-Wert, aus dem sich dann V, berechnet. 

In Abb. 14 sind einige a—V,-Kurven für B = 100 MeV und verschiedene 
Werte von ö, im relativistischen Fall gezeichnet. 


a? = 


— 
= 


4. Gegeben seien zwei Energien B, B und zwei zugehörige Phasen 6, ö. 
Es soll ein Potential a, Vj, das bei der Energie B die Phase 6, bei B die Phase ö 
liefert, bestimmt werden. 

a) Nichtrelativistisch. 

Wegen B/# = 2 h?/Ma® = Blx ist x/% = eonst,d.h. die Differenz zwischen 


den Abszissen # und x, die zu den verschiedenen Energien B und B, 
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aber zum gleichen Wert von a gehören, ist bei einer logarithmischen 
konstant. Wo der horizontale Abstand der Kurven 6 = 6 und 6 = 6 gleich 
dieser Differenz ist (dabei ist noch die richtige Zuordnung zu beachten: £ auf }, 


= y e 
x auf 6), dort sind 7 = =, V,a? und # = a... a? abzulesen, woraus sich ne 


Abb. 9. Die Funktion — a? = f («*) 
für 56) = const und kleine Werte 0m am un 
von «* [definiert durch Gl. (1.3)]; Abb.10. Die Funktion B? — a? = f(a?2) für 6, = cons 
ö, = 2 und 6,4. Bezeichnungen wie und kleine Werte von a? [definiert durch Gl. (1.4)} 
bei Abb.1 0,05 < 6, < 2,80. Bezeichnungen wie bei Abb. ! 


b) Relativistisch. 
Auch hier gilt > 
x/# = B/B = const er 2. 
für gleiches a. Daneben ist has get ae, 
Vy = — BY +(E"2+2BMe g/z)', 
N vo Vo — (E’2 + 2 BM«& y/z)*. 
Durch Vereinigung der beiden letzten Gleichungen entsteht mit Beriick- 
sichtigung von (1) 


(2) 
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Zu einem vorgegebenen ä£ gehört vermittels der Kurve 6 = 6 ein bestimm- 
tes 7, nach (1) ein bestimmtes % und dazu vermittels der Kurve 6 = 6 ein Y- 
Durch Wahl von # sind demnach 7 und y festgelegt, so daß die linke Seite 
von (3) als Funktion von £ gezeichnet werden kann. Ihre Nullstellen liefern 
die gesuchten #-Werte, aus denen sich nach 


= M Ba®/2 
a errechnen läßt. und der auf der Kurve 6 = 6 zugehörige Wert y legen 
nach (2) V, fest. 
. Es sollen Potentiale, die für eine gegebene Energie B eine Phase 
i’ ” ö< 6” und für eine andere gegebene Energie B eine Phase 6’ < 6 < 6” 
liefern, bestimmt werden. 6’, 6’, 6’, 6” liegen fest. 

Nach 4. ist der Abstand zwischen den Abszissen # (zu B gehörig) und 

7 (zu B gehörig) bei gemeinsamem a in logarithmischer x-Skala konstant: 
const. 
Im nichtrelativistischen Fall gehört ein Punkt %, 7 dann zu einem der gesuch- 
ten Potentiale, wenn gleichzeitig 

a) der Punkt #, 7 zwischen den Kurven 6 = 6’ und 6 = 6”, 

b) der Punkt x + A, y zwischen den Kurven 6 = 6’ und 6 = 6” liegt. 
Aus den so ausgewählten Z- und 7-Werten folgen mittels des gegebenen B 
die zugehörigen Größen von a und Vp. 

Im relativistischen Fall sind die Kurven ö = 6’ bis 6 = 6” am besten in 
a—V,-Kurven umzuzeichnen (siehe 3.). Die geeigneten Potentiale lassen sich 
dann sofort ablesen. Dasselbe Verfahren kann man natürlich auch bei nicht- 
relativistischer Betrachtung anwenden. 

6. Neue ö-Kurven lassen sich in der üblichen Weise durch Interpolation 
aus den Abb. 1 bis 13 gewinnen. In Abb. 1 ist die Kurve für 6, = 0,73 durch 
Interpolation bestimmt worden. 

7. Es soll die Abhängigkeit der Phase ö von der Energie B der einfallenden 
Teilchen bei gegebenen Potentialabmessungen a, V, untersucht werden. 


Wegen der größeren Übersichtlichkeit wird nur der nichtrelativistische 
Fall betrachtet. 


Aus der Entwicklung der Funktionen j, (x) und n, («) für kleine «!) Be; | 
a! 1-3-5---(21—1) 
folgt, daß die Anschlußbedingungen (1.1) bzw. (1.2) für x = 0, d.h. wegen 
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(Bo) 
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Die durch die Gln. (4) und (5) definierten Werte von V¥ sind diejenigen 
Potentialtiefen, bei denen bei fester Reichweite a gerade ein neuer gebundener 
Zustand (Bindungsenergie 0) mit der Drehimpulsquantenzahl / auftreten 
kann. 

Nun haben aber, wie sich zeigen läßt, die Gln. (4) genau dieselbe Form wie 
die Gleichungen, die die Minima der Kurven 6 = const bestimmen, nur daß 


dort ß’ an Stelle von {4 steht. Wenn 8’ eine Lösung einer der Minimums. 
bedingungen und ß, eine Lösung einer der Gln. (4) mit demselben Wert von! 
ist, so gilt 8’? = ßg* oder 

Vy + (B/2) = kinetische Energie im Schwerpunktssystem = V#, 
wo B die Energie (im Laborsystem), die zu dem betreffenden Minimum ge- 
hért, bedeutet. 

Bei Vergrößerung von B bewegt man sich bei festem a und V, auf einer 
Parallelen zur Abszissen- 
achse im Abstand y= 
on M V, a*/ii? = const nach 
334 rechts. Dabei wird der 
Moment kommen, wo zum 
ersten Male die Gerade 
y = const eine Kurve 6= 
6) in deren Minimum be- 
rührt. Dieser Stelle ent- 

spricht ein Maximum der 
Pr Phase als Funktion von B, 
i 2% wie aus den Abb. 1 bis6 
Mit den obigen Erkennt- 
382 nissen zusammen läßt sich 
das Ergebnis etwa so for- 
mulieren: 


Die Phasenverschie- 
bung 6 hat als Funktion der 
Energie B der einfallenden 


bei festen Potentialab- 
messungen a, V, bei den- 


kinetische Energie V, + 


(B/ 2) im Schwerpunkts- 
Abb. 11. Die Funktion — a? = f (a2) für 6, = const d 
und kleine Werte von a? [definiert durch Gl. (1.4)]; nes Teilchens ist, as von 
a < 6, < 6,0. Bezeichnungen wie bei Abb. 1 einem Potential mit dersel- 
ben Reichweite a gerade 
noch gebunden werden kann (mit der Bindungsenergie 0). Die Werte B = 8 
definieren eine Art von virtuellen Zuständen, in dem Sinne nämlich, daß die 
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Abb. 12. Die Funktion 6? — a? =f («*) 
für 6, = const und kleine Werte von 
& [definiert durch Gl. (1.4)]; 6, = 27 
und 6,4. Bezeichnungen wie bei 

bb.1 
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Abb. 13. Die Funktion 6? — a? = f (a?) fiir 6, = 
const und kleine Werte von a? [definiert durch 


Gl. (1.5)]; 0,10 < 4, 


< 0,50. Bezeichnungen wie 
bei Abb. 1 
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Abb. 14. Kastentiefe V, als Funktion der Reichweite a für die Teilchenenergie (im Labor- 
system) B = 100 MeV und verschiedene Werte von ö, 
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Wellenfunktion der Teilchen dort stärker verzerrt wird als in benachbarten 
Zuständen. Mit jedem weiteren Minimum, das die Parallele y = const bei 
einer Kurve 6 = ö(M berührt, d.h. mit jedem weiteren Phasenmaximum, 
wandert ein neuer Zustand in den ständig tiefer werdenden Kasten a, V* 
ein. Dazwischen liegen Minima der Phase, und zwar dort, wo die Gerade 
y= const eine ö-Kurve in deren Maximum berührt. Bei Erhöhung der 
Energie rücken Maxima und Minima einander immer näher. 


Ich bin Herrn Prof. Kockel für die Anregung dieser Untersuchungen und 
für sein förderndes Interesse während der Ausarbeitung zu großem Dank 
verpflichtet. 


Theoretisch-Physikalisches Institut der 


Bei der Redaktion eingegangen am 3. Dezember 1955. x 
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Inhaltsübersicht 


Es wird ein Überblick über die Ergebnisse und die phänomenologische 
Theorie von Proton—Proton-Streuexperimenten gegeben. Die Hypothese 
der Protonenisomere mit der weitergehenden Annahme, daß die einfallenden 
Protonen schon vor dem eigentlichen Streuakt als unterscheidbar von den 
ruhenden anzusehen sind, wird zur Berechnung des differentiellen Streuquer- 
schnitts bis einschließlich der D-Phasen für drei Krafttypen unter Berück- 
siehtigung der Coulombkräfte verwendet. Es zeigt sich, daß mit gewissen 
schmalen und sehr tiefen Kastenpotentialen, insbesondere solchen vom Serber- 
Typ, die Experimente zwischen 100 und 435 MeV befriedigend beschrieben 
werden können. 


§ 1. Streuversuche bei niederen Energien _ 

Bei der theoretischen Deutung des Aufbaus der Atomhülle ging man aus 
von der Coulombkraft, die die Wechselwirkung zwischen Kernen und Elek- 
tronen vermittelt, und suchte nach einer Mechanik, der die Teilchen gehorchen. 
Als Resultat ergab sich die Quantenmechanik. In der Physik der Atom- 
kerne gilt es, das umgekehrte Problem zu lösen. Man weiß oder kann mit 
ziemlicher Sicherheit annehmen, daß die Nukleonen quantenmechanisch be- 
schrieben werden können, während die Eigenschaften der Kernkräfte zu be- 
bestimmen sind. 

Eine der aufschlußreichsten Möglichkeiten, die Natur der Wechselwirkung 
zwischen Nukleonen zu untersuchen, stellen Streuexperimente dar. Unter 
diesen steht die Streuung von Neutronen oder Protonen an Protonen an hervor- 
ragender Stelle, da hier die einfachsten Nukleonensysteme, die neben dem 
Deuteron möglich sind, auftreten. 

Mangels stärkerer Beschleuniger begannen die Streuversuche im Gebiet 
niederer Teilchenenergien (bis 10 MeV). Ihre Auswertung ergab in Verbindung 
mit den Erkenntnissen, die auf anderen Gebieten der Kernphysik gewonnen 
wurden, folgende Resultate: 

1. Die Kernkräfte lassen sich durch verschiedene Potentiale (Kasten- 
potential, Yukawapotential u. a.), wenn diese nur von kurzer Reichweite 
vom phänomenologischen Standpunkt aus gut 


. Die Kernkräfte sind spinabhangig. 


arten 
t bei 
num, 
7 
t, Vy 
Trade 
der 
Er 
hac’, 
. 
: 
N 
= 


Annalen der Physik. 6. Folge. Band 18. 1956 


3. Die Kernkräfte sind nicht rein kugelsymmetrisch, sondern bestehen ; 
zu einem gewissen Teil aus Tensorkräften, die unter anderem Anlaß zu dem 
beobachteten Quadrupolmoment des Deuterons geben. 

4. Das Phänomen der Sättigung, das sich in der Konstanz der Bindung 
energie pro Nukleon bei schweren Kernen bemerkbar macht, kann dureh 
die Annahme, daß die Kernkräfte aus einer Mischung verschiedener Arten 
von Austauschkräften zusammengesetzt sind, erklärt werden. 

5. Die Kernkräfte können als ladungsunabhängig angesehen werden, d.h, 
die Kräfte zwischen einem Neutron und einem Proton sind die gleichen wie 
die zwischen zwei Protonen (abgesehen von der Coulombschen Abstoßung im 
letzteren Fall). 

$ 2. Streuung hochenergetischer Nukleonen 

Mit dem Bau starker Beschleunigungsmaschinen, die Teilchen mit meh- 
reren hundert MeV Energie liefern, in Berkeley, Rochester und an anderen 
Orten eröffneten sich neue Möglichkeiten. Insbesondere erhoffte man, näheren 
Aufschluß über die genaue Form des Kernpotentials zu erhalten, da die der 
Erhöhung der Energie entsprechende Verkleinerung der de Broglie-Wellen- 
länge es den Partikeln ermöglichen sollte, weiter in die Region der Kern- 
kräfte einzudringen. 

Wohl kam man der Beantwortung dieser Frage etwas näher, indem sich 
herausstellte, daß relativ langschwänzige Potentiale, wie das Yukawa- oder 
Exponentialpotential, etwas besser mit den experimentellen Daten überein- 
stimmen als kurzschwänzige (z. B. Kastenpotential), doch mußte man dafür 
eine Anzahl neuer ungelöster Probleme in Kauf nehmen. Die Ergebnisse 
waren nämlich recht unerwartet und zeigten nur eins sehr deutlich, daß eine 
Extrapolation der Schlußfolgerungen, die man aus den Versuchen bei niederen 
Energien gezogen hatte, auf die Verhältnisse bei hohen Energien nicht ohne 
weiteres gestattet ist. Bei der theoretischen Auswertung der n-p-Streu- 
experimente fanden Christian und Hart!), daß sich am besten eine Aus- 
tauschkraft vom Serber-Typ bewährt. Zu dieser gehört das Potential 

V=3(1+ V 
wobei PM den Majorana-Operator bedeutet, der die Ortskoordinaten der 
Partikeln vertauscht, ohne den Spin mitzunehmen. Diese Form der Wechsel- 
wirkung führt aber, wie Gerjuoy?) gezeigt hat, nicht zu einer Sättigung bei 
schwereren Kernen, vorausgesetzt, daß nur Zweikörper-Kräfte auftreten. 

Eine andere wesentliche Schwierigkeit, auf die das bei hohen Energien 
gesammelte Erfahrungsmaterial hinzudeuten scheint, besteht darin, daß 
die Kräfte zwischen einem n-p-Paar anders zu sein scheinen als die zwischen 
zwei Protonen, daß man also vielleicht die Hypothese der Ladungsunabhängig- 
keit der Kernkräfte aufgeben muß. 


$3. p—p-Streuung bei hohen Energien 

p—p-Streuversuche bei hohen Energien wurden durchgeführt von Birge 
und anderen?) bei 75 und 105 MeV, von Cassels und anderen*) bei 147 MeV, 

1) R.S. Christian u. E. W. Hart, Physic. Rev. 77, 441 (1950). 

2) E. Gerjuoy, Physic. Rev. 77, 568 (1950). 

3) R. W. Birge, U. E. Kruse u. N. F. Ramsey, Physic. Rev. 88, 274 (1951). 

4) G.M. Cassels, T. G. Pickavance u. G.H. Stafford, Proc. Roy. Soc. London 
(A) 214, 262 (1952). 
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von Chamberlain und anderen°)®) bei Energien zwischen 120 und 350 MeV, 
von Oxley und Schamberger’) bei 240 MeV, von Towler®) bei 240 MeV, 
von Mott und anderen?) bei 435 MeV und von Marshall und anderen”) bei 
144, 271 und 429 MeV. Wenn sich die Ergebnisse auch zum Teil in der abso- 
luten Größe des differentiellen Wirkungsquerschnitts widersprechen, so 
zeigen doch alle bis auf eine Ausnahme (Mott und andere®)) die völlig uner- 
wartete Eigenschaft, daß der Streuquerschnitt innerhalb der Fehlergrenzen 
und abgesehen von sehr kleinen Winkeln im Schwerpunktssystem nicht vom 
Streuwinkel abhängt. Man hätte vielmehr annehmen sollen, daß bei derart 
hohen Energien der Beitrag der P-, D- usw. Phasen zur Streuung so stark 
wird, daß eine merkliche Anisotropie der Winkelverteilung resultiert. 

Der erste Eindruck, den eine solche kugelsymmetrische Verteilung ver- 
mittelt, ist der, daß es sich um reine S-Streuung handeln muß. Der differen- 
tielle Wirkungsquerschnitt pro Raumwinkel im Schwerpunktssystem berechnet 
sieh in diesem Fall in bekannter Weise zu td 

a (0) = sind, = B, 
wo ö, die vom Kernpotential an der Partialwelle mit Z = 0 hervorgerufene 
Phasenverschiebung und B die Energie der einfallenden Protonen im Labor- 
system bedeuten. Von den Coulombkräften, die bei hohen Energien gegen- 
über den Kernkräften zurücktreten, ist hierbei abgesehen. Für B = 350 MeV 
ist der Maximalwert für o (©), der für 6, = 2/2 eintritt, gleich 1/k2 = 2,4. 10-27 
em? = 2,4 mbarn. Der experimentelle Wert dagegen liegt bei etwa 3,5 mbarn 
oder höher. Die einfache Annahme reiner S-Streuung ohne zusätzliche Hy- 
pothesen scheitert damit von vornherein. 

Es gibt prinzipiell zwei Möglichkeiten, die Isotropie der Streuung zu inter- 
pretieren. 


§ 4. Deutung der Isotropie durch Superposition 


Man kann die beobachtete Isotropie dadurch zu erklären versuchen, daß 
sich Partialwellen mit verschiedener Drehimpulsquantenzahl / so überlagern, 
daß gerade eine Winkelunabhängigkeit resultiert. Man wird diesen Versuch 
zunächst sogar mit der Absicht verbinden, die Hypothese der Ladungsunab- 
hängigkeit der Kernkräfte aufrechtzuerhalten und die Potentiale so zu 
wählen, daß auch die Daten bei niederen Energien damit richtig wiedergegeben 
werden. Zwar verhalten sich die Funktionen o,, (©) und o,, (©) bei hohen 
Energien ganz verschiedenartig: o,, ist deutlich anisotrop und dem Betrag 
nach relativ niedrig, beides den Erwartungen entsprechend; o,, ist dagegen 
isotrop und liegt recht hoch, beides den Erwartungen widersprechend. Dazu 
kommt, daß o,, (90°) nur sehr wenig mit der Energie variiert. Trotzdem hat 
der Versuch einige Aussichten auf Erfolg; denn beim n—p-System treten 
Zustände auf (38, 1P,...), die beim p—p-System verboten sind und die die 


5) O. Chamberlain, E. Segré u. C. Wiegand, Physic. Rev. 83, 923 (1951). 

6) O. Chamberlain, G. Pettengill, E. Segré u. C. Wiegand, Physic. Rev. 98, 
1424 (1954); O. Chamberlain u. John D. Garrison, Physic. Rev. 95, 1349 (1954); 
0. Chamberlain: Mitteilung in ™). 

7) C. L. Oxley u. R. D. Schamberger, Physic. Rev. 85, 416 (1952). 

8) O. A. Towler jr., Physic. Rev. 85, 1024 (1952). 

®) W. E. Mott, R. B. Sutton, J. G. Foxu. J. A. Kane, Physic. Rev. 90, 712 (1953). 

10) J. Marshall, L. Marshall u. V. A. Nedzel, Physic. Rev. 92, 834 (1953), 
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n—p-Streuung genügend stark beeinflussen können. Außerdem haben di 
Zustände des p—p-Systems gemäß dem Pauli-Prinzip andere Gewichte in 
n—p-System. Dies könnte genug 'Spielraum liefern, um beide Erfahrung: 
reihen durch eine Kraftsorte zu erklären. 

Derartige Versuche wurden unternommen von Christian und Noyesll, 
Case und Pais!?) und von Jastrow®). 

Christian und Noyes?) stellen fest, daß durch ein glattes, nichtsinge 
läres, monotones (d.h. überall anziehendes oder überall abstoßendes) Sip- 
gulettpotential, wie es zu den Experimenten bei niederen Energien paßt, 
stets die Streuung in Vorwärtsrichtung stark bevorzugt wird; denn bei Sin. 
gulettzuständen treten im p—p-System nur Partialwellen mit geradem / auf. 
Die zugehörigen Kugelfunktionen P, (cos @) [9 = Streuwinkel im Schwer- 
punktssystem] haben bei © = 0° alle den Wert 1, während sie für 0 = 9 
ihr Vorzeichen wechseln: P, (0) = 1; P, (0) = —0,5; P, (0) = 0,375.. Da 
bei einem monotonen Potential alle Phasen gleiches Vorzeichen haben, ent- 
steht bei der Überlagerung der einzelnen Wellen im Streuquerschnitt bei 0? 
konstruktive, bei 90° im wesentlichen destruktive Interferenz, vorausgesetzt, 
daß die Energie hoch genug ist, um Phasen mit höheren /-Werten beitragen 
zu lassen. Damit ist das Vorwärts-Maximum erklärt. 

Damit der angestrebte flache Verlauf der Funktion o,, (©) erreicht wird, 
müßte man ein nicht-monotones Potential so wählen, daß sin ö, für gerade | 
das Vorzeichen wechselt wie die zugehörigen Kugelfunktionen. Dadurch be- 
käme o,, (90°) einen hohen Wert. Wenn dann ein zentrales Triplettpotential, 
das wegen der Antisymmetrie der Ortsfunktion stets bei@ = 90° verschwindet, 
hinzugefügt wird, könnte o,, (9) = eonst annähernd eintreten. Nach 
Christian und Noyes scheint es aber nicht möglich zu sein, ein derartiges 
Potential, das auch noch die Daten bei niederen Energien richtig wiedergibt, 
zu finden. 

Es bleibt der Ausweg, ein Tensorpotential heranzuziehen. Ein solches 
führt nämlich nicht dazu, daß der daraus sich ergebende Wirkungsquer- 
schnitt in Triplettzuständen allgemein für © = 90° verschwindet; denn die 
Antisymmetrie der Ortsfunktion braucht hier nicht durch ©, sondern kann 
auch durch den azimutalen Winkel ® bewirkt werden. Es besteht demnach 
die Aussicht, die Bevorzugung der Vorwärtsstreuung durch das Singulett- 
potential auszugleichen. 

Christian und Noyes zeigen durch ein allgemeines Argument, daß die 
Interpretation der p—p-Streuung bei 340 MeV verlangt, das Potential im 
Gebiet von 0 bis 0,5 - 10-13 cm recht stark zu machen. Sie wählen infolgedessen 
eine singuläre Tensorwechselwirkung. Das Gesamtpotential hat die Gestalt 


Triplett: V = V, S,, rir)?’ 

_ 


R. Christian u. H. Noyes, Physic. 
12) K.M. Case u. A. Pais, Physic. Rev. 
18) R. Jastrow, Physic. Rev. 81, 165 (1951). 
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Die Rechnungen werden durchweg in der Bornschen Näherung durch- _ 


geführt und zeigen, daß sich mit den Parameterwerten 
r, = 2,6-10-% em; r, = 1,6-10"®em; V, = 13,27 MeV; V, = + 18 MeV 


die p—p-Streuung bei 32 MeV und bei 340 MeV erklären läßt, daß es aber — 
nicht möglich ist, mit demselben Potential die n—p-Streuung bei 40 und 280 MeV 
zu deuten. Die Verfasser schließen daraus, daß man bei Verwendung stati- _ 
scher Potentiale die Hypothese der Ladungsunabhängigkeit der Nukleonen- 
wechselwirkung aufgeben muß. 5 

Im Gegensatz dazu kommen Case und Pais") zu dem Ergebnis, daß man © 
mit einem nicht-statischen Zusatzpotential vom Typ 
wo 2 den Operator des relativen Bahndrehimpulses und S = 9, + ö, die Ma- — 
trix des Gesamtspins bedeuten, die Ladungsunabhängigkeit der Kernkräfte 
beibehalten kann. Fügt man einen solchen Term in der Wechselwirkungs- 
energie hinzu, so bleiben die Verhältnisse bei kleinen Energien ungeändert, 
weil dort nur S-Zustände an der Streuung beteiligt sind, in denen die Zusatz- 
kraft wegen L = 0 nicht auftritt. 

Case und Pais meinen übereinstimmend mit Christian und Noyes, © 
daß eine starke Wechselwirkung in kleinen Abständen zur Erklärung der — 
p—p-Streuung nötig ist, und halten es für möglich, daß die Zustände 18, 
3P,... empfindlicher für die Vorgänge in kleinen Entfernungen sind als die 
Zustände 38, 1P,... Da bei der n—p-Streuung der 3S-Zustand mit seinem — 
relativ großen statistischen Gewicht stark beteiligt ist, könnten die 


Unterschiede zur p—p-Streuung auf dieser Basis erklärbar sein. Durch ares 


Rechnungen in Bornscher Näherung wird nachgewiesen, daß die (%- S)- 


Terme den Streuquerschnitt bei @ = 90° heben, was auf die Unterschiede _ 
in den drei Phasen für ?P,, ®P, und *P, zurückgeführt wird. Der allgemeine _ 


Verlauf von o,, (©) bei 350 MeV wird sich damit erklären lassen. b 

Wie sich zeigt, hängt das Verhältnis o (30°)/o (90°) bei 350 MeV sehr — 
empfindlich von der Form des Radialteils V (r) ab. Als bester Wert — sich _ 
a (30°)/o (90°) = 1,6 für ein hochsinguläres Zusatzpotential 


1 d 
dr/r,) ( r/r, ) 


wo r, = const die Reichweite angibt. Es gilt hier V (r) ~ r- für kleine r. ‘ 
S—D-Interferenzterme verursachen eine Anisotropie im n—p-Streu- | 


querschnitt, die dazu zwingt, das Vorzeichen des (% - ©)-Termes entgegen- _ 
gesetzt zu dem beim Schalenmodell (Mayer™), Haxel und andere™)) ge- : 


bräuchlichen zu wählen. Aus demselben Grunde ergeben sich Schwierigkeiten, — 
die 30 MeV- und die 350 MeV-p—p-Experimente zu vereinbaren. Case und 
Pais nehmen an, daß sich auch die n—p-Daten durch ihren theoretischen — 
Ansatz deuten lassen. Die Begründung dafür sehen sie in dem Umstand, daß 
nach überschlägigen Berechnungen die (2 - G)-Krafte im n—p-System nurein | 
Viertel des Beitrages im p—p-System zu o (90°) liefern. Einen solehen Unter- — 
schied etwa verlangt auch das Experiment. 


1) M. Goeppert er. Physic. Rev. 75, 1969 (1949); 78, 16, 22 (1950). 2 
4) O. Haxel, J. H. D. Jensen u. H. E. Suess, Naturwise. 85, 376 — Physic. 
Rev. 75, 1766 (1949); Z. Physik 128, 295 (1960). 5 
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Das Modell von Case und Pais hat den Vorzug, daß nur feldtheoretisch 
sanktionierte Kräfte auftreten. Insbesondere ist zu bemerken, daß eine 
Spin-Bahn-Kopplung sicherlich existiert und daß die Verwendung des ( - &). 
Operators beim Schalenmodell erfolgreich gewesen ist. Die Änderung des 
Vorzeichens darf nach Breit und Gluckstern!®) nicht zu tragisch genommen 
werden, da in einem schwereren Kern energetisch wesentlich andere Verhält- 
nisse herrschen als bei einem Streuexperiment. Diese Autoren sind der An- 
sicht, daß sich eventuell mit einer verallgemeinerten Wechselwirkung 


A + W, 
wo X, und X, vom Bahnmoment abhängen, Erfolge erzielen ließen. 
Jastrow33) geht von dem Gedanken aus, daß in der klassischen Physik 
eine starre Kugel eine isotrope Streuverteilung erzeugt. Diesem Gebilde ent- 


spricht eine Kraft, die von einem gewissen Abstand r, an sehr plötzlich ab- 
stoßend zwischen den Teilchen wirkt. Jastrow wählt folgendes Potential 


co r< 1% 
singu ett: Vos . N (1 PM), r> 7 


‘Triplett: V = {a + (1—a) PX + [b+ (1—}) P”) yS,,) Vo, 


PM bedeutet den Ortsvertauschungsoperator. Die Anzahl der Parameter 
betragt acht. Die starre Kugel, deren Radius ry schlieBlich zu r, = 0,6 - 10-8 em 
gewählt wurde, tritt merkwiirdigerweise nur in Singulettzuständen auf. 
Die geringe Größe von ry garantiert, daß sich die Wirkung der Abstoßung bei 
kleinen Energien noch nicht bemerkbar macht, da sich die Teilchen auch in 
S-Zuständen wegen der relativ großen de Broglie-Wellenlänge einander 
nicht so stark nähern können. 

Ein anziehendes Potential bewirkt positive, ein abstoßendes negative 
Phasenverschiebungen. Von einer gewissen Energie an wird die Wellen- 
länge so klein, daß sich die Protonen in S-Zuständen genügend nahe kommen, 
um die Wirkung der abstoßenden starren Kugel überwiegen zu lassen. Es 
resultiert eine negative 1S-Phase. In D-Zuständen dagegen verhindert der 
zusätzliche Zentrifugalterm 6/r? eine ausreichende Annäherung. so daß D- 
Teilchen im wesentlichen dem anziehenden Teil des Potentials ausgesetzt 
sind und die Phasenverschiebung positiv bleibt. Dasselbe gilt a fortiori für 
höhere Drehimpulse (G, /,...; tP-, IF-,... Zustände treten im p—p-System 
nicht auf). Der S—D-Interferenzterm im Wirkungsquerschnitt ist mit 
P, (cos ©) multipliziert. P, (cos ©) wechselt aber beiO » 55° das Vorzeichen. 
Der Interferenzterm ist also bei kleinen Winkeln negativ, in der Nähe von 
© = 90° dagegen positiv, wie es zur Beseitigung des Vorwärtsmaximums be- 
nötigt wird. Es ist notwendig, den S—D-Interferenzterm mit heranzuziehen, 
da eine reine S-Streuung an der harten Kugel wohl eine isotrope Winkel- 
verteilung, aber einen zu niedrigen absoluten Wert des Streuquerschnitts 
ergäbe. 

Durch das Verschwinden von P, (cos ©) bei © » 55° entsteht ein Minimum 
des Singulett-Wirkungsquerschnitts in der Nähe dieses Winkels, das aber 
durch Triplett-Streuung aufgefüllt wird. 


16) G. Breitu. R. L. Gluckstern, Ann. Rev. of Nucl. Science 2, 365 (1952). 
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Durch die Forderungen, die an das Potential gestellt werden, auch die 
Experimente bei niederen und die n—p-Streuung bei hohen Energien richtig ee 
wiederzugeben, sind von vornherein die meisten Parameter festgelegt. EE 
bleibt eine Variation von r, möglich. Der Wert ry = 0,6 - 10-3 em erweist oe =; 
sich als der beste. Eine Vergrößerung führt zur Störung der Isotropie bei 
hohen Energien, eine Verkleinerung verbessert zwar die Isotropie bei 340 MeV, 
doch passen die Ergebnisse dann nicht mehr zu den 30 MeV-p—p-Daten. 

Die Berechnungen, die zum Teil durch numerische Integration und zum 
Teil in der Bornschen Näherung erfolgen, zeigen folgendes: Die Überein. 
stimmung mit p—p-Experimenten bei 30 MeV ist recht gut, während die 
Kurve für o,, (©) bei 340 MeV wohl bei größeren Winkeln einigermaßen iso- 
trop ist, aber unterhalb von © = 25° zu stark ansteigt. Die o,, (9)-Kurven : 
passen fiir 90 und 260 MeV besser als fiir 40 MeV. 

Als einen Vorzug des Jastrow-Modells kann man ansehen, daß mit einer mS 
starren Kugel die Sättigung der Kernkräfte in schweren Kernen erklärbar coal Re 
wie Parzen und Schiff!”) untersucht haben, wenn auch der Radius r, = 
0,6- 10-13 em für diese Anwendung zu klein ist. Andererseits ist schwer zu ver- 
stehen, daß bei Jastrow Tensor- und Zentralpotential im Triplett- 
die gleiche Reichweite haben, während Arbeiten von Feshbach und Schwin- _ R 
ger’) und von anderen Autoren (siehe Breit und Gluckstern”)) auf einen aa F 
beträchtlichen Unterschied beider Größen hindeuten. ; 

Die einzige Eigenschaft, die die Krafttypen von Christian und Noyes, 
von Case und Pais und von Jastrow gemeinsam haben, ist die starke — 
Wechselwirkung zweier Protonen in kleinen Entfernungen voneinander. Im 
Falle des Jastrow-Potentials ist die entsprechende Mesonenmasse von 
Breit!?) auf etwa 450 Elektronenmassen geschätzt worden. Mi 

Die Phasen werden in allen drei Fällen wenigstens teilweise mit Hilfe der _ EU 
Bornschen Näherung berechnet. Da deren Gültigkeit bei hohen Singu- = 
laritäten keineswegs gesichert ist, haben Goldfarb und Feldman”) exaktere Me 
Rechnungen, die auf Variationsprinzipien und numerischer Integration be- 
ruhen, angestellt. Sie benutzen folgende Potentiale: 


Chr.u.N. Jastrow C.u. 
v | r< To V ern 


8-6 d 


Singulett: 


ert 


V, 


12 (rr)? trlr, \ r/r, 
r, (in 10-13 cm) 2,6 0,40 1 des 
#4 (im 10-12 om) 1,6 0,75 
V, (in MeV) 13,27 375 45,8 ‚8 
(in MeV) +18 —50,8 +29.8 rod 


Es zeigt sich, daB das Jastrowsche Modell der starren Kugel den erhöhten _ 
Anforderungen einigermaßen standhält. Der p—p-Streuquerschnitt. bleibt — 


") G. Parzen u. L. I. Schiff, Physic. Rev. 74, 1564 (1948). 

8) H. Feshbach u. J. Schwinger, Physic. Rev. 84, 194 (1951). 
) 


8) G. Breit, Intern. Conf. Nucl. Physics and Physics Fundam. Particles 102 (1951). 
L. J. B. Goldfarb u. D. Feldman, Physic. Rev. 88, 1099 (1952). = 
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ziemlich flach, wenn er auch zu tief zu liegen kommt. Im Fall der Potentiale 
von Christian und Noyes und von Case und Pais dagegen ist die Born- 
sche Näherung völlig unzulänglich. Bei Christian und Noyes stimmt 
oy, (9) nur für Winkel O > 45° mäßig mit der Erfahrung überein, während 
das (2 - G)-Potential zu einer großen Anisotropie, also zu qualitativ falschen 
Verhältnissen führt. Um eine Lösung zu garantieren, werden die Singu- 
laritäten in einer Entfernung von der Größenordnung der Compton-Wellen- 
länge des Protons (0,2 - 10-13 cm) abgeschnitten. Die Verfasser halten es 
nicht für ausgeschlossen, daß bei einer anderen Lage des Schnittes bessere 
Ergebnisse zu erzielen seien. 

Zu ähnlichen Folgerungen kommt Swanson2) bezüglich des Potentials 
von Christian und Noyes. Er findet ebenfalls, daß die exakte Rechnung 
die Übereinstimmung mit dem Experiment verschlechtert, und zwar wird die 
Anisotropie um so größer, je mehr an Volumen vom Potential abgeschnitten wird. 

Noyes und Gamnitz™) unternehmen einen Rettungsversuch, indem sie 
die Interferenz von Coulomb- und Kernkräften mitberücksichtigen. Die In- 
terferenzterme reichen jedoch nicht aus, um das Vorwärtsmaximum, das sich 
nach Goldfarb und Feldman?®) herausgestellt hatte, zu beseitigen. Sie 
schließen daraus, daß ein Potential von der Christian und Noyes-Form 
durch das Experiment widerlegt sei. 

Gleichzeitig weisen Noyes und Gamnitz darauf hin, daß das Jastrow- 
Modell unter Umständen sehr empfindlich auf bisher vernachlässigte kinema- 
tische relativistische Effekte reagiert. Eine Erhöhung der reduzierten Masse 
auf 1,06 (M/2) bringt bereits qualitative Änderungen im Streuquerschnitt mit 
sich, während die reduzierte Masse 1,18 (M/2) zu völlig anderen Ergebnissen 
führt. Der Grund liegt darin, daß sich bei hohen Energien im Singulettzustand 
die negative Phase der starren Kugel und die positive Phase des anziehenden 
Potentialteils zunächst fast kompensieren, eine Vergrößerung der reduzierten 
Masse aber nur die positive Phase wesentlich beeinflußt, so daß eine relativ 
große Änderung der Gesamtphase resultiert. 

Jede phänomenologische Theorie hat naturgemäß den Mangel, eine konse- 
quente relativistische Behandlung nicht zuzulassen. Während sich kinemati- 
sche Effekte noch ziemlich leicht berücksichtigen lassen (Case und Pais”), 
ist es nur schwer möglich, den genauen Einfluß dynamischer Effekte zu über- 
sehen (Snyder und Marshak*), Siegel“). Es hat daher wenig Sinn, eine 
größere Genauigkeit bei den Rechnungen anzustreben. 


§ 5. Deutung der Isotropie durch Protonenisomere 


Alle Versuche, die Isotropie der Winkelverteilung bei hohen Energien 
durch Superposition von winkelabhängigen Partialwellen zu deuten, führen 
zu etwas unbefriedigenden Ergebnissen. Es gibt nun noch eine andere Möglich- 
keit der Erklärung, die mit den von niederen Energien her gewohnten Vor- 
stellungen bricht. Sie wurde zuerst von Breit25)2) vorgeschlagen. 


21) Don R. Swanson, Physic. Rev. 89, 740 (1953). 

2) H. P. Noyes u. H. G. Gamnitz, Physic. Rev. 88, 1206 (1952). i 
23) H. Snyder u. R. E. Marshak, Physic. Rev. 72, 1253 (1947). ir 
%) A. Siegel, Physic. Rev. 82, 194 (1951). oh 
25) G. Breit, Physic. Rev. 84, 1053 (1951). 


%) G. Breit, Physic. Rev. 87, 208 (1952). NENNE 
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Breit bemerkt, daß bei Energien von der Größenordnung der Ruhenergie _ 
des Pi-Mesons der Streuprozeß die Natur der Protonen unter Umständen 
derart ändern könnte, daß sie hinterher unterscheidbar würden. Der neue 
Zustand brauchte sich nicht unbedingt durch eine veränderte Masse auszu- 
zeichnen, sondern könnte durch irgendeine, anderweitig gar nicht beobachtbare | 
Variable hervorgerufen werden. Im Prinzip müßte es möglich sein, experi- _ 
mentell die Frage zu entscheiden, ob gestreute Protonen in jeder Hinsicht mit — 


anderen übereinstimmen, wozu sich am besten solche Versuche eigneten, deren = 


Ergebnisse von der Symmetrie der Wellenfunktion abhängen. Die ‘amp 
liegt nahe, daß die Isomere über Zwischenzustände des p—p-Systems mit 
hoher Anregungsenergie gebildet werden. 

Wenn man annimmt, daß die Isomere den Spin !/, und ungerade Parität — 
haben, so wären bei Erhaltung des Gesamtdrehimpulses und der Parität 
Übergänge von ®P, zu 3S, möglich. Da der Zustand ®P, das dreifache Gewicht 
des ohne Isomerenhypothese allein erlaubten S-Zustandes 1S, hat, resultiert 
ein statistischer Faktor 4 im Streuquerschnitt, der es erlauben müßte, die 
Isotropie auf das Vorherrschen reiner S-Streuung zurückzuführen. : 

Nun darf man sich aber wohl kaum vorstellen, daB alle Teilchen, die sich im _ 
3P,-Zustand befinden, in den 3S,-Zustand überführt werden. Der Faktor im _ 
Wirkungsquerschnitt wird demnach irgendwo zwischen 1 und 4 liegen und 
es dürfte schwer sein, seine genaue Größe anzugeben. Aus dieser Betrachtung 


ergibt sich die weitere Schwierigkeit, daß man wahrscheinlich einen beträcht- = 
lichen Anteil von ®P,-Streuung in Kauf nehmen muß, wenn man für eine 
genügend große S-Streuung sorgen will. Weiter werden auch entsprechende 
Übergänge *F, > 3D,, 3H, > °G, usw. auftreten, deren Intensitäten kaum 


vorauszusagen sind. 
Wollte man nun versuchen, vom phänomenologischen Standpunkt aus ~ 
ein Modell zu bauen, das den Mechanismus der Anregung von Protoneniso- 


meren über Zwischenzustände erklärt und außerdem der eben aufgezählten — 


Schwierigkeiten Herr wird und die experimentellen Daten richtig voraussagt, 
so dürfte man dabei wahrscheinlich nicht ohne verschiedene, physikalisch 
mehr oder weniger plausible Zusatzhypothesen auskommen. Unter diesen 


Umständen erhebt sich die Frage, ob esnicht möglich sei, eine einzige zusätzliche 5 Er; 
Annahme zu finden, die keine detaillierten Modellvorstellungen benötigt = 


und somit zwar auf die Einzelheiten der Entstehung isomerer Protonen ver- 
zichtet, dafür aber realtiv einfache quantitative Rechnungen erlaubt. | 

Eine solehe Annahme könnte darin bestehen, daß man die einfallenden 
Protonen schon vor dem eigentlichen Streuakt, also schon bevor sie in die _ 


Reichweite der Kernkräfte der Zielprotonen kommen, als unterscheidbar von 


ihnen ansieht. Man könnte zur Begründung dieser Annahme etwa die Vor- 


stellung bilden, daß bei der Annäherung zweier Protonen aneinander mit e 


hoher relativer Geschwindigkeit ein Mechanismus abläuft, der die Ununter- Re 
scheidbarkeit der beiden Protonen aufhebt. ii 
Freilich ist die provisorische Natur dieser Annahme nicht zu übersehen. 


Trotzdem läßt sich aber der Ansicht vielleicht nicht jede Berechtigung ab- HS 


sprechen, daß die Hypothesen, die einer tiefergehenden Theorie der Protonen- 
lsomere zugrunde gelegt werden müßten, in ihrer Gesamtheit unter Umständen 
nicht wesentlich plausibler sein würden. 
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Bedenkt man einerseits, daß unsere Kenntnisse von den Kork 
noch recht lückenhaft sind und daher der Wert phänomenologischer Modelle, 
die ins einzelne gehen, nicht recht sicher ist; bedenkt man andererseits, daß 
die vorgeschlagene Zusatzhypothese einfach ist und übersichtliche Rechnungen 
erlaubt, so wird der Versuch, die p—p-Streuung mit Hilfe dieser einfachen 
Hypothese zu deuten, vielleicht nicht ohne Interesse erscheinen. 

Dieser Versuch wird hier unternommen. 

Da sehr starke Potentiale auftreten, erschien eine Art relativistischer Be- 
handlung angezeigt, die allerdings nicht konsequent sein kann. Da keine 
exakte relativistische Lösung des Zweiteilchenproblems existiert und auch die 
von Breit?’)?5) angegebenen Korrekturen für sehr starke Potentiale ihren 
Sinn verlieren, wird neben der gewöhnlichen die relativistische Schrödinger- 
Gleichung verwendet, die ja sicher bessere Ergebnisse als eine völlig nicht- 
relativistische Behandlung liefert. 


§ 6. Beziehungen für den differentiellen Wirkungsquerschnitt 


Unter der Voraussetzung, daß die einfallenden Protonen von den ruhenden 
zu unterscheiden sind, gilt für den differentiellen Streuquerschnitt pro Raum- 
winkel im Schwerpunktssystem ohne Berücksichtigung der Coulombkräfte 


(9) = (O) = o,, (9) + 0,, (a — 9), 


wo unter @ der Streuwinkel in diesem System zu verstehen ist; denn auch die 
Rückstoßteilchen werden registriert. 
Mit dem bekannten Ausdruck für o,, bei spinunabhängigen Kräften wird 


(1) 
-(22+ 1) (22 + 1) sin 6, sin 6, cos (6,—46,)P, (cos ©) Py (cos ©) 


mit k? = MB/2 


wo M = M, = M,, die Nukleonenmasse, B die Energie der einfallenden 
Protonen im Laborsystem und P, Legendresche Polynome bedeuten. 
Bis einschlieBlich der D-Phasen lautet (1) 


(0) = g[sin? dy + 9sin? d, - P? (cos O) + 25sin? 6, - P3 (cos ©) 

(2) 
j 

Schreibt man die asymptotische Entwicklung der Wellenfunktion im Schwer- 
punktssystem in der Form : ER 


+ 10 sin ö, sin 6, cos — P, (cos ©) 


einfallende Welle + (@) “ 


i so wird der differentielle Wirkungsquerschnitt in unserem Falle 2°) 
2?) G. Breit, Physic. Rev. 51, 248, 778 (1937). 
28) G. Breit, Physic. Rev. 53, 153 (1938). 
*®) Bei nichtunterscheidbaren Protonen steht dafür |f(©) + f (a — 9)]?. 
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Wir berücksichtigen jetzt die Coulom bwechselwirkung mit und erhalten durch 
Einsetzen des bekannten f (©), das sich bei Verwendung der nichtrelativisti- 
schen Schrödinger-Gleichung ergibt (siehe z. B. Schiff®®), in (3): 


1 


4 
(21+ 1) sin 6,” P, (cos 8) 


l 


Fr 5 COS log sin? 0/2 + 2 (y, — yo) + 
+ 55 cos (n' log cos? 0/2 + 2 &°)] 
1)* sin? P? (cos ©) 


+ 21+ + 1) (c08@) Py (cos) 


. sin 6," sin 6,” cos [2 (vy, —vv) + 6? — 
mit = e? (M/2 h? B)t = (ehe) (M y, = arg (1+ 1+ in). 
Bis einschließlich 1 = 2 lautet dieser Ausdruck 


2 
o (0) =@' (0) + > [G;. (9) cos + G 0) sin sin + (9) (5) 
t=0 


1 +1) P, (cos ©) 2079 [n’ log sin? 0/2 + 2 (y, — yo)] 


+ cos [n’ log cos? @/2 + 2 ( ) \ 
Gi, (9) = BE T (9 1+ 1) P, {cos ©) amo -sin [n’ log sin? O/2 + 2 (y, — y9)] 


1 
+ sin log cos? 8/2 + 2 
— g (9) aus Formel (2) mit 6;” statt 4,. 


In dem Ausdruck fir o’'® (©) ist (vy, — yo) gegen (d)) — 6) vernach- 
lässigt worden, was bei hohen a erlaubt ist, da in diesem Fall Ye— Vo 
sehr klein wird. Die Phasen ö{” rühren im wesentlichen von den Kern-, 
zum Teil aber auch von den Coulombkräften her. Vernachlässigt man jedoch 
die Coulombsche Wechselwirkung innerhalb der Reichweite der Kernkräfte 
(also für r < a), was bei sehr starken Kernkräften den wirklichen Verhältnissen 
sicher sehr nahe kommt, und bedenkt man, daß dann bei hohen Energien die 
Coulomb-Energie e2/r (r >a) nur eine kleine Störung darstellt, so erkennt 
man, daß bei hohen Energien ö/” ~ 6, gesetzt werden darf, wobei ö, die Phasen- 
verschiebungen sind, die bei Abwesenheit der Coulombkräfte entstehen. 


%) L. I. Schiff, Quantum Mechanics, New York, McGraw-Hill Book Comp., 1949. 
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Bei Bien — und nicht zu kleinen Winkeln © beläuft sich re Ein- 
fluß der Coulombkräfte nur auf eine Korrektur und o’ (©) ist im wesent- 
lichen gegeben durch a’‘*) (9). 

Um die Massenveränderlichkeit bei hohen Energien zu berücksichtigen, 
verwenden wir die relativistische Schrödinger-Gleichung, deren Radialteil 
für die Streuung von Teilchen der Masse M, zwischen denen nur Coulomb- 
kräfte wirken, lautet 


+[#- r 
wobei 
=n (1+ B/2 Mc*)t 
ist. 

In der entsprechenden nichtrelativistischen Gleichung, aus der die Cou- 
lom banteile in den Formeln (4) und (5) hervorgehen, fehlt der letzte Term 
und steht 7’ statt n””. 

Wenn vorausgesetzt wird (s. o.), daß nur für r >a Coulombkräfte vor- 
handen sind, ist der letzte Term überall sehr klein und kann weggelassen werden. 
Daraus folgt sofort, daß der mit der relativistischen Schrödinger-Gleichung 
errechnete Streuquerschnitt o’’ (©) und G’’ (©) usw. genau dieselbe Form wie 
die entsprechenden einmal gestrichenen Größen haben, nur daß überall 7" 
statt n’ steht. Aus den gleichen Gründen wie bei der nichtrelativistischen Be- 


handlung darf man auch hier bei hohen Energien 6," durch 6, ersetzen. 

Die in den Ausdrücken für o, o’ und o”’ auftretenden Phasenverschiebungen 
6, sind für ein Kastenpotential aus den Formeln (1.1) bis (1.5) einer vorher- 
gehenden Arbeit?!) zu entnehmen. Dabei ist für o und o’ B = ß’, für o” da- 
gegen ß = 6” zu setzen. 


$ 7. Bestimmung geeigneter Potentiale aus experimentellen Daten 
und zugehörige Streuquerschnitte 32) 


Neuere Messungen von Chamberlain und anderen®) ergeben | ™ den 
differentiellen Streuquerschnitt im Schwerpunktssystem Er 


o (©) = 3,5 + 0,4 mbarn REN 


für alle (nicht zu kleinen) Winkel © und für Teilchenenergien PER Lab 
system) von 150 bis 350 MeV. 
Nach $6 gilt für reine S-Streuung Te 


Als geeignet sollen solche Kastenpotentiale (Reichweite a, Tiefe V,) aus- 
gewählt werden, bei denen das zugehörige ö, nach (1) für B = 150, 250 und 
350 MeV ein o’”’(& zwischen 3 und 4 mbarn liefert. Dieses Verfahren ist sinn- 
reich, wenn die Phasen 6, (1 > 0) klein sind. Die gesuchten Potentiale lassen 
sich graphisch bestimmen: Zu jeder der drei Energien werden zwei Kurven 

31) D. Geißler, Ann. Physik 18, 114 (1956); im folgenden als A zitiert. 

32) Im folgenden wird stets die relativistische Schrédinger-Gleichung zugrunde 
gelegt. Entsprechende Untersuchungen wurden mit Verwendung der nichtrelativisti- 
schen Gleichung angestellt. Das Ergebnis ist am Schluß der Arbeit erwähnt. 
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Y,= V, (a) mit dy = const gezeichnet, und zwar eine mit dem do, das zu o’'(B) 
= 3 mbarn, die andere mit dem 6,, das zu o’’(X) = 4mbarn gehört. Die einzelnen 
Kurvenpunkte werden mit den Beziehungen aus A berechnet. 

Das Ergebnis ist in Abb. 1 zu sehen. Alle Potentiale im schraffierten 
Gebiet sind brauchbar. 

Bei höheren Werten von a gibt es noch weitere Potentiale, die nach For- 
mel (1) einen Wirkungsquerschnitt innerhalb der experimentellen Fehler- 
grenzen liefern, die aber deswegen nicht betrachtet zu werden brauchen, weil 
sie zu große Werte von 6, und ö, ergeben, so daß der gesamte Streuquerschnitt 
deutlich anisotrop würde. 

Von den geeigneten Parameterwerten wurden die folgenden näher unter- 


Vo 


1687 
1713 
1556 
1582 
1608 
1474 95 
1502 26 
1392 
4 
1338 30 
1366 
1394 
310 | 8 
1338 
4969 
V, in MeV. Ale 
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Sie sind in Abb. 2 verzeichnet. a a 
Für diese Potentiale wurden für die Energien B = 100, 150, 250, 3500 
und 435 MeV die Phasenverschiebungen ö, (l = 0,1,2) nach den Formen 
_1-Z,-tga 

a 


cotg ß JE 
Bei 


a? (M = +), 
= 2Mc*(1+ B/2 M c*), AN 
berechnet, deren Gültigkeit aus den Anschlußbedingungen (1.3) bis (1.5) aus A 
folgt. Dabei mußte mit einiger numerischer Genauigkeit gearbeitet werden, 
da Z, tg « und Z, tg « meist in der Nähe von 1 liegen. ir 


Mit den so gewonnenen Werten für die ö, wurde mit Hilfe der Beziehungen 
des vorigen Paragraphen der gesamte Streuquerschnitt o’’ (©) einschließlich _ 


Z, = 


+ 


Ein- 
ent- 
gen, 
teil 
m b- 
sucht: 
| 
om Nr. = 1 Nr. | a Vo 
4 | 0,38 1280 eee 
vor- 0,38 1309 
0,38 1326 
0,48 1224 
Ln" 0,50 1281 
Be- | 0,52 1222 
0,52 1251 
0,52 | 1269 
gen 0,54 1204 
her- 0,54 1222 
0,56 1194 
0,56 1223 
19 
a in 10-13 cm, 
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and 
nn- 
sen 
ven 
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der fiir die fiinf verschiedenen Energien bestimmt, zwar 
eV für die folgenden drei Krafttypen: 
‘ a) Gewöhnliche Kräfte, d.h. es treten alle Zustände, in der hier ver- 
; of wendeten Näherung also S-, P- und D-Zustände, auf. 
ER b) Die potentielle Energie enthält einen Austauschoperator vom Serber- 
Typ 1 
ER : worin unter PM der Majorana-Operator zu verstehen ist. Bei der gewöhnlichen 
4 en ‘ Schrédinger-Gleichung, wo die potentielle Energie linear auftritt, werden 
durch den Serber-Operator alle Zustände mit ungerader Bahndrehimpuls- 
quantenzahl / (P, F,...) ausgeschieden. Wegen 
er 4 gilt dasselbe aber auch bei der relativistischen Schrödinger-Gleichung, die 
in der potentiellen Energie quadratisch ist. In unserem Fall kommen dann 

' nur S- und D-Zustände vor. 

c) Die potentielle Energie enthält einen Integraloperator, der nur Wechsel- 
wirkung in S-Zuständen zuläßt und dessen Quadrat in der gleichen Weise 
wirkt wie der Operator selbst). 

Bei Typ a) scheiden wegen zu großer P-Phase und daraus folgender starker 

_ Anisotropie bei höheren Energien Potentiale mit a > 0,38 - 10-13 cm (Nr. 23 

und höher) aus. Bei Typ b) liegt der Wirkungsquerschnitt bei 90° für alle 

5 Potentiale mit a < 0,48 - 10-13 em bei 435 MeV wesentlich zu tief, so daß hier 

die Nummern 4 bis 22 ausgeschlossen werden. Ein Teil der Ergebnisse ist in 
den Abb. 3 bis 12 dargestellt. 


§ 8. Vergleich mit dem Experiment 


In der folgenden Zusammenstellung der experimentellen Ergebnisse sind 

= die Fehler zum Teil geschätzt, da in manchen Arbeiten nur der von der Zähl- 

2 = statistik herrührende Fehleranteil angegeben wird, und ab und zu ein wenig 

nach oben abgerundet, was bei dem ohnehin mehr qualitativen Charakter des 

Vergleichs nichts schaden wird. Bezüglich der Winkelabhängigkeit sind nur 

ER bei den 435 bzw. 429 MeV-Messungen die Ergebnisse für die verschiedenen 

= Winkel gesondert aufgeführt und in die Abb. 3, 8 und 9 eingetragen. Bei den 

anderen Energien sind sich alle Autoren bezüglich der Isotropie bei nicht zu 
kleinen Winkeln einig. o ist stets in mbarn angegeben. 


IR 435 MeV: Mott und andere?) 

0 o 
25,2° 4,49 + 0,61 65,3° 3,59 + 0,32 
36,4° 4,16 + 0,39 90° 3,39 + 0,30 

429 MeV: Marshall und andere”) 

o 
28° 2,86 + 0,43 65° 3,11 + 0,38 
43° 3,18 + 0,40 80° 3,51 + 0,44 
54° 2,82 + 0,39 90° 3,42 + 0,34 


; 33) J. M. Blatt u. V. F. WeiBkopf, Theoretical Nuclear Physics, New York, John 
Wiley u. Sons, 1952. 
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350 bis 150 MeV: Chamberlain und andere®) 0... 

= 3,5 + 0,5 

40MeV: Oxley und andere’), Towler®) = 
° o = 4,81 + 0,38 (IV) 


147 MeV: Cassels und andere) 1%) 
105 MeV: Birge und andere?) 
o = 5,4+ 1,1. uly (VI) 
Außerdem liegen ältere Messungen von Chamberlain und anderen>)%) 
bei 120, 164, 250 und 345 MeV vor (VII), die ebenso wie die neueren (III) mit 
den Ergebnissen (IV) und (V) nicht vereinbar sind, während sie gut zu (II) 
ee Für a (90°) werden folgende Ergebnisse verwendet, die in den Abb. 10 bis 
12 eingezeichnet sind: 
429 MeV: Marshall und andere 1) o (90°) = 3,42 + 0,34 (VII) 
271 MeV: Marshall und andere”) o (90°) = 3,67 + 0,5 (IX) 
144 MeV: Marshall und andere?) o (90°) = 3,21 + 0,4 (X) 
350 bis 150 MeV: Chamberlain und andere®) o (90°) = 3,5 + 0,5 (XI) 
345 MeV: Chamberlain und andere®) o (90°) = 3,7 + 0,4 (XII) 
247 MeV: Chamberlain und andere) o (90°) = 3,4 + 0,5 (XII) 
164 MeV: Chamberlain und andere5) o (90°) = 3,8 + 0,6 (XIV) 
119 MeV: Chamberlain und andere®) o (90°) = 3,95 + 0,4 (XV) 
105 MeV: Birge und andere?) o (90°) -5,4+ 1,1. (XVI) 


Typ a) (Abb. 3 und 10): 

Nach den Messungen I liegen bei 90° alle sich aus den betrachteten Po- 
tentialen (Nr. 4 bis 22) ergebenden Streuquerschnitte zu tief, während bei 
den kleineren Winkeln nur die Nummern 19 bis 22 (© = 65°) und 17 und 18 
(9 = 36°, 25°) passen, die sich aber durch besonders große Anisotropie 
auszeichnen. Letzteres gälte erst recht für Nummern größer als 22 (a > 
0,38 - 10-13 em). 

Ähnlich ist die Situation beim Vergleich mit den Ergebnissen II, wo 
ebenfalls bei 90° und auch bei 80° alle Nummern deutlich zu tief liegen. Bei 
allen anderen Meßwinkeln passen Nr. 6, 7, 9 bis 13, außerdem bei 65° Nr. 14 
bis 18, bei 54° Nr. 4, 5, 8, 14 bis 16, bei 43° Nr. 14 bis 16 und bei 28° noch 
Nr. 4, 5 und 8. 

Mit den Messungen III stimmen bei 350 und 250 MeV die Nummern 4 
bis 18 bei allen Winkeln überein, bei 150 MeV dagegen liegen alle zu hoch. 
Andererseits zeigt sich beim Vergleich mit den älteren Messungen VII, daß 
die theoretischen Werte bei 350 MeV im allgemeinen etwas zu tief liegen, 


_™) E. Segré (Editor): Experimental Nuclear Physics, Vol.I, New York, John 
Wiley u. Sons, 1953. 


’er- 
en 
len 
lie 
nn 
el- 
se 

er 
3 
le 
er 
in 
d 
|- 
g 
3 

n # 
u Br: 


140 Annalen der Physik. 6. Folge. Band 18. 1956 


während sie alle bei 250 und 150 MeV bei allen Winkeln innerhalb 
mentellen Fehler bleiben. ; 

Verglichen mit den Ergebnissen IV und V liegen sämtliche Nummern fast 
stets zu tief, was nicht verwunderlich ist, da die Potentiale im Hinblick auf 
Übereinstimmung mit den bedeutend kleineren Werten III ausgesucht worden 
sind. 

Zu VI passen alle Nummern bei allen Winkeln. 

Betrachtet man o (90°) als Funktion der Energie, so zeigt sich, daB die 
Potentiale 4 bis 22 mit den Messungen IX, XIII, XIV, und XVI, die Poten: 
tiale 8, 16, 18, 21 und 22 auch mit XV übereinstimmen, während zu VIII, X 
und XII überhaupt keine der Nummern paßt. Von 200 MeV an sind alk 
Potentiale in Einklang mit dem Experiment XI, wobei die Nummern 13, 16, 
18, 21 und 22 besonders günstig liegen. 

Obwohl der flache Verlauf von o (90°) in Abhängigkeit von der Energie, 
der von Chamberlain und anderen®)*), Marshall und anderen”) und von 
Cassels und anderen‘) aus den Experimenten abgelesen wurde, nicht sehr 
gut wiedergegeben wird, bleiben doch die Wirkungsquerschnitte vieler Po- 
tentiale innerhalb der Fehlergrenzen der meisten Messungen. 

Insgesamt gesehen ist etwa Nr. 16 mit den Parametern . 


a = 0,34-10- em, V,=1894MeV 


ein relativ giinstiges Potential. 

Typ b) (Abb. 4 bis 8 und 11). 

Zu I passen nur bei 90° die Nummern 24 und 27 bis 34, zu den anderen 
Winkeln dagegen keins der Potentiale. 

Abgesehen von dem Wert bei 43°, wo das Meßergebnis etwas abseits liegt, 
stimmen mit II bei allen anderen Winkeln die Nummern 24, 25 und 27 bis 30, 
bei 54 bis 90° außerdem die Nummern 31 bis 34 und bei 28° zusätzlich noch 
die Nummern 23 und 26 überein. 

Im Gegensatz zu Typ a) sind hier die meisten der herausgegriffenen Po- 
tentiale in Einklang mit den Experimenten bei B > 420 MeV und großen 
Streuwinkeln zu bringen. Die leichte Abweichung von der Isotropie, die sich 
nach Marshall und anderen”) bei 429 und 345 MeV in Form eines Anstiegs 
von o mit wachsendem © beobachten läßt, wird richtig wiedergegeben. Alle 
Nummern kleiner als 23 liefern zu kleine Wirkungsquerschnitte bei 90° und 
435 MeV, während alle Potentiale mit a > 0,56 - 10-13 em wegen zu großer 
Anisotropie bei dieser Energie ausgeschlossen werden. 

Zu dem Experiment III passen bei 350 MeV die Nummern 23, 24, 25 und 
27 bei allen Winkeln >30° und ab 50° sämtliche Nummern, bei 250 MeV 
alle Potentiale für alle Winkel, während bei 150 MeV wieder alle höher als 
das Meßergebnis liegen. 

Für die übrigen Energien und Messungen (IV, V, VI, VII) gilt bezüglich 
der Abhängigkeit vom Streuwinkel qualitativ hier genau dasselbe für die 
Potentiale 23 bis 34, was bei Typ a) für die Potentiale 4 bis 22 zutrifft. 

Was die Abhängigkeit des Wirkungsquerschnitts bei 90° von der Energie 
anbetrifft, so wird jetzt der experimentell beobachtete flache Verlauf besser 
als bei Typ a) wiedergegeben. Im einzelnen passen zu Messung VIII die 
Nummern 24, 25 und 27 bis 34, zu IX Nr. 23 bis 34, zu XI Nr. 23 bis 34 ober- 
halb von 165 bis 220 MeV, zu XII Nr. 23 bis 34, zu XIII Nr. 23 bis 32 und 34, 
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Abb. 1. Auswahl geeigneter Potentiale. Abb. 2. Parameterwerte und Nummern 
Alle Kästen im schraffierten Gebiet liefern der untersuchten Kastenpotentiale 
bei reiner S-Streuung Streuquerschnitte 
zwischen 3 und 4 mbarn für Protonen- 

energien zwischen 150 und 350 MeV 
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60° 
Abb. 3. Differentieller Wirkungsquer- 
schnitt mit Coulombkorrektur für 43 
Kräfte vom Typ a) bei B = 435 MeV. ” 
Die arabischen Ziffern an den Kurven Abb. 4. Differentieller Wirkungsquer- 
geben die Nr. des zugrunde gelegten schnitt mit Coulom bkorrektur für Kräfte 
Potentials an (vgl. Abb. 2), die römischen vom Typ) bei B = 100 MeV. Bezeich- 
gehören zu den Meßwerten (vgl. § 8) nungen wie bei Abb. 3 
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zu XIV Nr. 23 bis 34, zu XV 
Nr. 25, 26, 28, 29, 32 und zu 
XVI Nr. 23 bis 34. Man sieht, 
daß hier die Ü bereinstimmung 
besser als bei Typ a) ist. 

Ein besonders günstige 
Potential ist Nr. 29 mit den 
Abmessungen 


a= 0,52 - 10-13 em, 
V, = 1269 MeV. 
90° Type) (Abb. 9 und 12). 


. . d . 
Abb. 5. Differentieller Wirkungsquerschnitt mit Cou- ti 
 lombkorrektur für Kräfte vom Typb) bei B=150 Ale 0 DIS ©, 19, 18 29, 

MeV. Bezeichnungen wie bei Abb. 3 27, 31, 32 und 34 behandelt 


worden. Im großen und gan 
o’(0) zen sind die Verhältnisse ähn- 
lich wie bei Typ a). 
Sämtliche Potentiale geben 
Ergebnisse, die nach den Mes- 
sungen I zu tief liegen. Das- 
selbe gilt für die Messungen II 
bei den Winkeln 80 bis 90°, 
während für die kleineren @ 
fast alle Nummern zu passen- 


T 
Abb. 6, Differentieller Wirkungsquerschnitt mit Cou. W ten führen. 
- lombkorrektur für Kräfte vom Typ b) bei B = 250 Hirte 
MeV. Bezeichnungen wie bei Abb. 3 ese 
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_ Abb. 7. Differentieller Wirkungsquerschnitt | Abb. 8. Differentieller Wirkungsquerschnitt 
mit Coulombkorrektur für Kräfte vom mit Coulombkorrektur fiir Krafte vom 
_ Typd) bei B = 350 MeV. Bezeichnungen wie Typ 6) bei B = 435 MeV. Bezeichnungen wie 

bei Abb. 3 bei Abb. 3 
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Mit den Experimenten III 
stimmen fiir 350 und 250 MeV 
alle Nummern fiir 30 bis 90° 
überein, bei 150 MeV da- 
gegen liegen die theoretischen 
Werte sämtlich zu hoch. 

Für die anderen Energien 
und Messungen (IV bis VII) 
gilt für alle Potentiale das- 
selbe wie bei den Typen a) 
und b). 

Beim Vergleich der theo- 
retischen Werte von o (90°) 
in Abhängigkeit von der 
Energie B der einfallenden 
Protonen mit den experimen- 
tellen zeigt sich, daß alle 
Potentiale zu den Messun- 
gen XI (von etwa 170 bis 
190 MeV an), XIII, XIV 
und XVI, die Potentiale 21, 
25, 28 und 32 zu der Messung 
XV passen. Der Verlauf ist 
weniger flach als bei Typ b). 


$9. Schlußfolgerungen 

Wenn man das Problem 
der p—p-Streuung bei Ener- 
gien von 100 bis 435MeV mit 
den Annehmen 

l.es existieren Protonen- 
isomere ; 

2. die einfallenden Pro- 
tonen sind schon vor dem 
Streuakt von den ruhenden 
Protonen verschieden; 

3. die  relativistische 
Schrédinger-Gleichung ist 
dem Problem angemessen: 

4. F- und höhere Phasen 
können vernachlässigt wer- 
den; 

5. das Kernpotential hat 
Kastenform 
und für die drei Typen von 
Kräften 

a) gewöhnliche Kräfte; 

b) Kräfte, die nur Wech- 
selwirkungen in Zuständen 
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Abb. 9. Differentieller Wirkungsquerschnitt mit Cou - 

lombkorrektur für Kräfte vom Typc) bei B = 435 
MeV. Bezeichnungen wie bei Abb. 3 
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Abb. 10. Differentieller Wirkungsquerschnitt mit 

Coulombkorrektur bei © = 90° als Funktion der 

Energie für Kräfte vom Typ a). Bezeichnungen wie 
bei Abb. 3 
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bei Abb. 3 
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mit Drehimpulsquantenzahl erlauben (Austauschkräfte vom 
Typ); 

e) Kräfte, die nur Wechselwirkungen” in S-Zuständen erlauben 
behandelt, dann erhält man folgende Ergebnisse: 

Um die Experimente richtig deuten zu können, muß man für alle drei 
Typen von Kräften sehr schmale und sehr tiefe Kastenpotentiale verwenden, 
und zwar mit den Reichweiten 


j —13 
0,26 bis 0,56 10 cm alle: 


1700 bis 1200 MeV. 


und den Tiefen 

Die Typen a) und c) widersprechen dem experimentellen Befund insofern, 
als sie beide bei Energien von mehr als 420 MeV zu kleine Werte für den 
Wirkungsquerschnitt bei Streuwinkeln von 80 bis 90° (im Schwerpunkts- 
system) liefern und der Ver- 
lauf von o (90°) in Abhän- 
gigkeit von der Energie nicht 
gerade sehr flach ist sowie 
auch zahlenmäßig mit einer 
ganzen Anzahl von Messun- 
gen nicht übereinstimmt. 
Der Typ e) hat den Vorzug 
größerer Isotropie gegenüber 
dem Typ a). 

Für Kräfte vom Typ b) 
liegen dagen die Verhält- 
nisse günstiger; denn hier 
können durch geeignete 
20+ N Potentiale die Werte von o 
bei Energien von mehr ak 
0 150 200 20 20 350 WO MeV 450 

Abb, 12, mit Sr Energie 

Energie für Kräfte vom Typ c). Bezeichnungen wie als bei den Typen a) und ¢), 

bei Abb. 3 und o (90°) stimmt außer- 

dem auch in der Größe mit 

den meisten Meßergebnissen überein. Als weiterer Vorzug ließe sich unter 

Umständen anführen, daß der eventuell aus den Experimenten ablesbare An- 

stieg von o mit wachsendem Streuwinkel bei Kräften vom Typ b), nicht aber 
bei'solehen der Typen a) und c) vorhanden ist. 

Eins der relativ günstigsten Potentiale vom Typ b) hat die Parameter- 
werte a = 0,52 - 10-13 em, V, = 1270 MeV. 

Legt man der Rechnung die nichtrelativistische Schrödinger-Gleichung 
zugrunde, so ändern sich die Verhältnisse nur darin, daß die hinsichtlich der 
experimentellen Ergebnisse brauchbaren Potentiale um mehrere hundert MeV 
größere Tiefen haben müssen (z. B. das obige Potential mit a = 0,52 - 10-% cm 
ein V, von 1780 MeV), während sonst alles im wesentlichen wie bei der relativi- 
stischen Rechnung bleibt. 
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Vielleicht ist es von einiger Bedeutung, daß sich auch bei der hier durch- 
geführten Art der Behandlung der p—p-Streuung bei hohen Energien heraus- 
stellt: 

1. Zwischen zwei Protonen existiert eine starke Wechselwirkung bei Ab- 
ständen von etwa 0,5 - 10-13 cm, was darauf hindeutet, daß schwerere Mesonen 
als das Pi-Meson mit für die Kräfte zwischen Nukleonen verantwortlich sind. 

2. Eine Austauschkraft vom Serber-Typ ist anderen Kräften überlegen. 

Anmerkung: In letzter Zeit wurden Ergebnisse von Smith und anderen*) 
und von Morris und anderen®) bekannt, die bei sehr hohen Protonenenergien 
(bis 1000 MeV) eine deutliche Bevorzugung der Vorwärtsstreuung zeigen. 


ber- 


Ich danke Herrn Prof. Kockel für die Anregung zu dieser Arbeit und für 
wertvolle Ratschläge bei ihrer Durchführung. 


%) L. W. Smith, A. W. McReynolds u. G. A. Snow, Physic. Rev. 98, 1167 (1955). 
s) T. W. Morris, J. D. Garrison, E. C. Fowler, W. B. Fowler, R. P. Shutt, 
A.M. Throndike u. W. L. Whittemore, Physik. Rev. 98, 1167 (1955). os 


Leipzig, Theoretisch-Physikalisches Institut der Universitat. 
Dr Bei der Redaktion eingegangen am 3. Dezember 1955. aid Rises. 
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falls. 
_ Die Elektronendichte und Geschwindigkeitsverteilung selbs 
bei der Pendelvervielfachung von Sekundärelektronen F 


tli 
Von K. Krebs und H. Meerbach 


ntla 
Mit 11 Abbildungen 


Die im stationären Zustand einer Pendelvervielfac '»hung zwischen den 
Elektroden bewegte Elektrizitätsmenge wird dadurch bestimmt, daß die von 
ihr an einer Elektrode erzeugte Temperatur gemessen wird. Dazu werden 
wohldefinierte Verhältnisse der Wärmeableitung geschaffen, so daß die Tem- 
peraturerhöhung direkt ein Maß für die Leistung des Konvektionsstromes 
der pendelnden Elektronen ist. Mit der im zweiten Teil ermittelten Energie- 
verteilung der Elektronen erhält man daraus die Pendelstromstärke zu einigen 
mA. — Mittels eines Systems von Käfigen, das an einer Elektrode angebracht 
ist, wird ferner nach einer Gegenfeldmethode die Energieverteilung der 
pendelnden Sekundärelektronen gemessen. Es ergibt sich eine Häufung der 
Energie um einen Wert herum, der zwar um durchschnittlich 20 % unter dem 
theoretischen Wert bleibt, jedoch den theoretisch erwarteten Gang mit der 
Feldlänge zeigt. Die Bevorzugung einer bestimmten Energie und ihr Gang 
mit der Feldlänge liefert eine Entscheidung für die seinerzeit entwickelte kine- 


matische Theorie der Pendelvervielfachung. 
Vee 


A. Allgemeines zur Pendelvervielfachung 


In einer früheren Arbeit!) wurde für die erstmalig von P. T. Farnsworth?) 
beschriebene Pendelvervielfachung von Sekundärelektronen eine kinematische 
Theorie angegeben und durch Messung der ‚„Einsatzspannung‘‘, d.h. der 
niedrigsten Hochfrequenzspannung, bei der sich unter fest gegebenen sonstigen 
Versuchsbedingungen eine Pendelvervielfachung ausbilden kann, experimentell 
bestätigt. In dem von Farnsworth ursprünglich ins Auge gefaßten An- 
wendungsgebiet, nämlich zur steuerbaren Verstärkung eines Elektronen- 
stromes, konnte sich die Pendelvervielfachung trotz der in neuerer Zeit er 
zielten Fortschritte?) offensichtlich nicht durchsetzen. Dagegen scheint sie 
eine allgemeine Bedeutung für die Probleme der Höchstfrequenztechnik zu 
gewinnen. So konnten in Laufzeitröhren auftretende Leistungsminderungen 
durch die Pendelvervielfachung erklärt werden*). Bei den von Smither und 


1) K. Krebs u. H. Meerbach, Ann. Physik 15, 189 (1965). ae or 2 

2) P. T. Farnsworth, J. Franklin Inst. 2, 411 (1934). FR ia 
“3 3) H. Beneking, Z. angew. Physik 4, 258 (1952), eet 
_ 4) K. Krebs, Z. angew. Physik 2, 400 (1950). ee en 
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Hoover’) kürzlich gemachten Beobachtungen scheint es sich, soweit aus dem 
kurzen Bericht zu ersehen, um ein ähnliches Problem zu handeln. 

Die von Posin 1946°) beschriebene Abhängigkeit der Zündspannung in 

füllten Zentimeter-Hohlleitungen wird von Kossel und Krebs’) eben- 
falls durch die Pendelvervielfachung erklärt, im Gegensatz zu der von Posin 
selbst gegebenen Erklärung durch Larmorpräzessionsresonanz. 

Franeis und v. Engel®) erklären mit Hilfe der Pendelvervielfachung den 
zeitlichen Verlauf des Anwachsens einer elektrodenlosen Hochfrequenz-Gas- 
entladung. Sie benutzen dazu eine von Gill und v. Engel?) entwickelte 
Theorie, auf die unten noch näher eingegangen wird und deren Voraussagen 
ineinigen Punkten im Widerspruch zu den experimentellen Ergebnissen dieser 
Arbeit stehen. 

Ebenfalls versuchen Hatch und Williams?!) die Ergebnisse ihrer Mes- 
sungen mit dieser Theorie zu erklären. Sie untersuchen die Zündspannung in 
Abhängigkeit von Frequenz und Feldlänge und stellen insbesondere auch eine 
Zündung fest, wenn die Hochfrequenzspannung, ausgehend von sehr großen 
Werten, langsam erniedrigt wird. Diese „obere Zündspannung‘“ liegt über 
der normalen Zündspannung und stellt sicher die in!) theoretisch geforderte 
und experimentell festgestellte ‚‚Aussetzspannung‘ dar. Hatch und Williams 
erwähnen auch kurz Versuche zur Messung der Geschwindigkeitsverteilung 
der Pendelelektronen, mit allerdings nur knappen Angaben über die Ergeb- 
nisse. Bei Feldlängen von 2cm ergaben sich Auftreffenergien um 100 Volt, 
was nicht im Widerspruch zu den Ergebnissen der vorliegenden Arbeit steht. 

Ferner besteht, wie verschiedentlich 1!) 12)3) ausgeführt, die Möglichkeit, 
die Pendelvervielfachung zur Spannungsstabilisierung in einem Höchstfrequenz- 
schwingkreis zu benutzen. Durch die einsetzende Vervielfachung wird dem 
Sender so viel Energie entzogen, daß die Spannung im Resonator nicht wesent- 
lich über die Einsatzspannung ansteigen kann. 

Schließlich sei noch erwähnt, daß es möglich sein muß!), Höchstfrequenz- 
spannungen in Hohlraumresonatoren zumindestens größenordnungsmäßig 
durch Beobachtung einer einsetzenden Pendelvervielfachung zu messen. 

Diese vielfache Beachtung der Pendelvervielfachung rechtfertigte die 
Fortsetzung der früher!) begonnenen Untersuchungen. Und zwar soll in der 
vorliegenden Arbeit die pendelnde Elektrizitätsmenge und im Zusammen- 
hang damit die Energieverteilung der Pendelelektronen bestimmt werden. 


B. Messung der Leistung des Pendelstromes 
a) Prinzip und Voraussetzungen des Meßverfahrens 


Zur Erläuterung des Meßverfahrens sei nochmals kurz an das Prinzip 
der Pendelvervielfachung erinnert. Die Elektronen pendeln zwischen zwei 
das Hochfrequenzfeld begrenzenden Elektroden derart, daß die Flugzeit von 


5) R.K. Smither u. C. W. Hoover, Physic. Rev. 98, 1149 (1955). 
*)D. Q. Posin, Physic. Rev. 69, 541 (1946). 

?) F. Kossel u. K. Krebs, Z. Physik 189, 189 (1954). 

‘) G. Francis u. A. v. Engel, Philos. Trans. Roy. Soc. (London) 246, 143 (1953). 
*) B. Gill u. A. v. Engel, Proc. Roy. Soc. (London) 192, 446 (1948). 

*) A. J. Hatch u. H. B. Williams, J. appl. Physics 25, 417 (1954). 

4) Patentschrift Nr. 882422 des Deutschen Patentamtes. = = 
M.H. Greenblatt, Physic. Rev. 75, 1824 (1949), 
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einer Elektrode zur anderen gerade gleich !/, Periode ist. Einzelheiten, sowohl 
der Theorie wie experimentelle, sind in oben zitierter Arbeit!) dargelegt, 
deren Ergebnisse hier im wesentlichen als bekannt vorausgesetzt werden müssen, 

Die pendelnde Elektronenwolke stellt einen hochfrequenten Konvek- 
tionsstrom dar. Wird seine Leistung — aus der Erwärmung der Elektroden — 
ermittelt, und außerdem die Geschwindigkeitsverteilung der Elektronen 
(im Augenblick des Aufpralls auf die Elektrode) bestimmt, so ist damit auch 
die den Konvektionsstrom bildende Ladungsmenge bekannt. Der erste 
Schritt besteht also in der Messung der Leistung des ,, Pendelstromes‘‘, worunter 
im folgenden die Elektrizitätsmenge in Coulomb verstanden werden soll, 
die pro Sekunde auf eine der Prallelektroden auftrifft. (Die Division durch 
die Frequenz ergibt dann aus dem Pendelstrom die wirklich zwischen den 
Elektroden pendelnde Elektrizitätsmenge.) 


Es ist natürlich klar, daß eine solche Messung nur für die speziellen Ver- 
hältnisse unserer Apparatur (Feldlänge 4 cm, Elektrodendurchmesser 2,6 cm, 
Material Al, Frequenz 80 MHz) Bedeutung haben kann und nur bedingt zu 
verallgemeinern ist. Dennoch halten wir ihre Mitteilung für gerechtfertigt, da 
eine derartige quantitative Bestimmung unseres Wissens bisher noch nicht 
vorliegt. 

Eine erste Abschätzung der Größenordnung der zu erwartenden Ströme 
konnte schon aus den Ergebnissen der Einsatzspannungsmessung!) ge- 
wonnen werden. Der durch das Loch in einer der Prallelektroden fließende 
Auffängerstrom®%) hatte die Größenordnung 105 A. Da der Lochdureh- 
messer etwa 1/,, des Elektrodendurchmessers betrug, konnte mit Pendel- 
strömen von einigen mA gerechnet werden. 

Der Konvektionspendelstrom setzt sich in den Zuführungen zu den Elek- 
troden als hochfrequenter Leitungsstrom fort und könnte dort gemessen 
werden. In den Zuführungen fließt aber auch der Leitungsstrom, der die 
Fortsetzung des zwischen den Prallelektroden fließenden Verschiebungs- 
stromes darstellt. Eine Trennung der sich überlagernden Teilströme ist nieht 
möglich. 

Der Verschiebungsstrom läßt sich aus den Abmessungen der Pendelstrecke und der 
Frequenz überschlägig berechnen. Mit den oben angegebenen Daten berechnet sich die 
Kapazität der Pendelstrecke zu 1,2 - 10-13 Farad und der Hochfrequenzwiderstand zu 
16,5k2. Da die angelegte Hochfrequenzspannung bei dieser Feldlänge etwa 800 Volt 
beträgt, liegt der Verschiebungsstrom mit rund 50 mA um eine Größenordnung über 
dem zu erwartenden Pendelstrom. Es ist also nicht mit hinreichender Genauigkeit mög- 
lich, den Konvektionsstrom durch eine Differenzmessung aus dem: Leitungsstrom zu be- 
stimmen. 

Anders liegen die Verhältnisse, wenn man die Leistung des Pendel- und 
des. Verschiebungsstromes an den Elektroden betrachtet. Ohmsche und 
dielektrische Verluste und somit die durch den Verschiebungsstrom hervor- 
gerufene Erwärmung der Elektroden sind gering. Die Energie der Kon 
vektionselektronen wird aber restlos an den Prallelektroden in Wärme umge- 
setzt. Die zu messende Leistung des Pendelstromes ist groß gegen die Lei- 
stung des Verschiebungsstromes und kann ohne Bedenken aus einer Differenz- 
messung ermittelt werden. 


13) Siehe Abb. 8 in!). Dieser Auffänger arbeitete, im Gegensatz zu dem hier in Teil 
beschriebenen, noch ohne Gegenfeld. 
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Der Meßvorgang und seine gedanklichen Voraussetzungen seien an Hand 
der Abb. 1 erläutert. Die eine Prallelektrode ist ein homogener Körper mit 
der Wärmekapazität C und wird über den Plattenträger T (Querschnitt q) 
mit einem gekühlten Körper K ver- 


bunden, der schematisch als großes Y 
Wärmereservoir konstanter Tempe- SN Yj; 
ratur betrachtet werden kann. Fer- Z, 


ki 


ner kann angenommen werden, daß N 9 Z 

die Temperatur in der Prallelek- 

trode örtlich nicht wesentlich ver- ee / YY, 
schieden ist. Da sich das Ganze im ii 2 
Vakuum befindet, eine  Abb.1. Prallelektrode zur Leistungsmessung, 
Wärmeableitung nur über den 

Plattenträger. In diesem herrscht 

dann ein linearer Temperaturabfall mit einem Temperaturgradienten 
(#+—®m)/l. Die Wärmeabstrahlung ist klein gegenüber der Wärmeableitung, da 
die Temperaturen relativ niedrig sind und das Elektrodenmaterial Aluminium 
nur ein geringes Emissionsvermögen hat. Außerdem geht die Abstrahlung in 
dem betrachteten Temperaturbereich noch in ausreichender Näherung linear 
mit der Temperaturdifferenz'), Wenn man von der Strahlung zunächst 
einmal absieht, so ergibt sich aus der bekannten Wärmeleitungsgleichung 


4 6; — m. (1) 


Dabei bedeutet k die Wärmeleitzahl, L die in der Prallelektrode durch Elek- 
tronenaufprall und L’ die störende, durch direkte Hochfrequenzwirkung er- 
zeugte Wärmeleistung. 

Diese Gleichung läßt zwei verschiedene Möglichkeiten zur Messung der 
Leistung L erkennen. 

1. Im stationären Zustand verschwinden die Ableitungen nach der Zeit 
und es bleibt: 

, keg k-q 
L+ L' (01 — On) = Ay — On) (2) 

Unter der oben begründeten Voraussetzung, daß die abgestrahlte Wärme- 
menge in erster Näherung auch der Temperaturdifferenz proportional ist, 
kann jetzt die Strahlung auch wieder in die Rechnung mit einbezogen werden. 
Mit einem experimentell zu bestimmenden Wärmeabstrahlungsfaktor A, 
folgt dann ganz analog: 


L+L' =A — dy) mit A= Ao + A, (3a) 


baw. bei Abschalten der Pendelvervielfachung und Fortdauer der Hoch- 
frequenzspannung, wenn 9; die sich dann einstellende Temperatur der Elek- 


trode I bedeutet, 

4 Ric 
L=4(# —Bı). (8e) 
4) Das ergibt sich, wenn in der Strahlungsgleichung S=oF (Tt — T$) die abso- 


luten (7) durch die Celsiusgrade (#) ersetzt (7 = 273 + #) und die Potenzen von @ ge- 
sammelt werden. Bis 100° C überwiegt dann das in ® lineare Glied. 


aus (3a) und (3b) folgt 
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Danach besteht die Leistungsmessung in der Messung einer Temperatur- 
differenz. Der Ausdruck 4, = k-q/l wäre aus den Materialkonstanten des 
Plattenträgers zu berechnen. 

2. Auch der Temperaturabfall nach dem Abschalten der Pendelverviel- 
fachung (nichtstationärer Zustand) gibt die Möglichkeit, die Pendelleistung L 
zu bestimmen. Im Augenblick des Abschaltens verschwindet die einströ- 
mende Wärmemenge / und die Temperatur der Prallelektrode sinkt, da ja 
weiterhin Wärme abgeführt wird. Es ist dann nach Gl. (1) und (3b) unter 
Berücksichtigung der Strahlung 


8, 
Cy = (1 — On) = (1 A). (4) 


Da sich nun # nicht sprunghaft ändern kann, müssen fiir den Abschalt- 
zeitpunkt sowohl (3a) als auch (4) gelten. Daraus folgt durch Vergleich: 
08 


0) 


Das heiBt, die Leistung ist bestimmt durch die Tangente an die Temperatur- 
abfallkurve im Augenblick t, des Abschaltens und die Wärmekapazität (, 
die in diesem Falle aus den Materialkonstanten der Prallelektrode zu berechnen 
ware. 

Diese Überlegungen zeigen grundsätzlich die Berechtigung, zur Leistungsmessung 
einen der genannten Wege einzuschlagen. Ihre praktische Anwendung aber empfiehlt 
sich nicht, da es weder möglich ist, A, und A, mit genügender Genauigkeit zu bestimmen, 
noch die Wärmekapazität C der Prallelektrode ohne weiteres in die Gl. (5) eingesetzt 
werden kann. Bei der Berechnung der Wärmeleitzahl wären einmal die Wärmeüber- 
gangszahlen an den vorhandenen Verbindung- 
stellen zu berücksichtigen, und außerdem dürften 
die Tabellenangaben von Wärmeleitzahlen reiner 
Substanzen kaum exakt für unsere technischen 
Materialien gelten. Das C in der Gl. (5) würde 
nur dann mit der Wärmekapazität der Prall- 
elektrode identisch sein, wenn der Plattenträger 
selbst keine Wärmekapazität hätte. Streng ge- 
nommen müßte die Wärmekapazität des Platten- 
trägers berücksichtigt werden, was bei der zeitlich 
veränderlichen Temperaturverteilung in diesem die 
Rechnung sehr komplizieren würde. 

Um von allen Rechnungen und beson- 
deren Voraussetzungen unabhängig zu 
werden, wurde die Meßapparatur unmittel- 
bar durch einen Elektronengleichstrom ge- 
eicht. Wie in Abb. 2schematisch angedeutet, 
Abb. 2. Zur Eichung der Leistungs- wurde die Gegenplatte durch eine Glüh- 

messung, schematisch elektrode, die aus einem direkt geheizten 
Wolframdraht bestand, ersetzt. Bei der 

Spannung U von etwa 800 Volt floß der Sättigungsstrom, der durch Varia 
tion der Heizleistung reguliert werden konnte. Die obere Prallelektrode 
erwärmte sich dabei in der gleichen Weise, wie bei dem Auftreffen der 
Pendelelektronen, und die Temperatur im stationären Zustand sowie der 
Temperaturabfal! nach 'Abschalten des Stromes konnten in der gleichen, 
weiter unten beschriebenen Weise und mit den gleichen Einrichtungen ge 
messen werden, wie bei der Erwärmung durch Pendelelektronen. Auch hier 
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wurde durch Überlagerung eines Magnetfeldes dafür gesorgt, daß alle von 
der Glühelektrode ausgehenden Elektronen auf die Vorderseite der Prallelek- 
trode — und nicht auf den Platten- er 
träger T oder den Körper K — gelang- 
ten. 1504 °C 

Zur Messung der Temperaturdiffe- 
renz (0; — wurde ein Paar Thermo- 
elemente aus Einschmelzlegierung 501 100+ 
und Kupfer benutzt, das bis zu Tem- 6 and BR 
peraturen von 90° praktisch linear ER 


arbeitete, wie experimentell geprüft re 
wurde®). Dadurch war es möglich, 


direkt am Instrument Um abzulesen. 
Die durch die Strahlung der Glüh- 
kathode hervorgerufene Erwärmung sad 
wurde durch Differenzmessung elimi- Abb. 3. Eichkurve des Thermoelements 
niert, ganz analog wie die oben mit L’ 
bezeichnete reine Hochfrequenzerwärmung. Dies geschah, indem zunächst 
nur die Heizung eingestellt und so lange gewartet wurde, bis sich der Gleich- 
gewichtszustand hergestellt hatte (Temperatur der Elektrode = 9). Dann 
erst wurde die Anodenspannung 
angelegt und wieder bis zum 
Gleichgewichtszustand gewartet. __ 
Die aus Stromstärke I und © 
Spannung U berechnete Lei- 
stung L wurde der Temperatur- ____Yochfrequenzzu- 
differenz (dr — bzw. der führung 
Thermospannungsdifferenz 
— U'y,) zwischen den beiden Be. 
tionären Zuständen zugeordnet. 

Das Ergebnis einer 
Reihe solcher Messungen ist- in 
Abb. 3 dargestellt. Die Ab- . 
weichungen von der theoretisch 
zu erwartenden Geraden erklären 
sich daraus, daß jede 


Y 


mehr als eine Stunde dauerte 
und es nicht ganz leicht war, 
über diese Zeitspanne alle | ays 
suchsbedingungen, insbesondere 

die Kühlwassertemperatur ( 
die Beschreibung der 
apparatur weiter unten) völli in 


konstant 
15) Dadurch, daß nach dem Am ra | 
bau des Systems die Elektrode E, 
ein Wasserbad veränderlicher, ge- Abb.4. Experimentelle Anordnung zur Lei- 
messener Temperatur getaucht wurde. stungsmessung 
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Als Mittelwert über alle Messungen ergab sich: 
Leistung des Konvektionsstromes = 1,52 + 0,03 100 6 
2 b) Durchführung der Messungen 
Das in Abb. 1 dargestellte Schema der Wärmeableitung wurde praktisch 
durch einen Elektrodenaufbau verwirklicht, der in Abb. 4 dargestellt ist, 
Als Plattenträger wurde ein Kupferrohr von 0,5 mm Wandstärke und 4,0 mm 
äußerem Durchmesser verwendet. Die „heiße‘‘ Lötstelle I wurde in diesem 
Rohr dort angebracht, wo dieses nach unten durch einen Kupferstopfen ver- 
schlossen ist. Dieser letztere hat außerhalb des Rohres ein Gewinde, auf das 
die Prallelektrode aus Aluminium aufgeschraubt ist. Durch eine Kupfer- 
scheibe, die mit Rohr und Stopfen verlötet ist und gut plan an der Aluminiun- 
platte anliegt, wird ein 
sir  withrleistet. Die Kühlung 
04° des anderen Endes de 
ferrohres erfolgt durch 
60:  @inen Kupferklotz, in den 
eg eingelötet ist. In einer 
Ausdrehung dieses Kupfer- 
klotzes liegt eine Kühl 
_ schlange aus Glasrohr, die 
mit dem Klotz durch ein 
Gemisch aus Wasserglas 
und Talkum innig verkittet 
ist. In dem Kupferklotz 
liegt die Lötstelle II, sodaß 
Sis eine direkte Erwärmung 
durch Hochfrequenz auch 
hier nicht zu befürchten ist. 
0 10 i daß bei den oben beschrie- 
Abb. 5. Tomperaturabfall bei Eichmessungen (Z) und Denen Eichungen und da 
Pendelvervielfachung (P) Messungen mit Pendel- 
strom die gleichen ther- 
mischen Verhältnisse vorlagen, wurde in beiden Fällen der vollständige Ver- 
lauf von # nach Aufhören der Leistung L beobachtet. Die Integration der 
Gl. (4) ergibt ja 


> 
A 
> 
7 
f 


In (Or — =] 
Bei logarithmischer Auftragung von (9, — 9) bzw. (Urn — Ut) über 
der linearen Zeitabszisse t ergeben sich also Geraden; ihre Neigung (s. Abb. 5) 
stimmte bei den verschiedenen Eichmessungen (E) und Pendelstrommessungen 
(P) gut überein!®). 

Wie oben beschrieben, muß die Pendelleistung aus einer Differenzmessung 
gewonnen werden. Es ist also dafür zu sorgen, daß bei den gleichen Hoch- 


t + const. (7) 


1) Aus im ganzen 19 solcher Temperaturabfall-Messungen wurde der Quotient 
4/C experimentell zu 0,00385 + 1 %, sec ermittelt. 
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frequenzbedingungen einmal ein Pendelstrom fließt und einmal nicht. Die 
Möglichkeit dazu ist dadurch gegeben, daß zur Ausbildung der Pendelver- 
vielfachung, jedenfalls bei den hier verwendeten Versuchsanordnungen, ein 
Magnetfeld von etwa 100 Oersted notwendig ist. Durch die Veränderung 
des Magnetfeldes konnte so die Pendelvervielfachung ermöglicht und auch 
wieder unterbunden werden. Die Messungen gingen nun so vor sich, wie es 
in der Abb. 6 als Beispiel dargestellt ist. Hier wurde über der Zeitachse 
aufgetragen: die Hochfrequenz- , 
spannung in einem willkir- ‚sch 
lichen Maßstab, der Auffängerstrom ?* / 
1,2) als Maß für die Pendelverviel- 


fachung und die auftretende Thermo- 
spannung Upp. U 


Die Schwierigkeit bei diesen 
Messungen bestand vor allem in | 


der Konstanthaltung des Pendel- 603 
stromes über hinreichend lange Zeit. J, 
Offenbar liegt das daran, daß der 130 
Sekundäremissionsfaktor der Prall- pA 
elektroden, der ja einen bestimmen- 1 20 
den Einfluß auf die Zahl der Pendel- 
elektronen hat, weitgehend von dem ’ 
Zustand der Elektrodenoberfläche q Un, 
abhängt. Dieser aber kann sicher . x 

durch die auftreffenden Elektronen © 
im Zusammenwirken mit den un- pee Pendelvervielfachung Hochfrequenz 


vermeidbaren spurenweisen Verun- | 
reinigungen in der Röhre verändert Abb. 6. Zeitlicher Verlauf von Auffänger- 


strom /,, Hochfrequenzspannung U,r und 


Pendelvervielfachung die Hochfrequenzspannung in der Röhre infolge des ver- 
mehrten Energieentzuges ein wenig absinkt. Gleichzeitig beginnt eine Erwär- 
mung der Elektrode, die sich aber wegen der thermischen Trägheit der Anlage 
über einen längeren Zeitraum erstreckt, bis endlich eine konstante Temperatur 
erreicht wird. Nach dem Abschalten der Pendelvervielfachung sinkt die Tempe- 
ratur wieder langsam ab und befolgt dabei, wie auch aus der Abb. 5 zu ersehen 
ist, sehr gut das Exponentialgesetz. Es wird ein konstanter Endwert der 
Temperatur erreicht, der dem Ausgangswert vor Beginn der Pendelverviel- 
fachung entspricht und durch Hochfrequenzerwärmung bedingt ist. Nach 
Abschalten der Hochfrequenzspannung gleicht sich auch diese Temperatur- 
differenz langsam aus. > 
ce) Ergebnisse 


Die durch den Pendelstrom hervorgerufene Temperaturerhöhung erzeugte 
in obigem Beispiel (Abb. 6) eine Thermospannung von 102 Skt. Etwa in 


17) Siehe diesbezüglich Fußnote!?) oben S. 148. ee 
_ 18) So bildet sich auf den Platten nach längerem Betrieb der Pendelvervielfachung 
ein Belag mit wechselnden Farben, sobald Spuren von Kohlenstoff in der Röhre an- 


wesend sind. Diese Beobachtung wird auch von König und Helwig, Z. Physik 129, 3 u 


491 (1951) beschrieben. 


: 
tisch 
It ist, 
mm 
liesem 
n ver- 
uf das ah 
upfer 
inium- | 
d ein 
Finheitep 
ng ge- 
ihlung 
> des 
durch 
3 
einer we 
ole 
Kühl 
T, die | 
h ein 
>f 
erglas 
kittet a 
= 
über, 
chrie- 
I den 
ndel- 
ther- 
Ver- 
n der 
(7) 
über 
b. 5) 
sung 
och- 


Ammalen der Physik. 6. Folge. Band 18.1956 


dieser Größe lager: auch die Ergebnisse der anderen Messungen. Die Hoch 
frequenzspannung wurde bei den Messungen so gewählt, daß sich möglichst 
stationäre Zustände der Pendelvervielfachung einstellten. Dies war ke 
Spannungen etwa in der Mitte zwischen Einsatz- und Aussetzspannung der 
Fall. Eine Spannungsmessung mit lose angekoppeltem Detektor genügt 
für diese Versuche, da ja auf Grund der Ergebnisse der Einsatzspannung: 
messung!) die relative Skala der Detektormessung mit der jeweiligen Ein- 
satzspannung geeicht werden konnte. Mit der obigen Beziehung (6) ergibt 
sich eine Pendelleistung von 1,57 Watt. Falls die Elektronen die von der 
Theorie geforderte einheitliche Auftreffenergie hätten, wäre damit die Pendel- 
stromstärke eindeutig gegeben. Im zweiten Teil der Arbeit wird gezeigt, daß 
diese Voraussetzung zwar nicht vollkommen erfüllt ist, jedoch ein deutliches 
Maximum der Energieverteilung existiert, welches den in Phase befindlichen 
Elektronen entspricht. Dieses Maximum liegt bei den größten überhaupt 
vorkommenden Elektronengeschwindigkeiten. Die schnellen Elektronen 
tragen aber ohnehin den Hauptteil der Pendelleistung. Also kann man ohne 
wesentlichen Fehler annehmen, daß die durchschnittliche Elektronenenergie 
diesem festgestellten Maximum entspricht. Dieses Maximum der Auftreff- 
energie liegt bei einer Feldlänge a = 4 cm und einer Hochfrequenzspannung 
von 1,2 x Einsatzspannung, den Versuchsbedingungen der in Abb. 6 dar- 
gestellten Messung, bei 490 Volt. Damit berechnet sich der Pendelstrom zu: 
1,57 Watt „A 
490 Vols 

und die in jeder Periode bewegte Ladungsmenge beträgt: se un 

3,2 mA 

80-108 sect 0,40 - 10-10 Coul. 
Es ist klar, daß eine solche Ladungsmenge schon deshalb nicht geschlossen 
an einer Elektrode ankommen kann, weil sie auf Grund der Coulombschen 
Abstoßung die Tendenz zu einer räumlichen Expansion haben muß. Die 
Ankunftsphase wird dadurch ,,defokussiert‘‘ und die Auftreffenergie kann not- 
wendigerweise nicht mehr einheitlich sein, sondern wird eine Verteilung um 
einen Mittelwert herum annehmen. Zu einer groben Abschätzung dieses 
Effektes sei — sehr vereinfachend — angenommen, daß die Ladung auf einer 
Kugel von dem Durchmesser der Prallelektroden konzentriert ist. Ein Elek- 
tron an der Oberfläche dieser Kugel würde dann durch die Ladung eine 
Beschleunigung von 3,9 - 101% cm/sec? erfahren und hätte in 0,6 - 10° see 
(der halben Hochfrequenzperiode) eine Energie von 16eV erreicht. Die 
Streuung der Auftreffenergie muß danach jedenfalls größer als + 16 Volt 


€. Die Energieverteilung der Pendelelektronen 
Aus der Theorie der Pendelvervielfachung!) ergibt sich, daß die pendelnden 
Elektronen auf eine bestimmte Startphase hin fokussiert werden sollten. Das 
würde bedeuten, daß in jeder Hochfrequenzperiode nur einmal in einem sehr 
kurzen Zeitraum Elektronen auf jede der Prallelektroden auftreffen dürften. 
Tatsächlich gibt es aber mindestens zwei Gründe für eine Defokussierung: 
1. Die Streuung der Austrittsgeschwindigkeit der Sekundärelektronen, die 
selbst bei völlig gleicher Startphase eine Variation der Ankunftsphase her 
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vorrufen muß und die in den theoretischen Ableitungen nicht berücksichtigt 
wurde, und 2. die bereits oben erwähnte gegenseitige Abstoßung der Elek- 
tronen. 

Diese beiden Faktoren bewirken in jeder Periode die Abwanderung eines 
bestimmten Teiles der pendelnden Elektronen in einen „unwirksamen Bereich“. 
Die von diesen Elektronen ausgelösten Sekundärelektronen fallen mit ge- 
ringer Auftreffenergie auf die Ausgangsplatte zurück und werden im folgenden 
abgekürzt als „Rückkehrelektronen‘“, im Gegensatz zu den die Gegenplatte 
erreichenden ‚„Ankunftselektronen‘‘, bezeichnet. 

Wie bereits hergeleitet!), ist die Auftreffenergie der Ankunftselektronen 
eindeutig von der Startphase abhängig, und zwar gilt: 


« 


U, 
= (cos Po — cos Pa) + V (8) 


Gl. (8) stellt die auf eine Hilfsgröße U* von der Dimension einer Spannung 
bezogene Auftreffenergie U, der Pendelelektronen als Funktion der Start- 
und Ankunftsphase (9, 94), der Anfangsgeschwindigkeit der Sekundär- 
elektronen (y) und der Hochfrequenzspannung (ß) dar”). 

Dabei ist noch zu berücksichtigen, daß die Ankunftsphase, die ja eine 
Funktion der sonstigen in Gl. (8) auftretenden Größen ist, in jedem be- 
sonderen Falle eliminiert werden kann. Bei fester Hochfrequenzspannung ist 
also, da die Austrittsenergie bzw. richtiger die Austrittsenergieverteilung ja 
feststeht, die Auftreffenergie nur eine Funktion der Startphase. Diese Be- 
ziehung ist umkehrbar, so daß aus der gemessenen Auftreffenergie Rück- 
schlüsse auf die Startphase möglich sind. Insbesondere bedeutet eine einheit- 
liche Auftreffenergie oder die Häufung der Auftreffenergie um einen bestimmten 
Wert herum, daß es eine bevorzugte Startphase gibt. 

Die Kenntnis der Auftreffenergie ist besonders aus folgendem Grunde von 
Bedeutung: Wie schon oben erwähnt, ist von Gill und v. Engel eine 
Theorie der Pendelvervielfachung aufgestellt und von Hatch und Williams 
erweitert worden, die in wesentlichen Punkten von der in!) gegebenen kine- 
matischen Theorie abweicht. Die genannten Autoren gehen von der Voraus- 
setzung aus, daß eine Pendelvervielfachung nur dann möglich ist, wenn die 
Auftreffenergie der Sekundärelektronen groß genug ist, um den Sekundär- 
emissionsfaktor 6, der ja bekanntlich eine Funktion der Auftreffenergie ist, 
über 1 anwachsen zu lassen. Bei Steigerung der Hochfrequenzamplitude 
wächst die Auftreffenergie der Elektronen, und die Pendelvervielfachung 
sollte nach dieser Theorie gerade dann einsetzen, wenn die obige Bedingung 
erfüllt ist. Der Einsatz der Pendelvervielfachung sollte demnach bei allen 
Feldlängen bei der gleichen Auftreffenergie der Pendelelektronen von etwa 
100 Volt erfolgen. Nach der in!) dargestellten Auffassung sollte dagegen bei 


19) Unter Start- bzw. Ankunftsphase ist der Zeitpunkt des Startes und der Ankunft 
der Pendelelektronen, bezogen auf die Phase der Hochfrequenzspannung, zu verstehen. 
% = ty (t, = Startzeitpunkt) und g,= ot, (t, = Ankunftszeitpunkt). In U*= 
m @® ate sind die Feldlänge a, die Frequenz ® sowie Ladung und Masse des Elektrons 
enthalten. Alle in den Gleichungen auftretenden Spannungen werden auf U* bezogen. 
Gl. (8) und folgende gelten dann ganz allgemein für alle Feldlängen und Frequenzen. 
So ist: 8 = U,/U* (U, = Amplitude der Hochfrequenzspannung U = U, cos wi) und 
y = U, U*(U, = Startenergie der Sekundärelektronen in eV). 
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größeren Feldlängen auch schon beim Einsetzen der Pendelvervielfachung 
eine größere Auftreffenergie vorhanden sein. Die Messung derselben kann algo 
eine experimentelle Entscheidung zwischen den beiden Theorien liefern. 


b) Aufbau der Apparatur und Messungen 


Zur Messung der Auftreffenergie wurde die Gegenfeldmethode nach 
Lenard verwendet, deren Prinzip in Abb. 7 dargestellt ist. Der Strom vom 
Elektronen verschiedener Geschwindigkeit trifft von links auf die Blende 
und wird durch Variation der Gegenspanung U, in seine Bestandteile zerlegt. 
Es können immer nur die Elektronen den Auffänger erreichen, deren Energie 

größer oder gleich der jeweiligen Gegen- 


8, 8, 8, spannung U, ist. Wird nun J, als Funk- 
n 9 22 tion von U, aufgenommen, so stellt diese 

4 y ; Kennlinie also eine Integralkurve dar, aus 
- A— der durch Differentiation die Energie 

N ; y verteilung gewonnen werden kann. Bei 
A; RER den hier beschriebenen Versuchen bestand 


ss ened die Blende B, aus einer Bohrung in der 
——ralllelektrode Z,. (Siehe Abb. 8.) Der Aut 
al ts fanger A ist mit einer zweiten Blende B, 
Hm versehen. Der Zweck der Blende B, wird 
ah‘ weiter unten bei der Beschreibung des 


=> + - Ya. vollständigen Aufbaus erläutert. Mögliche 
Abb. 7. Prinzipbild der Gegenfeld- Fehlerquellen sind in der Hauptsache 
methode Sekundärelektronenemission und Linsen- 


wirkung. Die Sekundäremission an der 
Blende B, wird durch möglichst scharfe Ausblendung des Elektronenstrahles 
an B, vermieden und die Sekundäremission im Auffänger selbst durch Be- 
rußen weitgehend herabgesetzt. Außerdem wird die Blende B, so klein ge- 
macht, wie es mit der ersten Bedingung irgend verträglich ist, so daß die im 
Auffänger ausgelösten Elektronen möglichst in demselben bleiben. Auch am 
Rand der Blende B, werden Sekundärelektronen ausgelöst. Dies läßt sich 
durch eine mechanisch saubere Gestaltung des Blendenrandes in gewissem 
Grade einschränken, aber nicht ganz unterbinden. Dafür stören die dort 
ausgelösten Elektronen auch verhältnismäßig wenig, denn ihre Energie ist 
nur klein (einige Volt; wenn man von den sogenannten ,,reflektierten‘‘ Elek- 
tronen absieht, die aber ihrerseits wieder fast die gleiche Energie wie die 
Primärelektronen haben und daher das Meßergebnis auch nicht wesentlich 
verfälschen können). Die langsamen ‚echten‘ Sekundärelektronen sind 
also nicht mit den wesentlich schnelleren Pendelelektronen zu verwechseln. 
Sie haben nur einen Einfluß auf den Verlauf der Gegenspannungskurve bei 
kleinsten Spannungen, der für die Bestimmung der Auftreffenergie nicht 
weiter interessiert. 


Was die elektronenoptische Beeinflussung des Elektronenstrahles anbe- 
trifft, so ist diese bekanntlich um so geringer, je enger sich dieser der Symme- 
trieachse des Systems anschmiegt. Daher mußte der Durchmesser der Blende 
B, gegenüber früheren Messungen!) auf 0,2 mm verkleinert werden. Der 
Auffängerstrom sank dadurch in die Größenordnung 10 bis 10° A. 
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Dadurch wurden besondere Abschirmungen notwendig, so daß die ganze 
Versuchsanordnung schließlich nach Abb. 8 gestaltet wurde. 

Die Hochfrequenzenergie wurde in gleicher Weise, wie bei früheren Mes- 
sungen!) vom Sender über ein abstimmbares Leitungssystem (in Abb. 8 
stark gezeichnet) der Pendelstrecke und der zur Spannungsmessung dienenden 
Diode LG 2 zugeleitet. Der Abstand a der Prallelektroden E,, E,, die Feld- 
länge, ist veränderlich. Eine Magnetspule M sorgt für die zur Ausbildung 
der Pendelvervielfachung notwendige Bündelung der Elektronen. Die Prall- 
elektrode HE, trägt die GegenfeldmeBanordnung. Sie besteht aus einer äußeren 


Zur fichung an Stelle von E, 


G sb Jum hier 
50 M2 


Abb.8. Versuchsanordnung zur Messung der Energieverteilung. Links unten maß- 
stabsgerechte Darstellung des Blendensystems 


Abschirmung S, welche die MeBstrecke gegen Hochfrequenzspannungen ab- 
schirmen soll. Ihre Vorderseite bildet die Prallelektrode E,. Eine in diese 
Prallelektrode eingelassene Silberscheibe (in der vergrößerten Nebenzeich- 
nung in Abb. 8 eng schraffiert gezeichnet) ist mit einer Bohrung von 0,2 mm 
Durchmesser versehen, die die Blende B, bildet. Der eigentliche Auffänger 
A ist nach dem Vorbild von Boersch™) von einem Schutzkäfig K mit der 
Blende B, umgeben. Ebenfalls ist die Zuleitung zum Auffänger abgeschirmt. 
Diese Abschirmung ist lückenlos vom Schutzkäfig K bis zu einem Abschirm- 
käfig, der in Abb. 8 dünn gestrichelt ist, durchgeführt. Innerhalb dieses Ab- 
schirmkäfigs befindet sich das stromanzeigende Galvanometer mit seinen 
Schaltelementen. Dabei ist, wie in Abb. 8 angedeutet, auch dafür gesorgt, 
daß die Abschirmung auch an der Stelle der Durchführung durch die Glas- 


*) H. Boersch, Z. Physik 189, 115 (1954). 
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wand der Röhre nicht unterbrochen ist. Diese ganze sorgfältige Abschirmung 
hat eine doppelte Aufgabe. Einmal soll sie die Kapazität zwischen dem Auf. 
fänger mit seiner Zuleitung und der Erde verringern, damit die bei Änderung 
der Auffängerspannung fließenden Aufladungsströme, die in der Größen- 
ordnung der zu erwartenden Auffängerströme liegen, die Messung nicht 
stören. Ebenso fließen bei der verwendeten Abschirmung eventuelle Isola- 
tionsströme nieht mit über das Galvanometer, denn zwischen Auffänger und 
der Abschirmung liegt so nur die kleine Spannung, die bei dem jeweiligen Auf- 
fängerstrom an dem Galvanometer + Vorwiderstand abfällt. Dieser Vor- 
widerstand von 6 MQ und der Überbrückungskondensator waren notwendig, 
um das Galvanometer von den Folgen von Gaseinbrüchen und vor Hoch- 
frequenzspannungen zu schiitzen, die unvorhergesehen in das System einge- 
koppelt werden kénnten. 

Um aus den in dieser Anordnung gemessenen Gegenspannungskurven 
zuverlässig auf die Energieverteilung der Pendelelektronen schließen zu können, 
wurde die Anlage noch geeicht. Diese Eichung geschah in der Weise, daß ein 
Elektronenstrom mit bekannter Energieverteilung auf die Elektrode E, ge- 
schickt und die Gegenspannungskurve aufgenommen wurde. Bei unseren 
Genauigkeitsanforderungen genügte dafür der von einer Glühkathode aus- 
gesandte Elektronenstrom. In ähnlicher Weise wie bei der Eichung der Lei- 
stungsmessung wurde auch 
jetzt die Gegenelektrode 
durch einen Glühdraht @ 
ersetzt und auf eine nega- 
tive Spannung gegen Erde 
je gebracht. In der Abb. 8 
ist die Eichschaltung links 
-- en (Die 
Hochfrequenzanlage fällt 
Uy~200 Volt bei der Eichung natürlich 

fort.) Es ergaben sich 
die Gegenspannungskurven 


U,*114 Volt der Abb. 9. Auf eine Diffe- 

\ rentiation der Kurven 

Up wurde verzichtet, da alles 

0 50 100 150 200 Volt 250 Interessante auch so zu 
Abb. 9. Kontrollmessungen zur Funktion der Gegen- sehen ist. Der Abfall bei 
feldanordnung kleinen Gegenspannungen 


entspricht den an der 
Blende B, ausgelösten Sekundärelektronen. Dann bleibt der Auffänger- 
strom über einen weiten Spannungsbereich konstant, das heißt, es sind 
keine Elektronen mit einer diesen Spannungen entsprechenden Energie 
vorhanden. Der Abfall des Auffängerstromes erfolgt, wie zu erwarten, 
erst bei einer der Zugspannung U, entsprechenden Gegenspannung. Der 
Abfall des Auffängerstromes erstreckt sich über einen gewissen Gegenspan- 
nungsbereich. Die Geschwindigkeitsverteilung der Glühelektronen ist dem- 
nach etwas breiter, als sie theoretisch auf Grund der thermischen Energie- 
verteilung und des Spannungsfalls längs des Glühdrahtes sein sollte. Die 
erstere beträgt etwa 1, der zweite rund 4 Volt, so daß mit einer Streuung 
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von 5 Volt zu rechnen wäre. Tatsächlich beträgt die Streuung im ungünstigsten 
Falle um 10 Volt. Die Verbreiterung könnte durch an der Blende B, mit 
geringem Energieverlust „reflektierte‘“ Elektronen hervorgerufen sein. Ein’ 
Fehler von 5 Volt ist aber bei der Messung der Energie der Pendelelektronen, 
die bei einigen 100 Volt liegt, durchaus zuzulassen, so daß die Apparatur in 
dieser Form als brauchbar anzusehen war. 


Die Ergebnisse der verschiedenen Meßreihen sind in den Abb. 10 und 11 
dargestellt. Den allgemeinen Verlauf einer Gegenspannungs-Auffängerstrom- 
Kennlinie zeigt Abb. 10a. In Abb. 10b ist der interessante Teil der Kurve, 


nämlich der Stromabfall, 
mit besonderer Sorgfalt 
gemessen, noch einmal in 
vergrößertem Maßstab dar- 
gestellt. Dazu ist die gra- 
phische Differentiation vor- 
genommen. In der Abb. 11 
sind die Messungen bei ver- 
schiedenen Feldlängen und 
Hochfrequenzspannungen 

zusammengestellt. Im all- 
gemeinen wurden für jede 
Feldlänge 3—4 (zwischen 
Ein- und Aussetzspannung 
liegende) verschiedene 
Hochfrequenzspannungen 

gewählt. Über den Ver- 
lauf der in Abb. 10a dar- 
gestellten Kennlinie, deren 
wesentliche Züge bei allen 
Feldlängen und Hochfre- 
quenzspannungen wieder- 
kehren, ist folgendes zu 


sagen: Der Stromabfall bei kleinen Gegenspannungen (in Abb. 10a bis etwa 
75 Volt) ist zum Teil durch die an der Blende B, ausgelésten Sekundärelek- 
tronen zu erklären. Wie aber ein Vergleich mit der Eichkurve Abb. 9 zeigt, 
sind außer diesen Sekundärelektronen noch weitere langsame Elektronen in der 
Energieverteilung vorhanden. Bei der Eichung, wo nur die an B, ausgelösten 
Sekundärelektronen eine Rolle spielen, beginnt der waagerechte Verlauf der 
Kennlinie schon bei 30 Volt, während in Abb. 10 bis 75 Volt schnelle Elek- 
tronen in beträchtlicher Zahl zu erkennen sind. 
Spannungen ist der Verlauf der Kennlinie in Abb. 10a nicht völlig waagerecht. 
Noch ausgeprägter ist diese Erscheinung bei einigen der Kennlinien in Abb. 11 


zu erkennen. 


Eine gewisse Streuung der Ankunftsenergie wird, wie schon oben erwähnt, 
hervorgerufen durch unterschiedliche Austrittsenergie und gegenseitige Ab- 
stoßung der Elektronen. Damit kann allerdings nur diejenige Streuung der 
Auftreffenergie erklärt werden, die im Endmaximum der Energieverteilung 
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Abb. 10. a) Beispiel einer J,—U,-Kennlinie bei Pendel- 
vervielfachung; b) Vergrößerter Ausschnitt aus a) und 
Differentialkurve 
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2) = 1010 Volt, Up = 747 Volt 
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Abb. 11. 7,—U,-Kennlinien bei verschiedenen Feldlängen und Hochfrequenzspannungea 
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auftritt und die aus der Differentialkurve der Abb. 10b zu etwa 50 Volt abge- 
lesen wird. Jene langsamen Elektronen, die außer den Sekundärelektronen 
vorhanden sind, müssen, da sie bei der Eichmessung fehlen, eine spezifische 
Eigenschaft der Pendelvervielfachung sein. Es können aber keine ,,Ankunfts- 
elektronen‘‘ sein, da so geringe Energien nach der Theorie!) Gl. (8) bei diesen 
nicht vorkommen können. 

Die Theorie läßt die Möglichkeit offen,-daß Elektronen auf die Ausgangs- 
platte zurückkehren. Die Auftreffenergie dieser ,,Rickkehrelektronen“ 
läßt sich berechnen und ergibt sich ganz analog zu Gl. eae et 

U —]2 
= (cos — cos Pr) + | 
wobei Yg = wtp die Phase der Elektronenrückkehr bedeutet. Aus den 
theoretischen Überlegungen!) folgt, wie auch schon die Anschauung quali- 
tativ bestätigt, daß die Ankunft eines Elektrons auf der Gegenelektrode ent- 
weder in der Halbperiode erfolgt, die an die Halbperiode des Startes an- 
schließt, oder aber (in Ausnahmefällen) in der 3., 5. usw. Halbperiode cos gy» 
und cos, haben daher stets entgegengesetztes Vorzeichen, so daß sich in 
Gl. (8) die Beträge der cos-Werte addieren. Die Rückkehr zur Ausgangs- 
elektrode erfolgt dagegen in der gleichen Halbperiode wie der Start, oder 
aber in der 2., 4. usw. Halbperiode danach. In diesem Falle haben aber cos 9, 
und cos px gleiches Vorzeichen, so daß sich in Gl. (9) die Beträge subtrahieren. 
Es ergibt sich also, da die Größe //y nur die Rolle einer kleinen Korrektur 
spielt, daß die Energie der Rückkehrelektronen stets beträchtlich kleiner sein 
muß, als die der Ankunftselektronen. Ein stetiger Übergang von Ankunfts- 
zu Rückkehrenergie wäre nur zu erwarten, wenn cos 9, bzw. cos pr klein 
wären, d.h. wenn Ankunfts, bzw. Rückkehrphase bei 2/2 liegen würden. 
In dieser Phasenlage ist aber, wie aus den kinematischen Überlegungen!) 
hervorgeht, eine stabile Pendelvervielfachung nicht möglich. Es wird also 
stets ein sprunghafter Energieunterschied zwischen den langsamsten An- 
kunfts- und den schnellsten Rückkehrelektronen bestehen. Dies zeigen auch 
die Messungen: die Deutung der langsamen Elektronen als ,,Riickkehrelek- 
tronen‘‘ möchten wir dadurch als bestätigt betrachten. Daß es immer eine 
beträchtliche Zahl von Rückkehrelektronen geben muß, ergibt sich schon 
aus der Überlegung, daß die räumliche Abwanderung der Elektronen aus dem 
Pendelraum durch das angelegte Magnetfeld weitgehend verhindert wird. 
Die „überzähligen‘‘ Elektronen können also nur durch eine zeitliche Defo- 
kussierung „außer Tritt fallen‘ und enden nach einigen Übergängen schließlich 
als Rückkehrelektronen. 

Es sei nun noch geprüft, inwieweit die aus Abb. 11 folgenden Auftreff- 
energien die Voraussage der Theorie erfüllen. Unter der Annahme einer be- 
stimmten Austrittsenergie der Sekundärelektronen (es wurde mit Ug = 5 Volt 
gerechnet) läßt sich für jede Feldlänge und jede Hochfrequenzspannung die 
phasenreine Starphase p, berechnen. Diese ist [siehe Gl. (20) in})]: 

1 (1—2/2y) 2 
9,= tare Via + are 
Bei Phasenreinheit 
Ann, Physik. 6. Folge, Bd. 18 
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vereinfacht sich die Beziehung (8) für die Auftreffenergie zu ine SER 

gs = [BV 2- cosy, + 
Die mit dieser Formel berechneten theoretischen Werte der Auftreff- 


energie U, sind als senkrechte Striche an der Achse der Gegenspannung in 
Abb. 11 eingetragen, wobei durch die beigefügten Buchstaben die Zugehörig- 


keit zu den entsprechenden /,—U,-Kennlinien gekennzeichnet ist. Theore- 


tisch wäre zu fordern, daß diese Striche an der Stelle des steilsten Abfalles 
der entsprechenden Kennlinie liegen sollten. Das ist nicht ganz der Fall. Der 
Gang der gemessenen Auftreffenergien mit Feldlange und Hochfrequenz- 
spannung stimmt zwar gut mit dem Gang der theoretischen Werte überein, 
die Meßwerte liegen aber durchschnittlich um 20 % tiefer als die theoretischen. 

Eine Abweichung in dem beobachteten Sinne ist zu erwarten. In den theoretischen 
Überlegungen wird das Hochfrequenzfeld stets als homogen angenommen. Das ist es 


in Wirklichkeit infolge der nur endlichen Ausdehnung der Elektroden nicht. Eine weitere 
_ Fehlerquelle ist in der systematischen Ungenauigkeit der Hochfrequenzspannungsmessung 


En zu suchen. An der zur Gleichrichtung benutzten Diode LG 2 (Abb. 8) tritt, wie aus den 


Elementen der Hochfrequenztechnik bekannt, eine Spannungserhöhung auf, die aus 
der gemessenen Eigenfrequenz des aus Pendelstrecke, Zuleitungen und Diode gebildeten 
Kreises abgeschätzt und rechnerisch berücksichtigt wurde. Allerdings ist hierbei ein 
systematischer Fehler von 10% durchaus möglich, so daß bei Berechnung von U, nach 
Gl. (12) ein Fehler von 20%, entstehen kann. 


Der qualitative Verlauf der Auftreffenergie als Funktion der Feldlänge ist 


aber im Hinblick auf die oben genannte Diskrepanz zwischen der Theorie 
von Gill und v. Engel und der kinematischen Theorie!) schon interessant - 
genug. Die kleinste HF-Spannung, bei der die Auftreffenergie bestimmt 
wurde, liegt bei allen Feldlängen um den gleichen Betrag von etwa 50 Volt 
über der jeweiligen Einsatzspannung. (Unmittelbar bei der Einsatzspannung 


läßt sich infolge der Instabilität der Verhältnisse eine einigermaßen repro- 


2“ duzierbare Messung nicht durchführen.) Falls, wie von Gill und v. Engel 
i SE angenommen, stets das Anwachsen des Sekundärfaktors über 1 bei steigender 


Auftreffenergie für das Einsetzen der Pendelvervielfachung mit wachsender 


_ HF-Spannung verantwortlich wäre, müßten die bei der kleinsten HF-Span- 


nung gemessenen Werte der Auftreffenergie bei allen Feldlängen in der Nähe 
von 100 Volt liegen, wie es bei 2,0 cm Feldlänge auch der Fall ist. Bei den 


größeren Feldlängen steigern sich aber die Auftreffenergien, wie nach der 


_ kinematischen Theorie auf Grund der erhöhten HF-Spannungen zu erwarten, 
auf über 600 Volt, ehe überhaupt eine Pendelvervielfachung einsetzt. Dieses 
wäre nach der Theorie von Gill und v. Engel keinesfalls zu verstehen. Auch 
in der von Hatch und Williams erweiterten Form kann diese Theorie nicht 
ohne ad hoc — Annahmen die Einzelheiten der Pendelvervielfachung erklären. 
Dieses vermag nur eine exakte kinematische Betrachtung, die ohne spezielle 


Annahmen nur von den Grundtatsachen des Bewegungsvorganges der pen- 


delnden Elektronen ausgeht. 


E Zum Schluß möchten wir der ‚Gesellschaft von Freunden der TU Berlin“ 
auch an dieser Stelle für finanzielle Unterstützung der Arbeit unseren besten 


Dank aussprechen. Auch den Zuwendungen aus ERP-Forschungsmitteln 
verdanken wir wesentliche Förderung unserer Untersuchungen. 


Berlin-Charlottenburg, II. Physikalisches Institut der Technischen 


Universität. 


Bei der Redaktion eingegangen am 8. März 1956. — 
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ig- Wechselstrommessungen an Kadmiumsulfid-Spaltzellen 

re- 

les Von E. Klier*) 

Inhaltsübersicht 

len Scheinwiderstandsmessungen an CdS-Kristallzellen ergeben, je nach dr 
hae Herstellungsart, recht verschiedene Resultate. Bei den mit Glimmentladung 
ed gereinigten Kristallen ist eine mit Raumladungen verbundene Kapazität in 
in der weiter unten beschriebenen Anordnung nicht feststellbar. Die Messungen _ 
1Us beweisen, daß die mit Glimmentladung behandelten Photozellen bei Wechsel- _ 
‚en strom das gleiche Verhalten wie bei Gleichstrom zeigen. 5 

en In einer Reihe von Veröffentlichungen!)?) wurde von einer Methode zur 
ak sperrschichtfreien Kontaktierung an CdS-Einkristallen berichtet. Nach 
mt dieser Methode werden die Kontaktflächen der Kristalle vor der Aufdampfung 


olt der Elektroden eine längere Zeit dem Einwirken einer Glimmentladung aus- 
gesetzt. Die auf diese Weise vorbehandelten Zellen zeichnen sich durch 


| folgende Merkmale aus: 

ei 1. Die Strom-Spannungscharakteristiken weisen eine strenge Giiltigkeit 
at des Ohmschen Gesetzes auf. 

hie 2. Es werden keine Formierungsprozesse, d.h. keine zeitliche Änderung 
ne der Eigenschaften bei Strombelastung, beobachtet. 


hip 3. Der Photostrom ist einer Potenz der Lichtintensität proportional, 
deren Exponent mit steigender Intensität von 1 auf 0,5 abnimmt. 


” i 4. eg Rauschstrom entspricht dem Rekombinationsrauschen nach 
isolf*). 
od 5. Bei Lichtabtastung findet keine Bevorzugung der Kathode statt. 
oli 6. Der Spannungsabfall zwischen den Elektroden verläuft linear ohne 
ht Spriinge an den Elektroden. 
we Diese Eigenschaften der Zellen werden von Buttler und Muscheid!) so 
lle gedeutet, daß durch Beschießen mit langsamen Elektronen dünne Adsorp- 
e tionsschichten beseitigt werden, so daß bei nachfolgendem Aufdampfen der 
Elektroden ein inniger Kontakt zwischen Metall und Halbleiter entsteht. 
“ Sie nehmen an, daß dabei unter Umständen mit abnehmendem Abstand 
n zwischen Metall- und Halbleiter-Oberfläche das Voltapotential und somit 
en — 
In *) Z.Z. Fysikälni ustav Karlovy university, Praha. 
1) W. M. Buttler, W. Muscheid, Ann. Physik (6) 14, 215 (1954); 15, 82 (1954); 
Arbeitstagung Festkörperphysik II. 1955, S. 102. 
en *) I. Poppe, Arbeitstagung Festkörperphysik II, 1955,S.114. 


®) J. H. Gisolf, Physica 15, 825 (1949). 
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die Randverbiegung des Leitungsbandes verschwinden könnten. Dadurch 
würden auch die Raumladungen und die für die Gleichrichterwirkung maß- 
gebende Verarmungsrandschicht verschwinden; der Kontakt wird ‚ohmisch“, 

Falls diese Deutung zutreffen sollte, müßte sich das Vorhandensein bzw. 
das Verschwinden der Raumladungen bei Wechselstrommessungen leicht 
feststellen lassen. 

Das Ziel dieser Arbeit ist, das Verhalten der Zellen mit ohmschen (0)- 
als auch mit nichtohmschen (NO)-Kontakten gegenüber Wechselstrom zu 
untersuchen ®). der 

2. Versuchsmaterial und 

standen drei Sorten von Spaltzellen zur Verfügung: arty 
1. Zellen mit zwei ,,NO“-Elektroden, 
zer 2. Zellen mit einer „O“- und einer ,,NO“-Elektrode, 2 

3. Zellen mit zwei „O“-Elektroden. 
Die Zellen wurden überschlägig auf Gleichrichtereffekt und ohmsches 
Verhalten bei Gleichstrom innerhalb zwei bis drei Zehnerpotenzen überprüft. 

Zur Messung des komplexen Wechselstromwiderstandes wurde eine ein- 
fache Brücke nach Rose®) aufgebaut (Abb. 1). Die festen Zweige der Brücke 


£ 


Abb. 1. Schaltung der Meßbrücke 


bildeten zwei Widerstände zu je 150 Ohm, das Nachahmungsglied bestand 
aus einem üblichen linearen Potentiometer zu 50 kQ2 und einem Drehkonden- 
sator von 400 pF in Parallelschaltung. Dadurch wurde die Auswahl der 
Kristalle auf diejenigen beschränkt, deren Widerstand bei der gegebenen 
Lichtanregung 50 kQ nicht überstieg. Der Meßbereich konnte auf Kosten 


*) Daß Sperrschichten auf unbehandelten Kontaktflächen entstehen, wurde durch 
Messungen von B. Perlwitz5) an CdS-Flächenzellen bewiesen. 

5) B. Perlwitz, Raumladungserscheinungen in Diplomarbeit, 
II. Phys. Inst. d. Humboldt-Universität, Berlin. 

6) F. Rose, Ann. Physik (6) 9, 97 (1951). tHe 
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der Meßgenauigkeit auf etwa 500 k2 erweitert werden, in dem man parallel 
zur Zelle einen bekannten Widerstand von etwa 50 kQ schaltete. Dieses 
Verfahren wurde nur bei Messungen mit niedrigen Intensitäten verwendet. 

Als MeBwert C,, wurde im weiteren die Differenz zwischen der gemessenen 
Kapazität und der Anfangskapazität der Schaltung, des Zellenhalters und der 
leeren Kristallunterlage bezeichnet. 

Die Gleichstrom-Vorspannung wurde der Zelle über die Ausgangsdiagonale 
der Brücke zugeführt. Dadurch wurde erreicht, daß die Kapazität des Gleich- _ 
stromkreises in die Messung nicht einging. Damit der Ausgang wechselstrom- __ 
mäßig nicht kurzgeschlossen wurde, wurden im Gleichstromkreise zwei 
Drosselspulen zu 0,1 H und 10 H eingeschaltet. 

Die Speisung der Brücke erfolgte durch einen Mittelfrequenzgenerator 
über einen Permalloy-Umspanner. Die Spannung an der Brücke betrug 
höchstens 100 mV , an der Probe also 50 mV. In keiner Messung wurde irgend- 
eine Abhängigkeit von der Wechselspannungsamplitude gefunden. Es wurde 
im Frequenzbereich 100 Hz bis 100 kHz gemessen. 

Als Null-Indikator wurden ein Meßverstärker und ein Röhren-Millivolt- 
meter benutzt. Bei Messungen unter 10 kHz wurde ein Kopfhörer verwendet. 
Bei Benutzung von Gleich-Vorspannung war es notwendig, den niederfre- 
quenten Rausch der Zelle durch ein zwischen Verstärker und Voltmeter einge- 
schaltetes Hochpaßfilter zu unterdrücken. (In Abb. 1 nicht eingezeichnet.) 

Die Eichung der Brücke wurde mit einem Satz von festen Widerständen 
und Kapazitäten, deren Werte bekannt waren, durchgeführt. Die Meß- 
genauigkeit der Anordnung betrug 1 bis 3%, die minimal noch feststellbare 
Differenz der Kapazitäten lag etwas unter 1 pF, die Frequenzabhängigkeit 
der Brücke im Bereiche 100 Hz bis 100 kHz dagegen innerhalb der ange- 
führten Meßfehler. 

Die Beleuchtungseinrichtung bestand aus einer Quecksilber-Höchstdruck- — 
lampe HBO 200 mit notwendigen Linsen und Graufiltern. Mit Interferenz- 
filtern wurden die Linien 436 mu und 578 mu ausgesondert. Die maximalen => 
Werte der Intensität betrugen 5 - 1016 bzw. 3,7. 1016 Quanten pro cm? sec. 


3. MeBergebnisse 


Aus umfangreicherem Versuchsmaterial werden hier die Meßergebnisse 
herausgegriffen, die einen Übergang von ausgeprägtern Gleichrichtereffekt 
zu ohmschem Verhalten der Elektroden veranschaulichen. 

Wie zu erwarten war, verhalten sich die Zellen der Gruppe 2., d. h. solche 
mit einer „O‘- und einer ,,NO‘‘-Elektrode, qualitativ wie Gleichrichter. Bei 
allen diesen Zellen ist ein kapazitiver Anteil des Wechselstromwiderstandes 
gut meßbar. Bezeichnet man die gesamte Impedanz der Zelle mit W 


Z=R-j-X, 
so kaun man den reellen und imaginären Teil aus den Meßwerten R,, und 
gemäß nachstehenden Formeln errechnen: 

Rn 


= 165 
ırch 
1aß- 
oh‘, 
ZW, 
icht 
0). 
zu 
hes 
üft, 
cke 
nd 
n- 
er 
ch 
it, 


von 


307 


26} 


22} 


In Abb. 2 ist der Verlauf der Werte R und X in Abhängigkeit von der 


‘ 
at 
a 
— 10° He 


Abb. 2. Abhängigkeit des Wirk- und Scheinwiderstandes 


der Frequenz, Zelle 16, ohne Vorspannung, bei An- 
regung im Grundgitter und im Ausläufer 


A= 


av 
10? n* 10° Hz 


Abb.3. Abhängigkeit des Wirkwiderstandes R von der 
Frequenz. Kurven I bis III, Zelle 15, Anregung mit 
weißem Licht, Kurven IV und V, Zelle 59, Anregung 


578 mu. Die Vorspannungen sind bei den Kurven 
angegeben 


BES Beim idealen Gleichrichter, bei dem eine Elekt 
sperrfrei ist, erhält man als Scheinwiderstandsdiagramm einen Halbkreis 
mit dem Mittelpunkt an der R-Achse. Bei unseren Zellen weichen die Dia- 


Frequenz halblogarithmisch aufgetragen. Man sieht, daß der Wirkwiderstand 
bei Anregung im Grundgitter (436 ma) und im Ausläufer (578 ma) mit der 


Frequenz monoton ab- 
nimmt, daß dagegen der 
Scheinwiderstand bei Aus- 
läuferanregung ein Maxi- 
mum hat. Da unzweifel- 
haft bei genügend hohen 
Frequenzen in allen Fällen 
der Scheinwiderstand ab- 
nehmen muß, ist anzu- 
nehmen, daß auch die X- 
Kurve bei Grundgitter- 
anregung durch ein Maxi- 
mum bei höheren Fre- 
quenzen hindurchgeht. 
In den Abb. 3 und 4 
ist der Frequenzgang der 
R- und X-Komponenten 
der Impedanz bei einer 
Zelle der Gruppe zwei bei 
Weißanregung (Kurven I 
bis III) und einer älteren 
ohmschen Zelle, deren 
Kontakte jedoch nicht 
ganz einwandfrei waren, 
aufgetragen (Kurven IV 
und V). Durch Gleich- 
Vorspannung kann man 
offensichtlich den Verlauf 
der Kurven beeinflussen. 
Es tritt eine Abflachung 
der Kurven auf, wenn 
die ,, NO‘‘-Elektrode posi- 
tiv ist. Bei der schwach 
nichtohmschen Zelle ohne 
Vorspannung erhält man 
einen ähnlichen Verlauf 
wie im vorhergehenden 


Falle der Grundgitter- 
anregung. Das deutet 
darauf hin, daß die 


Raumladungsschicht bei 
Grundgitteranregung we- 
niger ausgebildet ist. 


rode sperrend, die andere 


J 
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gramme sehr weit vom Halbkreis ab (Abb. 5), besonders bei Grundgitteran- 
regung. Etwas besser sehen die Diagramme bei starker Anregung mit weißem 


gegeniiber der Kreisform 
eine Erhöhung der Schein- 
komponente X und eine 
Verlegung des ersten 
Sehnittpunktes mit der 
R-Achse nach links mit 
steigender Frequenz sicht- 
bar (d. h. den kleinen 
Werten von R zu). 

Die Verringerung der 
beiden Komponenten bei Abb.4. Abhängigkeit des Scheinwiderstandes X von der 
positiver Vorspannung Frequenz. Bezeichnung wie in Abb. 3 
der „NO‘“-Elektrode deu- 
tet darauf hin, daß dabei die Kapazität der Sperrschicht vergrößert und 
der Widerstand verringert wird, was mit der Tatsache, daß CdS ein n-Leiter 
ist, in Einklang steht. 

Die Messungen des effek- 
tiven Widerstandes R,, und 
der effektiven Kapazität C re 
spannung bei einer bestimm- 
ten Frequenz geben weiter 


R 


Aufschluß darüber, wie der 0 20 70 42 
Übergang von einer sperren- pbb, 5. Scheinwiderstandsdiagramme, Zelle 16, 
den Zelle zur sperrfreien vor ohne Vorspannung 


sich geht. In Abb. 7 ist eine 
solche Messung an einer Zelle der Reihe 2 aufgetragen. Die Kurven zeigen 
deutlich das Zusammenschrumpfen der Sperrschicht an der ,,NO“‘-Elektrode, 
wenn diese positiv ist. Die Form 
der Kurven ist qualitativ völlig 
gleich mit solchen, die an Flächen- 
Gleichrichtern gemessen wurden 
(siehe z. B. H. Pfotzer’)). 

Die Kurven 3 und 3’ in 
Abb. 8 gehören der schon früher 
erwähnten Zelle mit gereinigten 
Kontaktflächen an, die sich im 
Laufe der Zeit verschlechtert 
Abb. 6. Scheinwiderstandsdiagramme, Zelle 15, 

‘ Anregung mit weißem Licht 
hergegangenen Beispiel, nur ist 
die Asymmetrie bei Vorspannungen verschiedener Polarität nicht so stark 
ausgeprägt. Die Sprünge in den Werten R, und C,,, bei Zurückkehren zur 


?) G. Pfotzer, Z. Naturforsch. 4a, 691 (1949). 
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mit höheren Gleichspannungen. 
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Vorspannung Null sind Folgen der Formierungsprozesse bei Messungen 


Die Kurven 1 und 2 gehören ebenfalls einer Zelle mit gereinigten Kontakten 
an, die jedoch von Anfang an nicht einwandfrei waren. Zuerst wurde die 


Kurve laufgenommen. Nach 
einem Formieren, das in 
langerer Belastung der Zelle 
mit 8 V Gleichspannung in 
beiden Richtungen bestand, 
konnte die Anderung des 


Widerstandes wesentlich ver- 


ringert werden. Vor dem 
Formieren nahm die Kapa- 
zitat mit steigender Span- 
nung in beiden Richtungen 
von 2,5 pF auf 1,5 pF ab, 
Nach dem Formieren war 


MeBfrequenz 100 kHz 
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Abb. 8. Abhängigkeit der ER R,„ und 

C,, von der Vorspannung bei den Zellen 3 und 

59. Anregung A = 578 mu, Meßfrequenz 100 kHz 
(Zelle 3) und 50 kHz (Zelle 59) 


innerhalb der Fehlergrenzen überein. 


nannten Punkte, 


stand der Zelle. 


20 20 10 02 04 06 V 


Abb. 7. Abhängigkeit der MeBwerte R,, und C,, von 
der Vorspannung bei Zelle 16. Anregung im Ausläufer, 


die Kapazität 1,5 pF unab- 
hängig von der Spannung. 

Bei dieser Zelle konnte 
man bereits eine angenäherte 


Übereinstimmung der mit Gleich- und Wechselstrom gemessenen Wider- 
standswerte feststellen. Die aus den Gleichstrom-Charakteristiken ermittel- 
ten Widerstände bei + 0, — und + 3 V Vorspannung wurden in Abb. 8 


(Kurve 1) eingetragen. Die Länge 
der senkrechten Striche zeigt da- 
bei ungefährdie Fehlergrenzen an, 
innerhalb derer sich der Wider- 
stand als ctga = R= T aus 
der Charakteristik heraus er- 
mitteln ließ. 

Die Wechselstrommessungen 
an den Zellen der Reihe 3 (mit 
zwei ,,O‘‘-Elektroden) können 
sehr kurz zusammengefaßt wer- 
den: 

1. Die Kapazität bei Be- 
leuchtung ist von der in Dunkel- 
heit gemessenen Kapazität inner- 
halb der Fehlergrenze (< 1 pF) 
nicht unterscheidbar. 


2. Der Wechselstromwiderstand ändert sich nicht mit der Frequenz 
zwischen 300 Hz und 100 kHz und stimmt mit dem Gleichstromwiderstand 


Von 11 ohmschen Kristallen erfüllten 8 Kristalle diese beiden oben ge- 


3. Eine Gleich-Vorspannung bis zu 8 V bleibt ohne Einfluß auf den Wider- 
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gen 4. Es werden keine Formierungsprozesse beobachtet. EN 
Die Punkte 3. und 4. wurden nur bei zwei von den erstgenannten 8 Kri- 
ten stallen erfüllt. 
die Zur Illustration sind in Abb. 9 die Intensitätscharakteristiken innerhalb 
ich von 4 Zehnerpotenzen doppelt-logarithmisch aufgetragen. Die Überein- 
in stimmung der Gleich- erg 
in sung ist besonders bei „ pace 
nd, höheren Intensitäten N ar ne 
er- chungen bei kleineren 1 2 
en stände > 50 k2 zurück- 
> 
der MeBergebnisse 
g. 
te Die Ergebnisse der 
te Messungen an Zellen 0° 
r- mit sperrenden Elek- 41, 9. Intensitätscharakteristik der ohmschen Zellen. — 
el- troden lassen infolge Punkte: Gleichstrommessungen, Karos: Wechselstrom- —__ 
8 derUnbeständigkeit, der messungen bei 100 kHz, Anregung A = 578 my, stärkste 
ge uniibersichtlichen Geo- Intensität 3,7 - 101% Quant/cm? sec x 
a- metrie und der kom- 
n, plizierten Spannungs- und Frequenzabhängigkeit keine einfache quantitative _ 
r- Deutung zu. Ein Vergleich der Meßkurven mit Ergebnissen an Trocken- _ 
gleichrichtern zeigt eine gewisse qualitative Ahnlichkeit, was darauf hindeutet, 
” daß eine Schottkysche Randschicht vorliegt. Das gemeinsame Merkmal 
T- der Scheinwiderstandsdiagramme ist, daß die Abweichung von der Kreis- __ 
form bei hohen Frequenzen eine Erhöhung der Reaktanz X und eine Ver- __ 
N legung des zu erwartenden Bahnwiderstandes den kleineren Werten zu an- 
it deutet. Ähnliche Abweichungen wurden auch bei Trockengleichrichtern 
n beobachtet und durch Schottky derart geklärt, daß sie durch die endliche 
r- Geschwindigkeit der Assoziations- und Dissoziationsreaktionen zwischen - 
freien Ladungsrägern und Störstellen bedingt sind. Nach dieser Auffassung 
e- werden bei steigender Frequenz durch die Trägheit des Auf- und Abbau- 
l- prozesses der Raumladungszone immer weitere Teile des ursprünglich neu- 
r- ralen Halbleiters in die veränderliche Reserverandschicht einbezogen, wo- 
7) durch die Dicke der Randschicht anwächst und die Länge der Strombahn 
im neutralen Halbleiter um denselben Betrag gekürzt wird. Dadurch werden 
2 bei steigender Frequenz sowohl die Kapazität der Randschicht als auch dr 
d Bahnwiderstand kleiner. Es ist wahrscheinlich, daß bei Kadmiumsulfid 
auch diese Vorgänge das Verhalten der Randschicht beeinflussen. Bei Grund- Ee 
2 gitteranregung kommt noch hinzu, daß die Elektronen nur durch Diffu- 
sion unter die Elektroden gelangen können, da die Lichtstreuung im Kristall 
re durch Absorption verhindert wird, so daB sich die en nur unvoll- — 


nen bei der Durchführung der Arbeit und bei der Auswertung der Messungen. 


Be ceabi bane Bei der Redaktion eingegangen am 26. Dezember 1955. 
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Die eigentliche Kapazität der Sperrschicht kann man aus dem Schein- 
widerstandsdiagramm nur grob abschätzen. Wenn man die Formeln anwendet, 
die für einen idealen Gleichrichter mit einem kreisférmigen Scheinwider- 
standsdiagramm gültig sind, erhält man Werte der Randschichtkapazität 
von einigen hundert pF. Wenn man weiter annimmt, daß sich die Randschicht 
infolge der Lichtstreuung und Diffusion der Elektronen unter der Elektrode 
auf einer Breite von etwa 0,2 bis 0,5 mm ausbildet, kommt man zu vernünftigen 
Schichtdicken von einigen 10° em. 

Die Messungen an Zellen mit gereinigten Kontakten besagen folgendes: 

1. Die Unabhängigkeit des Zellenwiderstandes von der Frequenz zwischen 
0 und 100 kHz beweist, daß sogar bei der oberen Grenzfrequenz keine wesent- 
lichen Teile des gesamten Widerstandes kapazitiv überbrückt werden. 

2. Der Befund 1. wird weiter dadurch bestätigt, daß die Kapazität der 
Zelle mit der beschriebenen Apparatur nicht meßbar ist, also höchstens 
einige Prozent der Kapazität der Zellen mit ungereinigten Kontakten be- 
tragen kann. 

3. Empfindlicher als die Kapazitätsmessungen ist die Messung der Ab- 
hängigkeit des Widerstandes von der Vorspannung. Eine schwach vom 
Ohmschen Gesetz abweichende Zelle braucht noch keine meßbare Kapazität 
zu haben und trotzdem kann man eine Änderung des differenziellen Wider- 
standes mit der Spannung beobachten. 

Abschließend läßt sich sagen, daß diejenigen Zellen, die nach der Vorbe- 
handlung der Kontaktflächen bei Gleichstrom ‚ohmisch‘ sind, sich auch dem 
Wechselstrom gegenüber wie rein ohmsche Widerstände verhalten. 

Über den eigentlichen Mechanismus der Entstehung des ohmschen Kon- 
taktes geben die vorliegenden Messungen keinen Aufschluß. Es kann daher 
nicht die Frage beantwortet werden, ob die Auffassung vom ‚innigen‘ Kon- 
takt nach Buttler und Muscheid oder ob Überlegungen, die von der Ent- 
stehung einer Oberflächenschicht anderer chemischer Zusammensetzung aus- 
gehen, mehr zutreffen. 


Herrn Professor Möglich danke ich für die freundliche Einladung, an 
seinem Institut zu arbeiten, sowie für sein förderndes Interesse, desgleichen 
den Herren Dr. Buttlerund Dr. Ammon fürihre Ratschläge und Diskussio- 


Berlin-Buch, Institut für Festkörperforschung der Deutschen Aka- 
demie der Wissenschaften. 
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Von Ernst Schmutzer 


Inhaltsiibersicht 


Falle ergibt sich nach Minkowski die Transfor- 
mation für die elektromagnetischen Materialgleichungen ziemlich leicht. _ a 
Komplizierter werden nun die Verhältnisse, wenn man beliebige Medien TE . 
zugrunde legt, d.h. wenn die Materialgrößen Tensorstrukturen annehmen. Es ie 
Die für diesen Fall zuständigen Transformationsbeziehuugen sollen hier 
abgeleitet werden. Anschließend wird dann die Ableitung der Maxwell- 
Gleichungen für beliebige Medien aus einem Variationsprinzip untersucht. 
In diesem Zusammenhang ergibt sich dann eine Unsymmetriebeziehung für 
die Materialtensoren hysteresisfreier Medien. Der bekannte Schluß auf die Ei 
generelle Symmetrie der Materialtensoren erscheint uns nicht stichhaltig. fe 
Anschließend wird durch mehrere Überlegungen gezeigt, daß nur der Min- oe 
kowski-Tensor der richtigen Abtrennung von Feld und Medium zu ent- 


sprec hen scheint. 


Die in den letzten Jahren durchgeführten Untersuchungen der relativisti- 
schen phänomenologischen Elektrodynamik befaßten sich meist mit dem | = A 
noch nicht endgiiltig gelésten Problem der exakten Form des elektromagne- Sa ee 
tischen Energietensors. Mit einer anderen Frage, nämlich der Formulierung Pr a 
der elektromagnetischen Materialgleichungen anisotroper Medien für beliebige 
zulässige Bezugssysteme befaßte sich kürzlich Marx"). Diese Behandlung as 
des Problems beruht im wesentlichen auf der Verallgemeinerung der vier 

elektromagnetischen Dreier-Feldvektoren auf Vierer-Feldvektoren, wodurch 
sich dann die Materialbeziehungen in Viererform schreiben lassen, wenn man 

auch die elektromagnetischen Dreier-Materialtensoren %;x auf 
tensoren &,;,, %,, verallgemeinert (die 4-dimensionalen Materialtensoren sollen ae? 
durch einen Querstrich gekennzeichnet werden). Man erhält die gewünschten 3 
Relationen dann ohne großen Aufwand für andere Bezugssysteme wegen pe x 
der tensoriellen Formulierung dieser Gesetze. Da nun die Dielektrizität 2 
und Permeabilität im Ruhsystem des Körpers feste Materialeigenschaften 
sind, so soll hier greased werden, wie die Materialgleichungen sich fiir ein a 


Translationsgeschwindigkeit ausdrücken lassen. Wir werden hier bei der Bi 
Ableitung dieser Transformationsformeln einen übersichtlichen Weg gehen. ee “a 3 
Bei den folgenden Betrachtungen wollen wir uns immer auf beliebige an- 
isotrope Medien, deren Materialgrößen Feldabhangigkeit, Hysteresis | sowie 


1) G. Marx, Acta physica Hung. 8, 75 (1954). inb 2 
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Inhomogenität zeigen können, beziehen, sofern nichts anderes gesagt wird, 


im Gegensatz zu Marx, der sich von vornherein auf symmetrische Material. 
tensoren beschränkt. 


Wichtigste Bezeichnungen: 
Die verwendeten Koordinaten seien: 2, = x, y, %3 = 2, 2 = ict, 
so daß wir nicht zwischen kovarianten und kontravarianten Tensorkompo- 


nenten zu unterscheiden brauchen. ar Ra 


A, Vierervektor des elektromagnetischen Potentials j 
0 H, —H, —iD, 
—H 0 H, —iD i 
(Fa) = 3 1 En 
hr H, —H, 0 D, regungstensor 
cr 0 B, —B, —ib, 
B, 0 B, —ik, 
ix) = 


elektromagnetischer Feldstärke- 


B, — B, 0 iE, Tensor 
iE, ik, 0 


F,, dualer Erregungstensor ni aus F;, durch Ersetzen 


id 
von durch ib) 
H,, dualer Feldstärke-Tensor (entsteht aus H,, durch Ersetzen 


Mit Hilfe des Levi-Civita-Symbols!) lassen sich diese dualen Ten- 
soren folgendermaßen schreiben: 


Aix = 4 Sines Hrn Fr = 4 

Dasjenige Bezugssystem, in welchem das Medium ruht, sei mit 2°, ein 
dagegen mit der Geschwindigkeit v in positiver x-Richtung bewegtes mit L 
bezeichnet. 

Im Falle isotroper Körper kann man im Anschluß an Minkowski die 
Relationen D° = ¢ ©, B° = u H® (e Dielektrizitätskonstante, u Permeabili- 


tät) in die folgende kovariante Schreibweise übertragen, um zu 2 übergehen 
zu können: 


ania se (1) 


1 
Poe =eH,, % == = ae 
Diese beiden Relationen lassen sich bekanntlich zur folgenden Beziehung zu- 
sammenfassen: 


F 


H 
> + (eu — 1) Hea — % (2) 
Ähnlich überträgt man das Ohmsche Gesetz j° = a &© auf das System 2: 


rob 3) 
sod 


(j, Vierervektor der Stromdichte, o Leitfähigkeit). 


Ju Ay; 
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§ 2. Transformation der elektromagnetischen Materialgleichungen 


Da auch im Falle isotroper Medien, abgesehen vom System 2°, die Ma- 
terialgrößen Tensorstruktur haben, soerscheint es sinnvoll, die Transformations- _ 
gesetze für allgemeine Materialtensoren abzuleiten. Das betrachtete Medium 
. sei durch den Dielektrizitätstensor e,,, den reziproken Permeabilitätstensor i 

%;, und den Leitfähigkeitstensor o,, charakterisiert. Die Verallgemeinerung oa 
von (1) lautet dann für 2: er 


Far = Ena Hay Fy» = Has (4) 
Die Gewinnung des Zusammenhangs von ®, §, E und ® für Z ist jetzt wesent- _ 
lich langwieriger als für isotrope Körper. Man muß nämlich (4) nach den = 
Komponenten von ® und § auflösen, wobei man aber außerdem die Tensoren i: 
&,, und %,, von 2° nach transformieren muß, da diese als reine Material- __ 
größen ja nur in 2° auftreten. Aus Gründen der in $ 3 durchgeführten ele- ae 
ganten Formulierung des Variationsproblems fir die Maxwell- 
beliebiger Medien gehen wir hier einen anderen Weg, der die gewiinschten 
Beziehungen liefert. 


durch den 4stufigen Materialtensor wobei also g auch 
noch feldstärkeabhängig sein kann: 
Daß sich dann F,,, = f,, (H,,) als eindeutige Funktion von H, ‚tatsächlich — Fr 
in der Form (5) Senn darstellen läßt, kann man durch Reihenentwicklung 4 
der F,, nach den H,, zeigen, wobei vorauszusetzen ist, daß die notwendigen. ee on 
1- mathematischen Forderungen für diese Reihenentwicklung erfüllt sind: N 


Fr = + Hor Hin + Yurorasiı Hor Hu Ha + 
=44,,..H 

=} uvikttixr 

mit 


n 
4 Megs 
Auvik = Aurik + Burorik Hor + Yurorasik Hor Hie + (7) 

e 4 
Dabei lassen sich die Entwicklungskoeffizienten in bekannter Weise durch die _ 
Ableitungen von F,, nach H,, ausdrücken. Die Entwicklungskoeffizienten, 
und damit die A,,.3 haben also Tensorstruktur. Somit brauchen wir uns 

) nicht auf homogene Funktionen zu beschränken). 5 


Wir merken hier noch an, daß Marx!) für symmetrische Materialtensorern _ 
in unsere Bezeichnungsweise übertragen, die folgende Relation für Ay. 
erhält: 

98 Oikrs B— 

= (Gera — % Sia). — 1 (9: — + 


deg 
wobei gz = +. ist und €;, und x;, die 4-dimensionale Verallgemeinerung der + 


ird, 
ial 
et 


w 


Ga me 2) M. v. Laue, Die Relativitätstheorie 1. Bd. Braunschweig 1952, S. 134. 
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Die Antimetrie der Feldtensoren führt zu den Relationen): | 


Auvap Auvaß Auvßa . (8) 


Außerdem gilt für Medien, deren Materialtensoren feldunabhängig sind, die 
Symmetrie der Materialtensoren, wie wir später noch zeigen werden, was 
sich in der folgenden Relation ausdrückt: 


Die folgenden Betrachtungen dieses Paragraphen mögen sich auf Medien 
beziehen, für die (9) gilt. Wegen (8) u. (9) besitzt dann A ‚as (36 — 6)/2 + 6 = 21 
unabhängige Komponenten. Gl. (5) lautet dann ausgeschrieben: 
Dy = Aygıa Ey + + Bs + (Arges By — Ba + Aysıa Bs) 
Dy = Aggıa Ey + + Agısa Es + (Ages By — Bz + Bs) 

Dg = By + Ey + + By — Ba + Asura Bs) 

A, = i(— Ey — Ey — E's) + By — Agsis Bz + Bs 
 (Ayzıa + + Ey) — Arges By + Aisis Be — Bs 

Hy = (— Ey — — + By — Ba + Aare 

Es besteht also zwischen ¢;,, x;, und A,5,„, folgender Zusammenhang: 

= Aya E12 = Araaw = 22 = = &33 = 

= %12 = — %13 = %22 = %23 = — %33 = 
Für das Ruhsystem 2° gilt A,a., = 0, wenn nur ein Index den Wert 4 annimmt, 
da in 2° erfahrungsgemäß ein magnetisches bzw. elektrisches Feld zu keiner 


elektrischen bzw. magnetischen Polaristation Anlaß gibt. Die Transforma- 
tion der A,3,s von 2° in XY erfolgt nun nach der bekannten Formel: 


(11) 


A 


vermöge der Lorentz-Transformation: 
= y (4 + x4), % ay + 3%); 
ke 1 v (13) 
Bei Benutzung der vereinfachenden Relationen (8) und (9) ergibt sich nach 
einiger Rechnung: 


= = Y Ela» £13 = Y Ean = — PR = (es + B? x23), 

Egg = (€33 — PR 29), = Ka = Y = Y 

= — BR Ego), = Y* (x23 + X33 = y? (x33 — (14) 
= 9, = — Y = Anges = IBY 


Assis = — €533); = — (ely + 
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Die gesuchten Beziehungen analog (10) nehmen also die folgende Form*) an: 
D, = Ey + y By + Es) — By (els B, — B,), 
Dy = Ey + By tes Ey) — B? (x33 — Bp 
— By? By + Bs) + (€23 By — Bs)], 
D, = y E, + [(egg Hy + Ey) — B? — x93 E,)] + By B 
+ By? By + B,) + (223 Bs — 33 
H, = 4 By + y (x42 By + m3 Bs) + By (x13 By — x42 Es), 
H, = y Bı + Y? By + 33 By) — B? (e33 By — 25 B)] + By eis 
+ By? Ey + £33 Ey) + Ey — E;)], 


(15) 


Hs = y By + Ba + x33 Bg) — B* By — B,)] —By ; 


— By? [legs Ey + By) + (x95 Bs — Ey)]. 
Man rt sich leicht, daß im Falle isotroper Körper (ei; = &6,,; «fe = 


x»d,= = 1, 2, 3) die folgenden Relationen resultieren: 


F 


\ 


Ähnlich könnte man nun me vorgehen, falls die Leitfähigkeit Tensorstruktur 
besitzt. Die Verallgemeinerung des Ohmschen Gesetzes würde dann lauten: 


Sy» 

q; H,;. 
Wir wollen aber hier auf die Behandlung dieses Falles verzichten. Wir be- 
merken lediglich, daß zum Problem der richtigen Formulierung des Ohm- 
schen Gesetzes in letzter Zeit eine interessante Arbeit von Schlomka?) 
erschienen ist. 


$3. Variationsprinzip für die Maxwell-Gleichungen beliebiger Medien 
Wir wählen hier ähnliche Bezeichnungen wie in einer früheren auf das 
Vakuum bezogene Arbeit*). Man geht in der Elektrodynamik des Vakuums 
bei der Aufstellung des Variationsprinzips bekanntlich aus von der folgenden 
Lagrangedichte: 
1 1 
A=—4,8,—-7 Ay Ain (Fix = (17) 
Deshalb ist offensichtlich im Falle beliebiger Medien die einfachste Verall- 
gemeinerung in der folgenden Form naheliegend°)®): 


4 (18) 


*) Herr Dr. M. Strauss teilte mir inzwischen mit, daß er das Problem durch Auf- 
bau von Aikaß aus Erzeugungsfunktionen zu lösen versuchte. Pyke 3 


s) T. Schlomka, Ann. Physik 8, 246 (1951). ER Be 


‘) E. Schmutzer, Z. Physik 143, 479 (1955). 
5) E. Henschke, Ann. Physik 40, 887 (1913). se 
*) J. Ishiwara, Ann. Physik, 42, 986 (1913). 
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Sommerfeld‘) ner dagegen, in unsere Schreibweise übersetzt, an Stelle 
von 1/, den Faktor !/,, wobei dann allerdings bei der Variation zwar H,, = 
> Su "a — A,, variiert, jedoch F,, konstant gehalten wird. Dies ist nach unserer 
we. Meinung infolge der Verknüpfungsgleichung (5) nicht korrekt, da dadurch eine 
7 es zwischen Medium und Vakuum entstehen würde. Man er. 
kennt nun den Vorzug der Beziehung (5) gegenüber der Formulierung (4), 
da in (5) die Auflösung des Gleichungssystems (4) nach den F;, bereits vorliegt, 
so daß vermöge (5) eine elegante Gewinnung der Maxwell-Gleichungen 
möglich ist, da sich die F,, sofort durch die Ableitungen der A, ausdrücken 
lassen: 


A == A, Hap (19 
is om Wir wollen nun sehen, was dieser Ansatz (18) bzw. (19) im Falle beliebiger 


2 ö: hysteresisfreier Medien zu leisten vermag. Durch Differentiation nach A, 
E und A,; erhält man aus (19): 


Wir fügen nun an dieser Stelle die später noch mies Unsymmetrie- 


Diese Beziehung gilt als Tensorbeziehung in allen im Rahmen der apes 
Relativitätstheorie zulässigen Koordinatensystemen. Für das a 
insbesondere lautet diese Relation, angewandt auf e;; und x), 


) 0 0 de; 
| 
ax? ae (Indizes laufen von 1—3). (22) 


(Man kann dies auch direkt in Dreier-Schreibweise durch die fiir diesen Fall 
analoge Integrabilitätsbedingung zeigen.) 

Aus (22) erkennt man nun, in welcher Weise die Feldabhängigkeit der 
Materialtensoren die Unsymmetrie dieser verursacht. Im Falle feldunab- 
hängiger Materialtensoren folgt aus (22) die Symmetrie dieser, während um- 
gekehrt im Falle symmetrischer Materialien (insbesondere isotroper Medien) 
die Feldunabhängigkeit nicht folgt. 

Mit Hilfe von (21) nimmt nun (20b) die folgende Gestalt an, wenn man 
die geeignet umbenennt: 


7) A. Sommerfeld, Elektrodynamik, Leipzig 1949,S. 281. 4 
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Aus (23) erkennt man nun die Leistungsfähigkeit des Ansatzes (18). 
man nämlich an der üblichen Form der Maxwell-Gleichungen 

OF 


fest, so folgen diese aus der Lagrange-Gleichung: 


wenn das Material den — vier Bedingungen geniigt: rend 
Chains 
(Ha Ay, oH, )- 0, 


insbesondere also, wenn man feldunabhängige Materialien voraussetzt. Damit e = 
ist für solehe Materialien auch im anisotropen Falle (gewöhnlich beschränkt 
man sich auf isotrope Körper) die gestellte Aufgabe gelöst. Nur unter dieser oa “4 
Bedingung ist also der Ansatz (18) fiir die Lagrangedichte mit der Form (24) 
der Maxwell-Gleichungen verträglich. Da aber für beliebige Medien die 
Voraussetzung (26) wohl kaum erfüllt ist, so fragt es sich, wie man in desem 
Falle (18) abzuändern hat, wenn man an der Richtigkeit der einfachen Form ' je 
der Maxwell-Gleichungen (24) festhalten will. 2 

Man sieht leicht ein, daß man dann in diesem Falle nicht von der La- : 
grangedichte direkt, sondern von deren Differential ausgehen muß, so daß 4 
man an Stelle von (18) anzusetzen ee 


Dieses Tensorintegral gi er am den Vektor A, und den Tensor 
H,, analog wie gewöhnliche Kurvenintegrale. Die explizite Form der La- 
grangedichte ist also erst durch Integration zu gewinnen. In (27) tretennun 
eigentlich keine Zahlenfaktoren mehr auf, denn der Faktor !/, verursacht 
lediglich, daß gleiche Produkte nicht doppelt auftreten. Außerdem ist es 
sinnvoll, daß in (27) s, nicht mit unter dem Differential steht. Gl. (18) führt 
nämlich im Gegensatz zu (27) zu einer weiteren Schwierigkeit, da nämlich , 
in Wirklichkeit auch feldabhängig sein kann, also eigentlich beim Differen- BR) 
zieren in (20b) lästige Zusatzglieder liefern würde. In (20b) haben wir nämlich ~ 
die Feldabhängigkeit der s, unterdrückt. 

Aus (27) folgen nun für beliebige Medien sofort mit Hilfe von (25) die Max- © 
well-Gleichungen in der Form (24), denn es ist: ie 

aA 1 
04, 

Damit haben wir für beliebige A inhomogene, feldabhängige, 2 
Hysteresis aufweisende) Materialgrößen die Max well-Gleichungen aus einem 
einfachen Variationsprinzip hergeleitet. Im Falle feldunabhängiger Material- — 
tensoren sind natürlich die Ansätze (18) und (27) äquivalent. 

Man könnte auch daran wee die für das Vakuum gültige Beziehung 
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gegeniiber (27) in einer anderen Form fiir Medien zu verallgemeinern, indem 
man nämlich den Vakuum -Feldanteil unberührt läßt, dafür aber den Einfluß 
des Mediums im Stromanteil eee 

wobei 8; = s, + c—_— und M,, = H,.— F,;, ist. Dies bedeutet, daß man zum 
Leitungs- und einen s, den elektrischen und magnetischen Po- 


larisationsstrom Se hinzuzufiigen hat, ahnlich wie dies durch die Dirac- 


sche Theorie des Elektrons nahe gelegt wird. 

Es soll nun gezeigt werden, daß diese Formulierung des Variationsprinzips 
mit (28) äquivalent ist, da sich der Unterschied zwischen beiden nur in einem 
Variationsintegral über eine Divergenz äußert. Es ist nämlich: 


J= Sf‘ (aA, dH,,) do = 5 Hin dH) do + J, 


wobei f das 4-dimensionale Volumintegral und f das Integral iiber die 
4 

Tensoren A, und H,, bedeutet und außerdem folgende Abkürzung eingeführt 

wurde: 

| | = | J (Mix dA) do — | J Mix (aA) do. 
4 € 4 4 6 

Auf der rechten Seite dieser Gleichung steht ein Integral über eine Divergenz» 


welches sich in ein verschwindendes Oberflichenintegral umwandeln lassen 
soll, so daß übrig bleibt: 


Somit entsteht der gewiinschte Ausdruck: 
[JU dA, — 5 Fy, dH, | do. 


Wir abe s nun den Beweis fiir die Relation (21) nachzuholen: Die aus den 
Feldgleichungen direkt folgende Energiebilanz ergibt auf der einen Seite 
der Gleichung Austrémung plus zeitliche Anderung der Energiedichte und 
auf der anderen Seite neben der Jouleschen Wärme dielektrische Hyste- 
resisanteile. Der bekannte Schluß auf die generelle Symmetrie der Material- 
tensoren entspricht dem Verschwinden der Hysteresisanteile, was aber nicht 
erfüllt zu sein braucht. Es dürften demnach in der Natur auch unsymme- 
trische Medien anzutreffen sein. Medien mit TER TN zeich- 


nen sich dadurch aus, daß der Feldanteil dA? = — — 43 „4H,, von dA in (27) 


kein vollständiges Differential bildet, d.h. A? heme ab von den durchlau- 
fenen Prozessen. Bei hysteresisfreien Medien ist dagegen dA# ein vollstän- 
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Die Integrabilitätsbedingung dafiir: 


Differential. 
vd ‚alas 


OH 
führt nun sofort bei Benutzung von (5) zur Unsymmetriebedingung (21) für 
hysteresisfreie Medien. 


§ 4. Entscheidung zwischen dem Minkowski-Tensor 
und dem Abraham-Tensor 


AbschlieBend sollen nun noch einige Bemerkungen zu der seit etwa einem 
halben Jahrhundert währenden Diskussion über die richtige Form des Energie- 
Tensors gemacht werden. Ohne auf die lange historische Entwicklung dieser 
Diskussion einzugehen, soll ein Kriterium angegeben werden, welches auf der 
Grundlage des Experiments zuungunsten des Abraham-Tensors ausfällt. 

Es werde eine im isotropen Medium in positive x-Richtung fortschreitende 
elektromagnetische Welle betrachtet, für die somit gilt: re 


Ho = 0, E2 =0. (30) 


Die Energieströmungsgeschwindigkeit: V = läßt sich nun leicht 
14 44 
bei Beachtung der ae Transformationsformeln für einen zunächs :hst be- 
liebigen Energietensor E,, 


E,y= 7" ip + BP + Ey + iB 
By = By B + Eu) + 
ausdrücken durch diejenige in 2°, welche mit V, bezeichnet werden möge: 


Ei Eu 


pe Eu rs) 


Diese Formel weist somit für V eine andere Struktur auf als die Formel für 
das Einsteinsche Additionstheorem der Geschwindigkeiten. Im Falle des 
Minkowski-Tensors und des mechanischen Energietensors: ©,, = —u™ 4; % 
entartet (32) in das Einsteinsche Additionstheorem, wie man leicht zeigen 
kann. 


Bei Beachtung von (30) erhält man für z= = —1. Zwischen Ej, und Ej, 


möge nun Proportionalität bestehen in der Form: E}, = « Ej,. Den Faktor « 
miissen wir nun ermitteln. Dadurch entsteht dann aus (32): 


Bewegt sich nun & gegenüber X° mit der Lichtgeschwindigkeit im isotropen 


Medium: v = Vy = °_, so muß für V folgen: V = 0. Dies ist nun tatsäch- 
Eu 


lich der Fall für x = eu. 
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Außerdem läßt sich & auch aus den Fizeauschen Mitführungsversuch 
ermitteln. Bei Vernachlässigung von Effekten höherer als erster Ordnung 
ergibt nämlich die Entwicklung von (33): 
1+a,, 
V = V2) = v(2-— 
Während im Falle des Abraham-Tensors (x = 1) die falsche Formel: 
V=V,—2 v(1 -=) folgt, ergibt sich im Falle des Minkowski-Tensors 


die richtige Mitführung: V = V, — v(1 = 2). Obwohl also der Abraham- 


Tensor einen anschaulich sehr vernünftigen Ausdruck*) für die Kraftdichte 
auf das Medium liefert, die auch eine Kraftwirkung auf die Polarisations- 
ströme enthält, was bekanntlich beim Minkowski-Tensor nicht der Fall 
ist, dürfte damit seine Unzulänglichkeit gezeigt sein. 

Für dieselbe Welle ergibt die Transformation der Energiedichte w (w, Ener- 
giedichte in 2°), bzw. der zugehörigen Photomasse m wegen wdV = 
waV 


n 


w = — + 1] = wor Q 
| 
Die quadratischen Form @ (5) ist nach Abraham (x = 1) positiv definit und 
nach Minkowski (x = eu) indefinit. Für v > — wird also die Energie- 
dichte negativ, eine zwangsläufige Konsequenz des Minkowski-Tensors?). 
Bei nichtabgeschlossenen Systemen dürfte also der Energietensor des 
Teilsystems nicht mit dem metrischen Energietensor (hier Abraham -Ten- 
sor) 77; identisch sein. Analog zum kanonischen Tensor 7; können wir aber 
den Minkowski-Tensor E}; in den La grange-Formalismus einfügen, 
wenn wir an Stelle von A, und A,, nur die physikalisch sinnvollen Größen 
A, und H,, variieren, d. h. 


—Ab,,=Fi, (8%) 


04, 
bilden. Im Falle des Minkowski-Tensors folgt dann für den Spin-Ener- 


gietensor Tj, = Ti, — Ei, = zz, (A; Fax) ein Divergenzausdruck, für die 
o 


Spindichte: ©” und für den Spin: A (w Kreisfrequenz der Welle). 


Herrn Prof. Dr. Falkenhagen danke ich herzlichst für die Möglichkeit, 
an seinem Institut die behandelten Fragen zu untersuchen, sowie für eine 
Reihe wertvoller Hinweise. 


8) G.Marx u. G. Györgyi, Acta physica Hung. 3, 213 (1954). 


” hs Rostock, Institut fiir theoretische Physik der Universitat. 
der Redaktion eingegangen am 14. Dezember 1955. 
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Die untere durch statistische Schwankungseffekte 
bedingte Meßgrenze einer Verstärkeranordnung 


Th. Ankel und W. Wintermeyer 


Inhaltsiibersicht 


Die prinzipiellen Anordnungen fiir phasenrichtige und quadratische 
Gleichrichtung werden hinsichtlich ihrer statistischen Schwankungen be- 
handelt und die dadurch lurch gesetzten unteren Meßgrenzen für Signale bestimmt. 

Ms 


Viele Probleme der modernen Experimentalphysik bedingen, daß Signale 
— in Verstärkeranlagen verstärkt — aus dem Rauschpegel der Verstärker 
herausgehoben werden. In der Praxis hat sich für periodische Signale, die 
im folgenden ausschließlich betrachtet werden, die Methode der phasen- 
richtigen Gleichriehtung besonders bewährt. Für diese Methode wird die 
untere Meßgrenze bestimmt und einer Anordnung mit quadratischer Gleich- 
richtung gegeniibergestellt. 
In Abb. 1 sind die beiden miteinander zu vergleichenden Me Meßanordnungen 
im Blockschema dargestellt. 


Fit) Ft) 


> >- 
> 


a 


R(w) 


b 
Abb. 1. Blockschema einer Anordnung mit a) phasenrichtiger, b) quadratischer Gleich- 


richtung. In der Abbildung bedeuten V einen Verstärker, 7’ einen Tiefpaß, ph.Gl.einen 


phasenrichtigen und gua. Gl.einen quadratischen Gleichrichter 


Die Meßgröße tritt in den Verstärker ein, dessen Übertragungseigen- 


schaften durch die komplexe Funktion P(w) gekennzeichnet sind. Der aus : : 


dem Verstärker austretende Strom wird im einen Falle phasenrichtig gleich- 


gerichtet, im anderen Fall durch den Gleichrichter quadratisch umgeformt _ 


und über einen Tiefpaß mit der Ubertragungsfunktion R(w) auf das Anzeige- 
instrument übertragen. Die Frequenzcharakteristik des Instrumentes wird 
durch die Funktion R(w) mitberücksichtigt. 
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Die störenden Rauscherscheinungen, die vom Schrot- und Johnson- 
effekt am Verstärkereingang verursacht werden, bewirken eine Schwankung 
des Ausschlages am Anzeigeinstrument. Die Größe dieser Schwankung 
wird durch das mittlere Schwankungsquadrat charakterisiert, das durch den 
Ausdruck 


(5 = (Dit) — = BH — OW 


gegeben ist, wobei ®(t) den Momentanwert der Anzeige und ®(t) deren zeit- 
lichen Mittelwert bedeuten. 

Ein Strom oder eine Spannung am Verstärkereingang soll dann als nicht 
mehr meßbar bezeichnet werden, wenn der durch die Meßgröße verursachte 
Instrumentausschlag im Schwanken der Instrumentanzeige untergeht, d.h. 
wenn 


Dyiess s | (2) 


ist. Dyes, bedeutet dabei den Ausschlag, der allein durch den Meßstrom oder 
die Meßspannung erzeugt wird. 

Zur Bestimmung der unteren Meßgrenze der Anordnungen wird zunächst 
das mittlere Schwankungsquadrat des allein von der Rauscherscheinung 
verursachten Ausschlages des Anzeigeinstrumentes bestimmt. Der Ver- 
gleich mit dem allein durch die Meßgröße hervorgerufenen Ausschlag liefert 
dann die untere Nachweisgrenze der Verstärkeranordnung. 


1. Phasenrichtige Gleichrichtung 


Zur Durchführung wird eine Methodik in Anlehnung an eine Arbeit von 
H. Dänzer!) verwandt. Der durch den Rauscheffekt verursachte Strom- 
bzw. Spannungsverlauf am Eingang des Verstärkers wird in einzelne Impulse 

von der Zeitdauer t zerlegt, die, für sich 


= allein betrachtet, am Verstärkerausgang 
= K—“——-~._ einen Stromverlauf / (t,) erzeugen, wenn 
1 +, die Zeit bedeutet, die der betreffende 

i a r + Impuls vor der MeBzeit t zeitlich zurück- 


liegt. Die Zeitachse hat man sich dabei 
Abb. 2. Aufteilung der Zeitachse in kleine Zeitabschnitte von der Dauer r 
unterteilt zu denken, so daß t, = k-t 
gesetzt werden kann, wobei k den Index des betrachteten Zeitabschnittes 
bedeutet (vgl. Abb. 2). 
Im Falle des Schroteffektes entspricht der betrachtete Impuls einer 
Elektrizitätsmenge q,, im Falle des Johnsoneffektes ist der Impuls durch 
das Integral 


Ik “a €, dt 


gegeben. 


Nach!) erhält man für den aus dem Verstärker austretenden Strom zur 


Meßzeit t als Überlagerungseffekt der von den einzelnen Zeitelementen r 
1) H. Danzer, Zur Theorie des Schrot- und Johnsoneffektes, Ann. Physik (6) 8, 
176 (1960). 
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= Qn J P (w) e% do, 
(1) 
eit- wobei P(w) die Frequenzcharakteristik des Verstärkers bedeutet (vgl. Abb. 1). 
Die phasenrichtige Gleichrichtung, die experimentell auf verschiedene 
cht Weise verwirklicht werden kann, läßt sich im Prinzip durch Multiplikation 
hte des Verstärkerausgangsstromes /(t) mit einer zweiten Größe darstellen, die 
Jh in Frequenz und Phase mit dem Meßstrom übereinstimmt. Damit ergibt sich 
für den Strom J, (t) hinter dem Gleichrichter: 
(2) I(t) = «+ I(t) sin + y) = % = P (t,) sin (Wy t + y). (5) 
‚der Dabei bedeuten « einen Proportionalitätsfaktor und y die Phase des Meß- 
stromes. 
shst Um die Wirkung des Tiefpasses zu beriicksichtigen, zerlegt man jedes 
ung einzelne Glied der Summe aus (5) nach Fourier und multipliziert die ein- 
Jer- zelnen Fourierkomponenten mit der Ubertragungsfunktion R(w) des Tief- 
fert passes. 
Zur Fourierdarstellung von (5) ist es zweckmäßig, die Zeit tin der Sinus- 
funktion unter Bezugnahme auf Abb. 2 zu ersetzen durch 
t=t + &, 
von wobei t; ebenfalls von dem Laufindex des betreffenden Zeitintervalles t 
abhängt. Unter der Fourieridentitat k 
ang 
a erhält man für die Fourierdarstellung eines Summengliedes aus (5) 
nde 
ck. sin (ot + y) 
, 
er T f [ei (ote?) T(w . P (w + @p)] dw. (7) 
ttes Aus (5) und (7) ergibt sich fiir den Ausschlag des Anzeigeinstrumentes, der 
1) dem aus dem Tiefpaß austretenden Strom proportional sein muß, wobei ß 
ner einen Proportionalitätsfaktor bedeutet: 
ırch 
5 +00 
= f P (@ — @) : R(w) dw 
nit 


e 


> 
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geht Gl. (8) über in 
k 


de bet Zur Bibung des mittleren Schwankungsquadrates müssen nach (1) die 
Größen und bestimmt werden. 


Be Für ® und damit auch für ®® ergibt sich, wenn man die Mittelung über 
| a Be Zeit t nach dem in!) angegebenen Mittelungsverfahren durchführt, der 
Wert Null. 
i oe Zur Bestimmung von ®2 muß der Ausdruck (10) mit seinem konjugierten 
Wert multipliziert und anschließend über t gemittelt werden. 
Man erhält für das mittlere Schwankungsquadrat der Rauscherscheinung 


“ 


Bedenkt man, daß t, = k-r ist, so kann bei kleinem r die Summe über k 
durch ein Integral über die Zeit k-r = T ersetzt werden, wobei das Zeit- 
differential dT den Wert r annimmt. Damit folgt aus (11): 


‘hea eine einfache Integralumformung geht (12) über in: MR a 


2.0p 2411p 2\ | 
8x (wm) (@ — Wp) | + | (@ + |”) (13) 
+00 

#4 (1) 

Gl. (14) stellt das mittlere Schwankungsquadrat der Rauscherscheinung bei 
einer Meßanordnung mit phasenrichtiger Gleichrichtung dar. 

Für die komplexen Übertragungsfunktionen P (w) und R (w) gilt: 
P (-o) = P*(+ 0), R(—w) = R*(+o) und damit auch |P (— o)?= 
IP(+ @)|%, |R(-o)|®=|R (+ 


Im Falle des Schroteffektes kann nach!) für qj/t gesetzt werden: 
qlt=e:i, 
wobei e die Elementarladung und i den mit Schroteffekt behafteten Gleich- 


strom bedeuten. Die Funktion P (w) stellt den Zusammenhang zwischen 
dem Eingangs- und dem Ausgangsstrom des Verstiirkers dr. 7 


4, 


Für den Johnsoneffekt gilt: ai &: 
= 2-k-T- W, 


dabei bedeuten k die Boltzmann-Konstante, T die absolute Temperatur 
und W den Ohmschen Widerstand am Verstärkereingang, der die Johnson- 
spannung erzeugt. P(w) gibt den Zusammenhang zwischen der Rausch- 

' spannung am Verstärkereingang und dem Verstärkerausgangsstrom an. 
Zur Bestimmung der unteren Meßgrenze wird der Ausschlag des Instru- 


u 7 Se mentes berechnet, der durch die Meßgröße allein hervorgerufen wird. Nimmt 
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man die zu messende Wechselspannung beziehungsweise den Wechselstrom — 
von der Form 
(10) u = Up: SIN Wet 


1) die an, so ergibt sich hinter dem phasenrichtigen Gleichrichter nach (5) ein 
Strom J, (t): 
über I, (t) = |P (@p)| sin? (wp t + y). 
t, der Bedenkt man, daß in der Praxis der Durchlaßbereich des Tiefpasses fast  —_ 
immer auf Frequenzen von Null bis zu einer Frequenz kleiner als 2 wo be 
ierten schränkt ist, so kommt von J, (t) nur die Gleichstromkomponente zur Wirkung. 
Man erhält für den poate des Anzeigeinstrumentes: 


(t) = P |-B-4-| R(0)). (16) 
Aus (2) folgt mit (14) und (15) für die kleinste noch nachweisbare Meßgröße: en 


inung 


(11) 


ra 1 2.] 2 
Diese Beziehung ist allgemein giiltig. es 
Idealisiert man die Übertragungsfunktionen von Verstärker und TiefpaB 
(12) als Rechtecke mit den Frequenzbreiten 2 Aw, für P (mw) und Aw, für R (w), 
so läßt sich für diesen speziellen, für die Praxis meist ausreichenden Fall das 
Integral in (16) leicht auswerten. Liegt insbesondere der Durchlaßbereeich 
von P (w) symmetrisch um die Frequenz wp, so ergibt P (w — w,) eine Fre- 
quenzverteilung, die symmetrisch um die Frequenz Null liegt. Das Integral _ 
wird also durch die jeweils kleinste Frequenzbreite bestimmt und ist dieser 
proportional. Man erhält für die untere Meßgrenze: | 
wenn Aw die jeweils kleinste der beiden Bandbreiten Aw, und Aw, von Ver- — 
stärker und Tiefpaß bedeutet. ; 


2. Quadratische Gleichrichtung 


Bei der quadratischen Gleichrichtung ist zu berücksichtigen, daß die 2 
Rauscherscheinung infolge der Quadrierung eine 
erhält, so daß der Ausschlag des Instrumentes um einen konstanten Wert 
Schwankungen ausführt. Das mittlere Schwankungsquadrat der störenden 
Rauscherscheinung ist nach Gl. (1) gegeben: 


(6 () =H — Oy. 


Die mathematische Behandlung geht auf die Gln. (3) und (4) zurück. 
ih Aus dem quadratischen Gleichrichter tritt ein Strom J, (t) 


I, (t) = y Pi) 


isch- 


aus, mit y als Proportionalitätsfaktor, der für die Schwankungseffekte in der 


Form 
stru- 


nmt 2 (t) Te Ue (ty) P* (ty) 


— 

eich- 

chen 
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dargestellt wird. Die Wirkung des Tiefpasses läßt sich wieder in der Weise 
berücksichtigen, daß jeder einzelne Summand von (18) nach Fourier zer 
legt wird und die einzelnen Fourierkomponenten mit R (w) multipliziert 
werden. 

Zur Fourierzerlegung des Produktes, das von den beiden Zeiten t, und t, 
abhängt, wird eine der beiden Zeiten, z. B. t,, ersetzt durch 


t, = t, + At,,- (siehe Abb. 3). 


Man führt also die Zeitskala für t,- durch eine lineare Transformation in die 
für t, über. Die Größe At,,, die nur von den Indizes k und k’ abhängt, ist 
für die Fourierzerlegung eine Kon- 


OtKK: tin stante, die einen konstanten Phasen- 
. 
— faktor ergibt. 
t+=0 t In Abb. 3 ist die Zeit t, willkürlich 
größer gewählt als t,; die Überlegung 
: Ig' gilt aber auch für ¢, kleiner als t,, wo- 
Abb. 3. Aufteilung der Zeitachse bei sich nur das Vorzeichen von Ahy 
2 ändert. 
Für die Fourierdarstellung erhält man leicht: 
Phar 
+00 


J P* Q) P(w dd. (20) 


Der e-Faktor in (20) stellt keine Zeitabhängigkeit der Fouriertransfor- 
mierten B,, (w) dar, sondern bedeutet — wie schon erwähnt — einen kon- 
stanten Faktor, der nur von den beiden Indizes k und k’ abhängt. Für ver- 
“ hiedene k und k’ ergeben sich also auch verschiedene Amplitudenfunktionen 


x (w). : 
Mit (19) folgt aus (18) für den Ausschlag des Instrumentes: = = 


Taw 5% Buy (0) R(w) ekdo, (2) 
kk’ — 00 


wobei y’ ein Pröportionalitätsfaktor ist, oder mit der Definition os 
» 
hy (4) = Qn J By (w) R (w) e'®'x dw (22) 
auch 


Zur Berechnung des mittleren Schwankungsquadrates werden nach Gl. (1) 
wieder die beiden Größen © und ®2 bestimmt. Für ® (t)” ergibt sich mit (23) 


Dt)” =y’? = Fer (ty) fii (tr) (24) 


und für ®2 (t) 4 


(25) 
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Setzt man =kr=T, t,=1-t=T" und At, =(-Wr=T-T, 


= y Sf IP (Q)2- |P(Q |B dQdeo. (29) 
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Die Mittelung erfolgt analog zum Vorhergehenden bei festem k, k’, 1 und I’ 

fir jedes einzelne Glied der vierfachen Summe. Wegen der statistischen 

Verteilung der einzelnen Impulse bleiben von der Summe nur die Glieder _ 

übrig, bei denen entweder alle vier oder jeweils zwei Indizes gleich sind. Da- 
mit folgt: 

(t) | (ty) |? + ye de 


(tx) fii (t,) fer ) fkı (t,) + fet fi 


Die Mittelwerte der Schwankungsgrößen bestimmen sich nach der Poisson- Er ! 
verteilung näherungsweise zu: pes 


4=3 (4); eat = (Gi)”. (27) 
Faßt man die erste Summe von (26) entsprechend der Beziehung (27) jeweils — 
mit einer der drei übrigen Summen zusammen, so entstehen.drei vollständige 
Doppelsummen über k und /, die auch die Glieder mit k = / enthalten. Mit 
(24), (26) und (27) erhält man nach (1) für das mittlere Schwankungsquadrat: 


so lassen sich analog zum Vorhergehenden die Doppelsummen über k und/] | 
als zeitliche Doppelintegrale schreiben, wobei die Zeitdifferentiale dT und ==> 
aT’ jeweils wieder den Wert r annehmen. Die Auswertung dieser Doppel- _ ; 
integrale ist im Anhang durchgeführt und ergibt für das mittlere Schwan- _ 
kungsquadrat der Rauscherscheinung bei einer Meßanordnung mit qua- — 
ther Gleichrichtung: 


Der Instrumentausschlag, der durch die Meßgröße 
allein hervorgerufen wird, hat in Analogie zu den bei der phasenrichtigen 
Gleichrichtung bereits Re Uberlegungen den Wert: 
() = y' ub |P (cp)? - |R(0)}. (30) 
Für die untere ag erhält man aus (29) und (30) unter Benutzung 
von Gl. (2): 


/% Va 


(31) 
ip | PQ anneal 


Auch diese allgemein giiltige Beziehung soll fiir den Fall rechteckiger E23 2 
Frequenzcharakteristiken ausgewertet werden. Nimmt man wieder ls 
Frequenzbreiten 2 Aw, für P (w) und Aw, für R (w) an, so folgt für: 


Weise 
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a) Ads < 2 Aa, 
und speziell für Aw, < 2 Aa, AN 
Homin = - i 32a) 
Upmin = | * (33) 


Ein Vergleich der phasenrichtigen Gleichrichtung mit der quadratischen 
in bezug auf die untere Meßgrenze führt im Fall rechteckiger Frequenz- 
charakteristiken zu folgenden Ergebnissen ?): 

a) Aw, < 2 

Für den Fall, der auch in der Praxis fast immer vorliegt, daß die Frequenz- 

breite des Tiefpasses klein ist gegenüber der des Verstärkers, erscheint die _ 
phasenrichtige Gleichrichtung günstiger als die quadratische. Bei der 

quadratischen Gleichrichtung hängt die untere Meßgrenze noch von der 

Bandbreite des Verstärkers ab, wie aus Gl. (32a) hervorgeht. Nimmt man 

beipielsweise an, daß ein handelsüblicher Verstärker mit einer Bandbreite bis 

zu 10 kHz Verwendung findet, so unterscheidet sich die untere Meßgrenze 

von der bei phasenrichtiger Gleichrichtung um den Faktor 10, wie ein Ver- 

gleich von (17) und (32a) ergibt. 

b) do, > 2 Aa, . 

Besitzt der Verstärker die kleinere Bandbreite, so sind beide Meßanord- 
nungen einander gleichwertig. In diesem Fall hängt die untere Meßgrenze bei 
beiden Methoden nur von der Frequenzcharakteristik des Verstärkers ab. 


Physikalisch bedeutet dies, daß der Tiefpaß keinen Beitrag zur unteren Meß- 
grenze liefert. 


3. Anhang 


Zur Berechnung der Doppelsummen aus Gl. (28) setzt man die Beziehungen 
(20) und (22) in (28) ein. Man erhält für die erste Summe von ( 


if P* (Q) P (m + Q) R 07 a0 da] (34) 


+0 


SS PQ) Pr +2) Be (m 3 


2) Die i in liek Gin. (32) und (33) angegebenen Beziehungen stimmen mit den Resultaten 
einer von H. Müser auf andere Weise durchgeführten Rechnung überein, die anläßlich 
eines Kolloquiums vorgetragen wurde. 
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Führt man die Multiplikation aus, so ergibt sich für die rechte Seite von (34): 


© +00 
(wm + 2) P* (@ + 92’) etia ) etia 


R (o) R* (w') T dQ dQ! dw dw’ aT aT’. 
Durch geeignete Umstellung geht (35) über in: ues ER 


1 co +00 1 co +00 a 
P (Q’) P* (0 + Q') ei? T. | 
.ei@T P* (Q) P(w+ 2): R (w) - R* (w') T dQdwdo' aT. 


Unter Berücksichtigung der Fourieridentität Gl. (6) reduziert sich das 
sechsfache Integral (36) auf ein vierfaches: 


(36) 


oo +00 


\P (Q)|2 P* +2) P (w + Q) R (w) R* (w') 7 


(2 
Eine weitere Umstellung führt auf 


| P (Q) |? P* (w’ + Q) R*(w’) dQ do’. 
Durch nochmalige Anwendung von Gl. (6) erhält man das Ergebnis: 
(ty) fei (ty) taal | P (Q)|? | P (w+ 2)|? | R (co) |? dQ dw. (39) 
— oo : 


Die Auswertung der zweiten Doppelsumme in Gl. (28) führt in analoger Weise 
zu dem gleichen Resultat: 


= SJ |P + MP |R de. (40) 


Frankfurt a. M., Institut für angewandte Physik der Universität und 


2 Bei der Redaktion eingegangen | am 4. Februar 1956. 
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Das Reflexionsvermégen und die Durchlässigkeit 
eines Schichtsystems visko-elastischer Medien bei Einfall 
einer ebenen Schallwelle unter beliebigem Winkel. |') 


Vv on Siegfried Kästner 


Inhaltsübersicht 


In der Arbeit wird zunächst kurz die Wellenausbreitung in einem unbe- 
grenzten visko-elastischen Medium betrachtet. Es zeigt sich hierbei zum 
Unterschied von einem rein-elastischen Medium neben der auftretenden 
Schallabsorption, daß sich im allgemeinen Fall inhomogene Wellen ausbreiten. 
Für den Verschiebungsvektor folgt hieraus, sowohl bei einer Kompressions- 
welle als auch bei einer Scherwelle, eine elliptisch polarisierte Schwingung. 

Anschließend wird die Wellenausbreitung in einem geschichteten Medium 
behandelt, wobei sich unter Verwendung von sechsreihigen Matrizen eine 
übersichtliche Darstellung ermöglicht. 

Hierzu wird zunächst die Wellenmatrix ® einer Schicht aufgestellt, 
welche die Verschiebungs- und Spannungskomponenten in einer zu der Ober- 
fläche parallelen Ebene als Funktion der Wellenamplituden der sich innerhalb 
der Schicht ausbreitenden Kompressions- und Scherwellen liefert. Außerdem 
wird speziell für x = 0 die reziproke Wellenmatrix B-! (0) benötigt, die um- 
gekehrt die Wellenamplituden als Funktion der Verschiebungs- und Span- 
nungskomponenten in der Ebene x = 0 ergibt. Als Produkt dieser beiden 
Matrizen erhält man schließlich die Kettenmatrix m einer Schicht, durch die 
der funktionelle Zusammenhang zwischen den Verschiebungs- und Spannungs 
komponenten der rechten Grenzfläche der Schicht und den entsprechenden 

Größen an der linken Grenzfläche hergestellt wird. Die Kettenmatrix eines 
beliebigen Schichtsystems ergibt sich danach auf Grund der Stetigkeits- 
 bedingung für den Verschiebungsvektor und die Normal- und Schubspan- 
nungen an einer Trennfläche zweier Medien aus dem Produkt der Ketten- 
 matrizen der einzelnen Schichten. 

Ist die Kettenmatrix des Systems bekannt, so lassen sich in einfacher 
Weise die reflektierten und durchgehenden Wellen als Funktion der einfal- 

lenden Welle bestimmen. 

x Zahlreiche Probleme der Akustik lassen sich in idealisierter Darstellung 
als Durchgang einer ebenen Welle durch ein System von parallelen, ebenen 


1) Teilergebnisse einer an der math.-nat. Fakultät der Humboldt-Universität Berlin 
1954 ne Dissertation mit gleichem Titel. 
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Schichten behandeln. Dabei sei nur an Probleme der Ultraschall-Meßtechnik 
und -Therapie, der Bau- und Raumakustik sowie der Schalldämmung erin- 
nert. Zur Lösung der verschiedenartigen praktischen Aufgaben in den ein- 
einen Gebieten interessiert bei einem geschichteten Medium vor allem die 
Frage, welcher Anteil der einfallenden Welle reflektiert wird und welcher 
Anteil das Schichtsystem durchsetzt. Von dieser Aufgabe sind bisher drei 
spezielle Fälle gelöst worden, der senkrechte Einfall der Welle, der schräge 
Einfall für eine Schicht und der schräge Einfall für ein Schichtsystem rein- 
elastischer Medien. 


Der senkrechte Einfall der Welle wurde von P. J. Ernst?) und M. Nuovo?) 
behandelt. Eine zusammenfassende Darstellung über Wellenausbreitung 
in geschichteten Medien bei senkrechtem Einfall unter Verwendung von 
Analogien wurde von K. Altenburg und S. Kastner‘) gegeben. 


Komplizierter werden die Verhältnisse bei beliebigem Einfallswinkel der 
Welle. Betrachtet man zunächst speziell eine Schicht, und zwar ein rein- 
elastisches Medium, so hat für diesen Fall bereits Lord Rayleigh®) eine 
Formel für eine schubspannungsfreie Platte angegeben, was physikalisch der 
Durchdringung einer Flüssigkeits- oder Gasschicht entspricht. Den Fall des 
Schalldurchgangs durch eine feste Platte hat in exakter Form H. Reissner®) 
berechnet, wobei seine Endformel allgemein in der von R. Bär”) umgeschrie- 
benen Form verwendet wird. F. Levi und N. S. Nagendra Nath*) haben 
eine andere Ableitung der Reissnerschen Formel gegeben, die der strah- 
lungsakustischen Betrachtung entspricht und dementsprechend anschaulicher 
ist. Den Speziallfall dünner Platten hat L. Cremer?) behandelt und gezeigt, 
daß ein Maximum der Schalldurchlässigkeit eintritt, wenn die Spurgeschwin- 
digkeit der einfallenden Welle mit der Geschwindigkeit einer freien Biegewelle 
der Platte übereinstimmt und diesen Effekt als Koinzidenz bezeichnet. 
Anschließend hat J. Götz!P) gezeigt, daß auch bei dicken Platten die Maxima 
des Schalldurchgangs als Koinzidenzen gedeutet werden können. Schließlich 
hat A. Schoch!!) eine Ableitung der Durchlässigkeit gegeben, die den Zu- 
sammenhang zwischen Totaldurchgängen und freien Plattenwellen unmittelbar 
erkennen läßt. 


Der Fall einer visko-elastischen Platte bei beliebigem Einfallswinkel ist 
von Y. Torikai!2) berechnet worden. 


Bei einem Schichtsystem wurde der schräge Einfall einer Welle bisher 
nur für rein-elastische Verhältnisse behandelt. Die Schalldurchlässigkeit 
eines geschichteten Mediums, das abwechselnd aus zwei verschiedenen Ma- 
terialien besteht, wovon das eine fest und das andere flüssig ist, hat R. B. 

2) P. J. Ernst, J. sci. Instrum. 22, 238 (1945). cae % 

3) M. Nuovo, Ric. scient. 16, 88 (1946). rl: 

‘) K. Altenburg u. 8. Kästner, Ann. Physik (6) 18, 1 (1953). ar 

5) Lord Rayleigh, The theory of sound II, $271, London 1896. __ 

*) H. Reissner, Helv. phys. Acta 11, 140 (1938). at 

’) R. Bär, Helv. phys. Acta 11, 397 (1938). RT 

*) F. Levi u. N. 8. Nagendra Nath, Helv. phys. Acta 11, 408 (193 

*) L. Cremer, Akust. Z. 7, 81 (1942). 

“) J. Götz, Akust. Z. 8, 145 (1948). 

") A. Schoch, Acustica 2, 1 (1952). 

*) Y. Torikai, J. phys. Soc. Japan 8, 234 (1953). 
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Lindsay) berechnet. Der Fall des Schalldurchgangs durch ein System au 
beliebig vielen, verschiedenen Schichten wurde von W. T. Thomson!) 
und Y. Torikai) gelöst. 

Im folgenden soll die Lösung des Problems für den allgemeinen Fall dureh- 
geführt werden, daß eine ebene Schallwelle unter beliebigen Winkel auf ein 
System aus beliebig vielen visko-elastischen Schichten trifft. Das Reflexions- 
vermögen und die Durchlässigkeit des Schichtsystems werden berechnet}), 

Dazu ist erforderlich, zunächst einmal die Wellenausbreitung im unbe- 
grenzten visko-elastischen Medium etwas näher zu betrachten. Zum Unter. 
schied von rein-elastischen Medien zeigt sich hierbei neben der auftretenden 
Schallabsorption eine Inhomogenität der Welle sowie eine elliptisch polari- 
sierte Schwingung des Verschiebungsvektors. = 

Ge 


2. Die Grundgleichungen 


Das dynamische Verhalten eines visko-elastischen Mediums wird durch 
das Newtonsche Kraftgesetz fiir Kontinua und eine Beziehung zwischen 
dem Spannungs- und Deformationstensor beschrieben. Bei Beschränkung 
auf kleine Deformationen, wie sie bei akustischen Wellen auftreten, hat bei 
letzterer das erweiterte Hookesche Gesetz Gültigkeit, bei dem ein linearer 
Zusammenhang zwischen den Spannungen und Verschiebungen sowie deren 
Ableitung nach der Zeit besteht. Unter Verwendung kartesischer Ortskoor- 
dinaten x, y, z und der Zeitkoordinate t lautet das Hookesche Gesetz für 
homogene, isotrope Körper!‘), auf die sich die folgenden Ausführungen be- 
schränken: 

av, a av, bv, 
2 (u, + (4, + A, tye = (it. + Me 5) 
(A, + A, te= (u. + My 5) + (2,1) 


2 (ite + +(4, +4, 5 tev = (tte + ag) (Ge + 


Hierbei werden mit v,, und t,,„ die Komponenten des Verschiebungs- 
vektors » bzw. des Spannungstensors t bezeichnet. Neben den Laméschen } 
Elastizitätskonstanten A,, u, treten die entsprechenden Viskositätskon- 
stanten A,, a, auf. Außerdem wurde die Divergenz u des Verschiebungs- 
vektors » eingeführt (u = div dv). 

Von den Elastizitätskonstanten 2,, u, und den Viskositätskonstanten 4, 
u, kénnen bei einem speziellen Stoff einzelne Elemente den Wert null haben. 
Hierauf läßt sich eine Einteilung der Medien aufbauen, die in Tabelle 1 dar- 
gestellt ist, wobei das Kreuz bedeutet, daß der entsprechende Koeffizient von ! 
verschieden ist. In der Gruppe A sind die Medien zusammengefaßt, in denen 


13) R. B. Lindsay, J. Acoust. Soc. Amer. 11, 178 (1939). digi oitndl tam 

14) W. T. Thomson, J. appl. Phys. 21, 89 (1950). ABBEnEn ad 

15) Y. Torikai, J. Acoust. Soc. Japan 8, 21 (1952). 

16) Damit wird auf Grund des hier gültigen Überlagerungsprinzips die Möglichkeit 
geschaffen, auch kompliziertere Wellenvorgänge in einem System geschichteter Medien, 
wie z. B. zeitlich begrenzte Vorgänge oder den Einfall eines räumlich begrenzten Biindels 
zu behandeln. Fe 

17) Vgl. H. L. Oestreicher, J. Acoust. Soc. Amer. 28, 707 (1951). 
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Tabelle 1 
Einteilung der Medien 


| 4, Me | A, 
Allg. visko-elastisches Medium x x x x 
Al | baie abe 
in- i i h 
Rein-elastisches Medium x | war 
tmn =0 | Ideale Flüssigkeit x | x 
B (m = n) 
Rein-elastische ideale Flüssigkeit x a 


Schubspannungen auftreten können, während bei den Medien der Gruppe B 
sämtliche Schubspannungen verschwinden. Von den 15 Kombinations- 
möglichkeiten haben im Rahmen dieser Arbeit nur die vier praktische Be- 
deutung, die in Tabelle 1 mit Namen aufgeführt sind, doch gelten die folgenden 
Ableitungen ganz allgemein für alle Stoffe. Auf die Tabelle wird bei Betrach- 
tung der Wellenausbreitung in geschichteten Medien noch Bezug genommen 
werden. 

Das Newtonsche Kraftgesetz für Kontinua lautet, falls man, wie hier 
erlaubt ist, die Konvektionsglieder streicht: 


wobei o die Dichte ist. 

I. 


3. Lösung für zeitlich periodische Vorginge $= en. 


Speziell für periodische Vorgänge macht man für die Komponenten des 
Verschiebungsvektors und des Spannungstensors folgenden Ansatz: 


K (Vin e~ tat) = (Da e-iat) 18), (3,1) 


Hierbei kennzeichnet R den Realteil, i = V —1 die imaginäre Einheit, 
o=2nv die Kreisfrequenz und v die Frequenz. 


18) Prinzipiell kann beim Ansatz (3,1) der Exponent positiv oder negativ gewählt 
werden. Die einzelnen Autoren benutzen zum Teil den einen zum Teil den anderen Weg, 
was beim Vergleich der Arbeiten verschiedener Autoren zu beachten ist. In der dieser 
Arbeit zugrunde liegenden Dissertation wurde mit positiven Exponenten angesetzt, 
während hier zum negativen Exponenten übergegangen wurde. 


Ann, Physik. 6. Folge, Bd. 18 13 


aus 
ein 
ons- 
16 
nbe- 
ter- 
den 
lari- 
arer 
ren 
oor: 
für 
be- 
hen 
<on- 
ngs- 
ben. 
dar- 
nen 
ikeit 
dels 
. 


Annalen der Physik. 6. Folge. Band 18. 1956 


Für die komplexen Amplituden der Verschiebungs- und Spannungs- 
komponenten V,,, T,,,, erhält man damit unter Einführung der komplexen 
Materialkonstanten A, u 


fuk —tek 


die Beziehungen 


wobei für die Divergenz des Verschiebungsvektors » ganz analog zu Gl. (3,1) 
u = KR (U eiet) (3,5) 


gesetzt worden ist. 
Eliminiert man aus den Gin. (3,3) und (3,4) die Spannungen, so ergibt sich 
für den zeitfreien Verschiebungsvektor ® die Beziehung a ae 


wo B+ (A+ 2) graddiv 8 — uw rotrot = 0. , 


Zur weiteren Vereinfachung führt man außer der Divergenz den Rotor des 
Verschiebungsvektors ein: 


div 8 = U rot 8 = BW. (3,7) 


Unter Verwendung bekannter Beziehungen der Vektoranalysis erhält man 
danach durch Divergenz- bzw. Rotorbildung der Gl. (3,6) für U und ® die 
bekannten Wellengleichungen 


‘mit den Ausbreitungskonstanten Ver Ve 


w 9 
Vr A Vs u ° (3,10) 
19) Bei periodischen Vorgängen erweisen sich die Elastizitäts- und Viskositätskon- 
stanten im allgemeinen frequenzabhängig. Dies läßt sich auch theoretisch zeigen, falls 
man anstatt von den Gin. (2,1) auszugehen, lineare Beziehungen benutzt, die außerdem 
beliebig hohe zeitliche Ableitungen der Spannungs- und Verschiebungskomponenten ent- 
halten [vgl. T. Alfrey, jr. Mechanical Behaviour of High Polymers (NewYork 1948) 
S. 537]. Diese besonders zum Studium von Relaxationserscheinungen verwendeten 
Gleichungen führen beim Übergang auf komplexe Amplituden zu den gleichen Beziehungen 
(3,3), nur mit frequenzabhängigen Konstanten. Da im Rahmen dieser Arbeit die Fre 
quenzabhängigkeit eine untergeordnete Rolle spielt, wurde hier jedoch von den anschau- 
lichen Grundgleichungen (2,1) ausgegangen. 
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Die erste Wellenart entsprechend Gl. (3,8) (gekennzeichnet durch den Index k) 
nennt man Kompressionswelle, da hierbei nur eine reine Dehnung des Mediums 
eintritt. Die zweite Wellenart (3,9) (gekennzeichnet durch den Index s) 
heißt Scherwelle. Hierbei erfolgt nur eine reine Scherbeanspruchung. 

Die Ausbreitungskonstanten y,, y, sind im allgemeinen komplexe Größen, 
da A und u ebenfalls komplex sind. Setzt man Gl. (3,2) in (3,10) ein, so ist 
eine Trennung in Real- und Imaginärteil möglich: 


y=%,+tiß, j= (3,11) 

Nach H. T. O’Neil?!) wird « als die charakteristische Phasenkonstante 
und ß als die charakteristische Dämpfungskonstante bezeichnet, beides sind 
nach Definition Materialkonstanten. Ihre Bedeutung zeigt sich bei Betrach- 
tung der speziellen Lösung der Gin. (3,8) und (3,9) für eine homogene ebene 
Welle. Hierbei findet man, daß bei diesem Spezialfall x die Phasenlage der 
Welle bestimmt, während für das Abklingen der Amplitude 6 maßgebend ist. 
Im allgemeinen Fall ist der Zusammenhang aber komplizierter. 

Ganz allgemein gilt jedoch, daß eine Absorption der Welle nur auftritt, 
falls 8 von null verschieden ist. Eine notwendige und hinreichende Bedingung 
hierfür ist wiederum die Tatsache, daß die entsprechenden Viskositätskon- 
stanten von null verschieden sind. Beachtet man außerdem die Bedingung, 
daß sich Scherwellen nur in Stoffen ausbreiten können, bei denen u + 0 ist, 
was sich sofort aus Gl. (3,6) ablesen läßt, so kann man hierauf eine Einteilung 
der Medien vom akustischen Standpunkt aus aufbauen, die in Tabelle 2 dar- 
gestellt ist. Hierbei soll der Ausdruck Ausbreitung „gedämpfter Wellen‘ in 
der Tabelle in dem Sinne verstanden werden, daß in den entsprechenden 


Tabelle 2 
Einteilung der Medien vom akustischen Standpunkt 


| A, Me 


= 


und Scherwellen 


xXXXXXXX 


Gedämpfte Kompressions- und ungedämpfte 
Scherwellen 


Ungedämpfte Kompressions- und Scherwellen 


Gedämpfte Kompressionswellen 


| 
| Ungedämpfte Kompressionswellen ee 


20) An die dieser Stelle sei bemerkt, daß bei der Verwendung von Indizes in dieser Arbeit 
folgendermaßen verfahren wurde. Der Index j kann die Werte k (Kompressionswelle) 
oder 8 (Scherwelle) annehmen, während die Indizes m, n die Werte x, y, z durchlaufen 


können. ; Die — noch eingeführten Indizes p, q dagegen haben als Wertebereich die 
Zahlen 1,2... 


2) H. T. O'Neil, Physic, Rev. (2) 75, 928 (1949). 
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Stoffen „Schallabsorption‘ auftritt. Die Ergebnisse des nächsten Kapitels 
werden zeigen, daß zwischen beiden Größen noch ein kleiner Unterschied 
besteht. 

Die Wellengleichungen (3,8) und (3,9) stehen im Mittelpunkt der weiteren 
Betrachtungen. Sind sie gelöst, so erhält man den Verschiebungsvektor 
durch einfache Differentiation. Aus Gl. (3,6) gewinnt man 


B = U + rot B (3,12 
k > 


1 aU 1/aW, ow, iy 


Werden die Gln. (3,13) in die Beziehungen (3,3) eingesetzt, so folgen densi 
die entsprechenden Formeln fiir die Spannungskomponenten 


2 dy? yi \éxdy dy a 


Oy Gz 
4 2W, , ®W, 
yz Ox dz Oydz Ox dy 
Auf die Gin. (3,13) und (3,14) wird ebenfalls bei Betrachtung der Wellenaus- 
breitung im geschichteten Medium Bezug genommen werden. 


(3,14) 


ot, 


4. Lösung für ebene Wellen 
Die Gln. (3,8) und (3,9) sind beide von der Form ~~ 


4Q+y°Q=0, 
= Q eine skalare oder vektorielle Größe sein soll, was aber hier von ‘unter: 
geordneter Bedeutung ist, da in diesem Kapitel das Gemeinsame beider 
Wellenarten behandelt werden soll. Die Ausbreitungskonstante y ist, wie 
schon im letzten Kapitel gezeigt wurde, im allgemeinen komplex. 
Eine einfache Lösung der Gl. (4,1) ist die ebene Welle 


(4,2) 
in der neben dem Ortsvektor r ein Einheitsvektor n mit den Komponeten n, 


auftritt. Ist dieser Einheitsvektor reell, so breitet sich eine homogene ebene 
Welle aus, während bei komplexen n eine inhomogene ebene Welle auftritt. 
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Unter einer inhomogenen ebenen Welle”) versteht man bekanntlich eine _ 
Wellenform, die durch zwei Richtungen gekennzeichnet ist. Hierbei muß 
man unterscheiden zwischen dem Normalenvektor auf den Flächen kon- | 
stanter Phase, dem Phasenvektor n,, und dem Normalenvektor auf den _ 
Flächen konstanter Amplitude, dem Amplitudenvektor ng. Diese Tatsache — 
wird durch folgende Gleichung beschrieben: es 
Hierbei ist im Exponenten eine Trennung in Real- und Imaginärteil durchge- _ 
führt worden. Wie man aus Gl. (4,3) ablesen kann, wirkt in Richtung n, die 
Phasenkonstante «’ und in Richtung ng die Dämpfung ß’. In einer beliebigen © 
Richtung u, die mit n, den Winkel u, und mit ng den Winkel u, einschließt, er 
wirken a’ cos u, und ß’ cos us. Die Phase der Welle ändert sich also am _ 
stärksten in Richtung n, und die Amplitude in Richtung ny. Die reellen — 
Konstanten x’, ß’ werden zweckmäßig nach H. T. O’Neil?!) als die effektive 
Phasen- und Dämpfungskonstante bezeichnet. Zum Unterschied von den 
charakteristischen Phasen- und Dämpfungskonstanten, die Materialkon- — 
stanten sind, nehmen die effektiven Konstanten, wie im einzelnen noch gezeigt 
wird, mit dem Grade der Inhomogenität der Welle zu. Im folgenden soll kurz 
der Zusammenhang zwischen den Gl. (4,2) und (4,3) hergest lit werden. en 
Setzt man für den komplexen Vektor n Kar Re: 


n=a+ib, Ad (4,4) 


wobei a und b reelle Vektoren sind, so gewinnt man auf Grund der Tatsache, 
daß n ein Einheitsvektor ist, für a und b die Beziehungen = 
a-b=0 a? — b? = 1. (4,5) 
Die beiden Vektoren stehen also aufeinander senkrecht. Die Bedingung fir — 
ihre Beträge kann man erfüllen, falls man setzt = 
a = Coj pe, b = Sing ey (4,6) 
mit den Einheitsvektoren e,, e, in Richtung von a, 6 und der reellen Kon- © 
stanten p, die ein Maß für die Inhomogenität der Welle ist. ete 
Setzt man Gl. (4,4) unter Berücksichtigung von (4,6) in (4,2) ein, so liefert 
eine einfache Rechnung die effektiven Ausbreitungskonstanten und Aus- 
breitungsrichtungen als Funktion der charakteristischen Ausbreitungskon- 
stanten und von a und b: a 
1 


(x Coj ye, — B Sin ep) (4,9) 


"Var + (a2 + Sin? p 
V 


1 
— (B Coj pe, + « Sing e). 


~ + (a2 + BA) Gin? 


22) Die Notwendigkeit auch inhomogene Wellen in den Kreis unserer Betrachtung 
zu ziehen, zeigt sich schon bei einfachen Problemen der Brechung einer Welle an dr 
Grenzfläche zweier Medien, von denen das eine den Schall absorbiert, während imanderen 
keine Absorption auftreten soll. Eine einfallende homogene Welle erzeugt in diesem Fall 4 
an gebrochene inhomogene Welle, ganz gleich aus welcher Richtung der Schallstrahl _ 

ommt. 

2) Die verschiedenen Vorzeichen innerhalb der Klammer sind durch den negativen | 
Exponenten im Ansatz (3,1) bedingt. 2 
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Aus den Gin. (4,7) und (4,8) ersieht man, daß die effektive Phasen- und 
Dämpfungskonstante immer größer als die entsprechende charakteristische 
Konstante ist. Nach den Gin. (4,9) und (4,10) ergibt sich, daß die Vektoren 
1x, 1, in der Ebene liegen, die von e, und e, aufgespannt wird. Für den Winkels, 
den die beiden Richtungen n,, ng einschließen, erhält man bei skalarer Multi- 
plikation der letzten beiden Gleichungen die Beziehung 

a B 
V + (a? + Sin® [B® + (+) Sin? gy] 

Für die in Richtung n, wirkende Dämpfungskonstante ß’ cos e ergibt 
sich nach Gl. (4,11) ß’ cos e = x B/a’ <ß, sie ist also im Gegensatz zur effek- 
tiven Dämpfungskonstanten kleiner als die charakteristische Dämpfungs- 
konstante. 

Wie schon erwähnt ist @ ein Maß für die Inhomogenität der Welle. Hier 
seien die zwei Grenzfälle diskutiert. Setzt man m = 0, so wird nach Gl. (4,6) 
auch 6 = 0 und es ergibt sich «’ = «, p’ = B, n, = ng = a, € = 0. Es breitet 
sich dann eine homogene Welle aus. Geht andererseits 9 — 00, so laufen 
x’, ß’ auch gegen unendlich und im Grenzfall stehen die beiden Richtungen 
1x, 1g aufeinander senkrecht. Ist p endlich, so liegt der Winkel e zwischen 
0 und 2/2 und wächst monoton mit 9. Eine Ausnahme macht nur der Spezial- 
fall 6 = 0. Während für g = 0 sich hierbei eine ungedämpfte Welle aus- 
breitet, wird für jeden endlichen Wert von @ bereits e=n/2. An diesem 
Spezialfall läßt sich außerdem der kleine Unterschied zwischen Dämpfung 
und Absorption einer ebenen Schallwelle erkennen. Für endliches @ breitet 
sich nämlich eine gedämpfte Welle aus, es tritt aber trotzdem keine Absorp- 
tion auf. Während man von Ausbreitung gedämpfter Wellen spricht, sobald 
im Exponenten der Gin. (4,2) und (4,3) ein Realteil auftritt, gilt für den Fall 
der Schallabsorption noch die Zusatzbedingung, daß e + 2/2 sein muß, da 
andernfalls die in Richtung des Phasenvektors wirkende Komponente der 

Dämpfung null ist. 


cos € = 


(4,11) 


arg +2 
5. Die Schwingung des Verschiebungsvektors 
Die Betrachtungen des letzten Kapitels waren sowohl fiir ei 
wellen als auch für Scherwellen gültig. Bei der Untersuchung der Schwingung 
des Verschiebungsvektors bei den beiden Wellentypen muß jedoch eine 
Trennung erfolgen, da hierbei der physikalische Unterschied beider Wellen in 
Erscheinung tritt. 
I. Kompressionswellen 


Nach Gl. (3,12) gilt bei einer reinen Kompressionswelle für den Ver- 


schiebungsvektor ®, folgende 


. Vi 
7 . 
Unter Verwendung der Lösung einer ebenen Welle (4,2) RR 
U = 
mit der komplexen Konstanten U, 


U,=U,+il, 


Be 
| erhäl 
Um 
geht 
für 
mit 
‘A 
4 
Aus 
Seu 
tore 
man 
Br: 
zu 
Dan 
il 
Für 
baw 
b, = 
tun, 
der 
fall, 
Seh 
. 
Glie 
(9,1) 
; 
Sch 
(5,2) 
= - 


erhält man für ®, die einfache Formel TER OR Orr 


Um daraus die Schwingungsrichtung des Verschiebungsvektors zu bestimmen, 
geht man zweckmäßig auf die vollständige, reelle Lösung BURN: Dabei er- 


gibt sich nach (3,1) 


für den Verschiebungsvektor v, die Beziehung al 


= A, (P, cos 6, + sin (5,6) 
mit den Abkiirzungen ben 


= Oy (Uy x, + Ug B,) + br (U, By — Ue xy) (5,8) 


(5,5) 


= Oy — By r—wt. (5,10) 
r™ GI. (5,6) ersieht man, daß der Verschiebungsvek tor eine elliptische 
Schwingung durchführt. Die Ellipse liegt in der Fläche, die durch die Vek- 
toren a,, b, aufgespannt wird. Zur Bestimmung der Hauptachsen bildet 
man »?, wobei sich nach einfacher Umformung ergibt: 


vt = (U2 + U) (a + BD +, (5, (6,11) 


2 (U, + U, Bx) (U, By — Us xx) 
(Uy % + Us By)? (UB, Oty)? 


Damit kann man die Hauptachsen »,,, 0; der Schwingungsellipse angeben 
Dey = A, (Pr COS + Sin Y,) (5,13) 


= Ap (— Pe +4, (5,14) 
Für gewisse Spezialfälle entartet diese elliptische Schwingung in eine lineare 
bzw. in eine Kreisschwingung. Setzt man gy, = 0, so wird nach (4,6) auch 
6, = 0. In diesem Fall haben nach Gin. (5,7) und (5,8) p, und q, die Rich- 
tung von a,, es tritt also eine linear polarisierte Schwingung auf. Dabei liegt 
der Verschiebungsvektor in Richtung des Phasenvektors™). Für den Grenz- 
fall, daß gp, + oo geht, ergibt sich andererseits eine zirkular polarisierte 
Schwingung, da dann der periodische Anteil von v} gegenüber dem konstanten 
Glied nach (5,11) vernachlässigt werden kann. ee 
Bei der Ausbreitung einer reinen Scherwelle gilt zur Bestimmung der 
Schwingungsrichtung des Verschiebungsvektors nach Gl. (3,12) die Beziehung 
dB, == rot B. (5,15) 
fe 
2 og homogene Kompressionswelle wird deshalb auch als Longitudinalwelle be- 
ne 
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Die Lösung für W lautet nach (4,2 
Setzt man die Lésung (5,16) in 5 15) ein, erhält man ar, Flesiöt 
8,=— * (n, (5,17) 
Vs 
Hierbei ist zu bemerken, daß der komplexe Vektor WW, 


B, = B, + iB, 
nicht beliebig wählbar ist. Unter Beachtung der Tatsache, daß W der Rotor 
von ® ist und damit div® = 0 sein muß, ergibt sich für W, die einschrän- 
kende Bedingung 

n,=). (5,19) 
Nach einer analogen Rechnung wie bei der Ausbreitung von Kompressions- 
wellen erhält man für den reellen Verschiebungsvektor v, bei einer Scher- 
welle wieder eine elliptische Schwingung 


v, = A, (p, cos 6, + q, sin ö,) re) (5,20) 
3 
mit den Abkürzungen 
p, = o, (W, x b, + W,xa) — B, (Wx a, — Wy x b,) (5,21) 
= B, (BX b, + Wax a,) +0, (WX a, — BW, x b,) (5,22) 
6, = («, a, — 8, 6,) — (5,24) 


Zur Bestimmung der Hauptachsen der Ellipse bildet man wieder v?. Dabei 
erhält man ganz analog zu Gl. (5,11) 


|= + Bs) x a, — We Xx b,)? + (Wy b, + Wy x a,)? 


(5,25) 
+ 4 (B, W.)? cos 2 (6, — 
mit 
o, B, (Bi — Ws) — CH B) 
und die entsprechenden Gleichungen fiir die Hauptachsen 
ame A, (Ps COS Ys + qs sin Ys) (5,27) 
Veg = A, (— P, Sin y, + 9, COS Y,). (5,28) 


Diese elliptische Schwingung entartet für gewisse Spezialfälle wieder in eine 
lineare oder Kreisschwingung. Eine linear polarisierte Schwingung erhält 
man, wie man auch den Gin. (5,21), (5,22) ablesen kann, für den Fall, daß 
W, und ®, in der Ebene liegen, die von a, und b, aufgespannt wird. Die 
Vektoren p, und q, haben dann beide die Richtung a, x b,. 

Von Interesse sind außerdem die Grenzfälle 9, = 0 und 9, — oo. Für 
9, = 0 wird nach (4,6) auch b, = 0 und nach (5,21), (5,22) steht dann der 
Phasenvektor n,, senkrecht auf der Schwingungsellipse®). Für 9, > 0 


25) Eine homogene Scherwelle wird deshalb auch als Transversalwelle bezeichnet. 
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a Abb. 1. Das visko-elastische Schichtsystem 


26) Einschließlich der Halbräume I und II. 
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geht Formel (5,24) gegen ken 
A; 2 2 

= + (Wi cos* u, + Wy uy) | 

(5,29) 


+] ‘i — - + 4( ®,)? cos 2 (6, — y,)}; 


wobei mit w,, %, die Winkel zwischen ®, bzw. YW, und der Richtung a, xb, 
bezeichnet sind. Für uw, + 2/2 und nur solche Fälle sind hier von Interesse, - 
da der Fall v, = 2/2, wie oben gezeigt, auf eine linear polarisierte Schwingung . a 
führt, überwiegt bei v? das konstante Glied gegenüber dem periodischen — 
Anteil. Für 9, > co erhält man also im Grenzfall eine zirkular polarisierte 
Schwingung. 
Nach den bisherigen allgemeinen Betrachtungen über ebene Wellen n 
visko-elastischen Medien wird im folgenden die Wellenausbreitung in einem = 
System behandelt, daß aus n parallelen, ebenen Schichten besteht (Abb. 1). 
Jede Schicht 7 ist durch die Dicke d, und die komplexen Seniieibestianten, 
A, 4; gekennzeichnet. Die Oberflächennormale liegt in x-Richtung, in y-— 
und z-Richtung sind die Schiehten unbegrenzt. Vor der ersten und hinter der & 
n-ten Schicht befinden sich die Halbräume I und II. fr, 


6. Die akustischen Verhältnisse in geschichteten Medien 


on groBer Bedeutung fiir die sich ausbildenden akustischen Wellen- Rae 
alias in diesem Schichtsystem ist die Frage, aus welchen Stoffen die 
einzelnen Schichten bestehen. Unter Beachtung der Tabellen 1 und 2 lassen 
sich hierbei folgende drei Fälle unterscheiden: 

I. Sämtliche Schichten des Systems?) gehören der Gruppe A an. Au | 
der Trennfläche benachbarter Schichten ist dann die Stetigkeitsbedingung vn 
für sechs Größen zu erfüllen, und zwar für die drei Versobiebungskompo- ees 
nenten, sowie fiir die Normalspannung und die Schubspannungen. Dr 

II. Sämtliche Schichten des Systems2) gehören der Gruppe B an. 
der Trennfläche benachbarter Schichten ist die Stetigkeitsbedingung für 
zwei Größen zu erfüllen, und zwar für die Normalkomponente der Verschie- pe 
bung und die Normalspannung. Schubspannungen treten keine auf. ee 

III. Die Schichten des Systems?®) gehören teils der Gruppe A teils der 
Gruppe Ban. Als Stetigkeitsbedingungen gelten, falls benachbarte Schichten 
der gleichen Gruppe angehören, die entsprechenden Bedingungen dieser 
a Stoßen Medien aus verschiedenen Gruppen aufeinander, so müssen 
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die Bedingungen der Gruppe B erfüllt sein und außerdem auf der Grenzfläche 
die Schubspannungen verschwinden. 

Am einfachsten läßt sich der Fall II behandeln, da hier nur Kompressions- 
wellen auftreten können. Da dieser Fall aber wenig praktische Bedeutung hat, 
wurde auf eine Darstellung verzichtet. 

Komplizierter ist der Fall I, da sich hier Kompressions- und Scherwellen 
ausbreiten. Die Lösung für diesen Fall anzugeben, ist der Inhalt der folgenden 
Ausführungen, während der Fall III, der die größten Schwierigkeiten bereitet, 
in einer anschließenden Arbeit kurz behandelt werden soll. 

Im folgenden werde also angenommen, daß sämtliche Medien des Schicht- 
systems der Gruppe A angehören. Von Interesse sind nun die akustischen 
Wellenverhältnisse bei Einfall einer ebenen Kompressions- oder Scherwelle 
unter beliebigen Winkel aus dem linken Halbraum I. Bekanntlich erzeugt 
eine Kompressions- oder Scherwelle im allgemeinen bei Brechung und Re- 
flexion beide Wellenarten. Es entstehen also in den n Schichten je eine hin- 
und rücklaufende, außerdem noch im Medium I je eine reflektierte und im 
Medium II je eine auslaufende Kompressions- und Scherwelle. Im Mittel- 
punkt dieser Arbeit steht nun die Aufgabe, die im Halbraum I reflektierten 
und die im Halbraum II austretenden Wellen als Funktion der einfallenden 
Welle zu bestimmen. Die dazu notwendige Rechnung wird so geführt, daß 
man in jeder Schicht eine hin- und rücklaufende Kompressions- und Scher- 
welle ansetzt. Aus den Stetigkeitsbedingungen für den Verschiebungsvektor 
und die Normal- und Sehubspannungen an jeder Trennfläche erhält man 
gerade so viele Gleichungen als notwendig sind, um die Wellenamplituden 
in den n Schichten und dem rechten Halbraum sowie die reflektierten Wellen 
im linken Halbraum als Funktion der einfallenden Welle zu berechnen. 

Diese Rechnung läßt sich nun in übersichtlicher Weise gestalten, falls 
man eine Matrizenmethode anwendet, die aus der Elektrotechnik übernommen 
wurde. Mit Hilfe der Vierpoltheorie?’) behandelt man dort Probleme der 
Stromverteilung in Leiterketten, bei denen an einer Trennfläche Stetigkeit 
für zwei physikalische Größen (Strom und Spannung) verlangt wird, unter 
Benutzung von zweireihigen Matrizen. Auf Grund analoger Grundgleichungen 
und Randbedingungen läßt sich diese Methode unmittelbar auch in anderen 
Gebieten, wie z.B. Akustik, Optik und Wellenmechanik zur Lösung von 
Wellenproblemen in geschichteten Medien für den Spezialfall des senk- 
rechten Einfalls*) verwenden. Für den allgemeinen Fall der Ausbreitung 
von akustischen Wellen in geschichteten Medien bei beliebigem Einfalls- 
winkel muß die Methode erweitert werden. Da hier an einer Trennfläche die 
Stetigkeitsforderungen für sechs Komponenten (drei Verschiebungs- und drei 
Spannungskomponenten) zu erfüllen sind, erhält man hier im allgemeinen 
sechsreihige Matrizen). 

Im einzelnen wird die Rechnung so geführt, daß man zunächst nur eine 
beliebige Schicht betrachtet, in der man eine hin- und rücklaufende Kom- 

27) R. Feldtkeller, Einführung in die Vierpoltheorie der elektrischen Nachrichten- 
technik (Leipzig 1944) 4. Auflage. 

23) Für den Spezialfall eines ebenen Problems zerfallen die sechsreihigen Matrizen, 
wie noch im einzelnen gezeigt wird, in je eine vier- und eine zweireihige. Die Wellenaus- 
breitung in einem System rein-elastischer Medien bei Einfall einer homogenen Welle, 


die sich immer als ebenes Problem betrachten läßt, wurde von W. T. Thomson") und 
Y. Torikai!5) mit vierreihigen Matrizen behandelt. 
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pressions- und Scherwelle ansetzt. Für diese Schicht wird die Wellenmatrix ® 
aufgestellt. Diese liefert die Verschiebungs- und Spannungskomponenten 
in einer Ebene x = const als lineare Funktionen der Wellenamplituden der 
Kompressions- und Scherwellen. Anschließend wird speziell für x = 0 
die reziproke Matrix ®-! (0) gebildet. Diese beiden Matrizen werden später 
zur Berechnung des Reflexionsvermögens und der Durchlässigkeit des Schicht- 
systems in den Halbräumen I und II benötigt. Für die eigentlichen Schichten 
1 bis n wird dann schließlich die Kettenmatrix m?®) aufgestellt, die durch das 
Produkt B (d) B-1(0) entsteht. Diese bestimmt die Verschiebungs- und 
Spannungskomponenten der rechten Grenzfläche einer Schicht als Funktion 
der entsprechenden Größen an der linken Grenzfläche. Sind diese drei Ma- 
trizen ®, B-1(0) und m bekannt, so kann man in einfacher Weise das Re- 
flexionsvermögen und die Durchlässigkeit des Schichtsystems angeben. 


7. Die Wellenmatrix einer Schicht 


In diesem Kapitel wird die Wellenmatrix einer Schicht aufgestellt. Sie 
liefert, wie schon erwähnt, die Verschiebungs- 
und Spannungskomponenten in einer Ebene 


x=const als Funktion der Wellenampli- Aw 


tuden der Kompressions- und Scherwellen. 


Dazu betrachtet man eine beliebige visko- fay 


elastische Schicht i, die Grenzflächen seien 
durch die Koordinaten x, und z,festgeegt 
Im allgemeinen Fall breitet sich in dieser 
Schicht unter den hier vorliegenden Vor- ’ 2 
aussetzungen, je eine hin- und riicklaufende Abb.2. Visko-elastische Schicht 
ebene Kompressions- und Scherwelle aus. 
(7,1) 
W= + BW”. (7,2) 
Hierbei wird mit einem Strich die in positiver x-Richtung laufende und mit 
zwei Strichen die reflektierte Welle bezeichnet. Im einzelnen gilt nach Gl. (4,2) 


U’ = Up ete U" = Ug (7,3) 
= Wo ef y+ WW" = We et (7,4) 


mit den Wellenamplituden Uj, Uj’, Wo, 

Die komplexen Vektorkomponenten n,,, im Exponenten können hierin bei _ 
vorgegebenen Materialkonstanten der einzelnen Schichten als bekannt voraus- — 
gesetzt werden. Sie sind mit dem Ausbreitungsvektor der einfallenden Welle © 
im Halbraum I durch das Brechungsgesetz verkoppelt, daß sich aus den schon | 
mehrmals erwähnten Stetigkeitsbedingungen ergibt. Danach muß bekanntlich — 
falls wie hier die Flächennormale in x-Richtung zeigt — an einer Trennfläche | 
benachbarter Medien die Größe y, n, ; links und rechts gleich sein und außerdem _ 
für Kompressions- und Scherwellen denselben Wert annehmen, und ent- 
sprechendes gilt für y;n,;. Da y,n,, und y,n,, innerhalb einer Schicht kon- 
stant sind, werden sie damit zu Konstanten für das gesamte Schichtsystem. 
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Die x-Komponenten n,, ergeben sich anschließend aus der Nebenbedingung, 
daß der Ausbreitungsvektor ein Einheitsvektor ist. 

Um die Verschiebungs- und Spannungskomponenten als Funktion der 
Wellenamplituden anzugeben, benutzen wir die Beziehungen (3,13) und 
(3,14), die zunächst die Verschiebungs- und Spanungskomponenten als Funk- 
tion von U und ® liefern. Dabei sind von den Spannungskomponenten nur 
die von Interesse, die an einer Trennfläche x = const stetig sein müssen, 


das sind T,,, T,, und T,,. Die erforderlichen Gleichungen lauten: 


002 = \Oxdz Oudy 
(- Yi Oxdy dydz oy? ) 


Zur Berechnung der Verschiebungs- und Spannungskomponenten nach 
Gl. (7,5) müssen Ableitungen nach x, y und z gebildet werden. Führt man 
die Differentiation durch, so ergeben sich die Verschiebungs- und Spannungs- 
komponenten als lineare Funktionen der Wellenfunktionen®) 
+ + Ay + We’ + Us Wy + Wy’ 
+ WE + U WY 
HU + + Wag We + 
+ Wy, We + Au We 
Vy = U" + + We + We + Wy + Wi 
+ Ws; Wz + We 
Tey = Uy + Ug We + Ay Wi + Wh + Ay, WY 
+ WE + | 
V, = My + U WE HF 
+ Us, We + We 


+ Us, We + We 
Die Elemente &,, sind Funktionen der Ausbreitungskonstanten y,, sowie 
der komplexen Vektorkomponenten n,, ;, sie sollen im einzelnen nicht aufge- 
führt werden. Zu einer übersichtlichen Darstellung der Gleichungssysteme 
(7,6) gelangt man durch Übergang zur Matrizenschreibweise. Unter Ein- 


) Der Vorteil der in dem Gleichungssystem gewählten Reihenfolge der Zeilen und 
der Elemente innerhalb einer Zeile zeigt sich erst bei der Betrachtung des Spezialfalls 
ebener Probleme. 
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Eliminiert man in der Gl. (7,8) Wi, W/, so geht die rechteckige Matrix % in u 
eine quadratische Matrix W’ über Se 
Führt man an Stelle p die Matrix p’ 
u" 
Wy 
ein, so ergibt sich statt Gl. (7,8) die Matrixbeziehung 
O = p’. (2, 16) 


führung der Matrizen Q, und p 


U’ 

V, 
T W: 
| - - Ww" 
> Mee we 
Te we 
x 


geht das Gleichungssystem (7,6) in die folgende Beziehung über 

Q == A p. 
Da die Vektoren ®’ und ®’ nicht beliebig wählbar sind, kann man von den _ = 
Wellenfunktionen W;, W/, W; und WY, W/, WY je eine eliminieren. Nach 
Gl. (5,19) gilt bei einer fortschreitenden Scherwelle für die Wellenfunktion ®@’ 
die einschränkende Bedingung 


B’'-n, = 0 

oder in Komponenten 
Wi Nes + Wy ny, + Win, = 9. (7,10) 

Da diese Beziehung für jede Teilwelle erfüllt sein muß, erhält man für die > 
reflektierte Welle eine analoge Gleichung £ 


- + Wi n,, + Wi n,, = 0. (7,11) 
Aufgelöst nach Wi, und W2’ ergibt sich aus den Gin. (7,10) und (7,11) 


Um schlieBlich die Verschiebungs- und Spannungskomponenten als Funktion — 
der Wellenamplituden U5, Ug, Wio, Wio, anzugeben, benutzt man 
die Gln. (7,1) bis (7,4). Hierbei kann man unter Beachtung des Brechungs- 
gesetzes einen Schicht gleichen Faktor C 
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abspalten. Nach dem Einführen der Matrix der Wellenamplituden ® 
Us 
U, 


” 


z0 
Wyo 
Wee 
gewinnt man schließlich die gewünschte Beziehung, welche die Verschiebungs- 
und Spannungskomponenten Q als Funktion der Wellenamplituden ® liefert 


SPC]. (7,19) 


Für die Matrixelemente ®,, der Wellenmatrix ® ergeben sich unter Benutzung 
der Abkürzungen 


folgende Werte: 


(7,20) 


. 
—in R th... 
Bn e Bie = e 
Vr Vs 
Bis —# oris Bis 28 


e 
Vs Ys 
By = (A+ 2um Be = = —2yn,, n,, e* 
Ba, = Ny, Ba = 2 =—2un,,n,, 
= 2 ung, my, e* By = — 2 By nis) eis 
ul-—-2n)ei =2un,, Bas 2u ny, 
—i i INN; 
Vr Bsa - Vs Nzs 
Ys Mrs Vs zs 
= 2un,, n,, = —2 my, 
Ns ä Ba — 1)e* 
enes Problem: 
Von besonderem Interesse ist der Spezialfall, daß bei der im Malbeonsn: I 
einfallenden Welle der Amplitudenvektor n, und der Phasenvektor n,; mit 
dem Normalenvektor der Trennfläche np in einer Ebene liegen*!). Diese 
Ebene bleibt dann im gesamten Schichtsystem erhalten. Wählt man das 


31) Diese Bedingung ist z. B. bei Einfall einer homogenen Welle immer erfüllt, da 
dann n, und nt, die gleiche Richtung haben. 
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Koordinatensystem geeignet, so erhält man in diesem Fall nur eine a Ss 
keit von x und y. Da n,, = 0 wird, verschwinden sämtliche Matrixelemente, i 
die n,, als Faktor enthalten. 
= lautet die Matrix ® 
Brn Bie Bis Bu 
Ile Bs, Ban Bas Bau 
Bar Biz Bas Baa 0 
0 0 0 0 Bey Be, 
Die sechsreihige Matrix zerfällt also in eine vier- und eine zweireihige. Statt 
Gl. (7,19) erhält man in diesem Fall zwei Matrixbeziehungen. Die erste lautet _ 


= B, B,C 


mit den Abkürzungen 


Byn 


Bir 


Für die Matrixelemente 8 


findet man 
Vr 


upa 

Bn = - Byi2 = 

Bd,a=—-2un, n,,e* By og = 2M Nes Mys 


— tiny, 


Bya = Mg, By ag = — 2M 
Die zweite Beziehung ergibt sich ganz analog zu (7,28) ss 
B, C 


mit den Abkiirzungen 


0. = 


8. Die reziproke Wellenmatrix einer Schicht 


a Durch die reziproke Wellenmatrix ®-1 einer Schicht werden die Wellen- 
amplituden der Kompressions- und Scherwellen ® als Funläon der Ver- 
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schiebungs- und Spannungskomponenten © in einer Ebene x = const gegeben 


1 
= 8107|. | (8,1) 


Sie wird im Rahmen dieser Arbeit nur speziell fiir x = 0 benötigt. Die Ele- 


mente 8," (0) dieser Matrix werden nach bekannten siege der Matrizen- 
rechnung**) durch folgende Gleichung bestimmt: 


wobei mit U {®,,} die Unterdeterminante von ®,, und mit |®] die Deter- 
minante von ® bezeichnet wird. Hier seien nur die Ergebnisse dieser ein- 
fachen Rechnung mitgeteilt. Hierbei erhält man für die Elemente %;} (0) 
der reziproken Wellenmatrix folgende Werte®): 


iy, (1 -2 nz) R71 Vie 91 Nys | Vr Mys 

iv: 2 2 

Bis Vs 2uy, Nek a = 2 uy: 


-1__ Mys vs 


Vs 2 Mek 
— Nys 


-1 
1 —ty, Bso 


8 


= 


= == „m. — nz.) 
-1 . -1 
Bsi= VY, B52 = 


B55 2 (n 


en ein ebenes Problem zerfällt auch die reziproke Wellenmatrix im eine 

vier- und eine zweireihige. Statt Gl. (8,1) treten wieder zwei Beziehungen auf. 
Die erste lautet ». 


a) F. Neiss, Determinanten und Matrizen; Springer Verlag Berlin-Göttingen-Heidel- 
berg 1948. ae 
%#) Bei den Matrixelementen wurde das Argument 0 weggelassen. vr 1 
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mit den Elementen ®) 


_ ty s) i 


tat! 


(8,5) 


2 


— 


Für die zweite Begehung ergibt sich otk 


B= 820,05 (8,6) 


mit den Elementen %) 


—ty,n —1 —1 
$i, = Bu = Tu S745 = ws = 
A BEER 9. Die Kettenmatrix einer Schicht 


Die Kettenmatrix einer Schicht liefert die Verschiebungs- und en, 
komponenten der rechten Grenzfläche der Schicht als Funktion der ent- 
sprechenden Größen der linken Grenzfläche. Um sie aufzustellen betrachtet 
man wieder eine beliebige visko-elastische Schicht der Dicke d, deren Genz- 
flächen durch die Koordinaten x, und x, gegeben sind (Abb. 2). Für die Ab- 
leitung der Kettenmatrix werden nur die Wellenmatrix der Schicht und die 
entsprechende reziproke Matrix benötigt, die in den letzten beiden Kapiteln 
aufgestellt wurden. Im einzelnen geht man hierbei folgendermaßen vor: 
Zunächst kann man die Verschiebungs- und Spannungskomponenten an den 
Grenzflächen der Schicht mit Hilfe der Wellenmatrix als Funktion der Wellen- 
amplituden angeben. Für die linke Grenzfläche ergibt sich hierbei nach 
Gl. (7,19) 

Q (1) = B(x) BC, (9,1) 


wobei die Matrix der Verschiebungs- und Spannungskomponenten am linken 
Rand mit © (I) bezeichnet wird. Ganz entsprechend erhält man am rechten 
Rand © (II) 

(IL) = B (x) BC. (9,2) 


Unter Einführung der Plattendicke d kann man hierbei für die letzte Beziehung 
auch setzen 


© (ID) = B (x, + d) BC. (9,3) 


Um nun die Verschiebungs- und Spannungskomponenten am rechten Rand 
als Funktion der entsprechenden Größen am linken Rand, d.h. also Q (II) 
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als Funktion von © (I), zu erhalten, muß Beziehung (9,1) nach x aufgelöst 
werden. Unter Verwendung der reziproken Matrix lautet diese Umkehrung 


B= B(x) (9, 


Kombination der letzten Gleichung mit (9,2) liefert schließlich die gewünschte 
Beziehung 


mit der Kettenmatrix er 
m = B (x, + d) B- (x). 

Hierbei läßt sich das Produkt (9,6) noch etwas vereinfachen. Wie man ohne 
weiteres einsieht, ist die Kettenmatrix einer Schicht unabhangig von der Lage 
des Koordinatensystems. Man erhält also bei Gl. (9,6) für beliebiges x, immer 
den gleichen Wert. Da kann man aber ohne weiteres 2, auch gleich null 
setzen, so daB sich schlieBlich die Kettenmatrix einer Schicht aus der ein- 
facheren Beziehung 


m = B(d) 8! (0) (9,7) 


berechnen läßt. Führt man die Multiplikation (9,7) durch, so ergeben sich 


09 
für die Matrixelemente m,, der Kettenmatrix folgende Werte: are j 


Ya Mor - . 
2 Ny, n,, ny, (1 — 2 — 2n},) 


dent i E Nek Nys Mes sin x + — Nz, (4 — Giz 2 Mes) sing 
v 


n 


“(cos # — cos 0) 
Ys 


ye . 
Mes) sin x — MY, (1 2 ) sina 


Fer ola Mee 


= — 2iuy,n,,(1— 2n2,) (cos x — cos a) 
(1 — 2n2,) . 

sin x + 2n,,n,,sin o] 

Nek hank 


— 2iuy,n n,,(1—2n},) (cos x — cos 0) 
2n Er 
= —i 


sinx + 2n,, m,, sin o| 
ak 
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= My 
= 2 n? ) cosa 
ys n?, 


Nyx Bye 
BY,” 

=2n,, 

Mek 

Ms 

MY, (nis + {1 — 4n2,]) 

Nys 


sino 


= — Max My, — 
= 


= — AMY, Mys Sin — 


Nys Nes (1 — 4n3,) sin 


35 
2 n?, cos x + (1— 
Nee Nes a= 
= 
LYs "es + LY, % 
2. 
wy, (m3, + 


Vp Ny, „sin sin o 


= Wigs. 
2: 


Beim ebenen Problem zerfallen, wie in den letzten Kapiteln gezeigt wurde, 
die Wellenmatrix und die reziproke Wellenmatrix im gleichen Sinne. Nach 
einem bekannten Satz der Matrizenrechnung*®) tritt dann bei der Ketten- 
matrix, die aus dem Produkt beider gebildet wird, das gleiche auf. Statt 
Gl. (9,5) ergeben sich beim ebenen Problem zwei ere Die erste lautet 


u O, (I) = m, (I ) 


mit der Kettenmatrix 


m, = 8, (d) (0) 
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mn = — (1 — 2 n2,) cos x + 2 n?, cosa = 


Vr Mak 
= uy sinx + 
[2 Mer Ms my, (1—2 
= (cos x — Cos a) 
1-2n},).. 
M,sı =) sinx—4uy, N,,n,,sino 
1—2n2 
= — 2tpy, (1 — 2 n;,) (cos x — cos a) 
1-2 
inc + 2 sin o 


M32 = My 


Ma = 222, 


n N, sin . 


Ma = 


Mya. = 
My, 43 = 4uy, ‘ak n? , sin 
M, 44 = Myg3- 


Ys 
Nes 

ne Für die zweite Beziehung erhält man ganz analog zu (9,11) 
(Il) = m, ©, (1) 


omit der Kettenmatrix 


m, = 8B, (d) 8," (0) 
den Matrixelementen 


M,5; = cos 0 = 


= 


= — Meg SN Migs = 


10. Die Kettenmatrix eines Schichtsystems 

In diesem Kapitel soll der Fall behandelt werden, daß n Schichten be- 
liebiger Dicke und beliebiger Materialkonstanten hintereinander liegen. 

; Dazu werde zunächst der Fall von zwei Schichten betrachtet. In einer 
beliebigen Schicht gilt nach den Ableitungen des letzten Kapitels die Ma- 

 trixbeaiehung 


(IIT) = mQ (I). (10,1) 


(IT) = m, Q, (1) (10,2) 
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und analog in der zweiten Schicht 

(11) = m, O, (II). (10,3) 
Da die beiden Schichten hintereinander liegen, gilt auf Grund der Stetigkeits- 
bedingung 


Qu (11) = OQ, (1) (10,4) 
und man kann Gl. (10,2) in (10,3) einsetzen. Dabei ergibt sich 
© (II) = m, [m, O, (D]. (10,5) 


Da bei einem Produkt von Matrizen das assoziative Gesetz gilt, kann man 


dafür auch schreiben 
Oy (I) = (m, Q, (I). (10,6) 


Damit hat man wieder die Werte der Verschiebungs- und Spannungskom- 
ponenten am rechten Rand der Doppelschicht als Funktion der entsprechenden 
Werte am linken Rand erhalten. Setzt man zur Abkürzung 


M = m, m,, whisk (10,9) 


so folgt bei einer Doppelschicht eine Beziehung, die ganz analog Gl. (10,1) ge- 
baut ist: 
© (II) = MOM. (10,10) 
Bei dem Matrizenprodukt (10,9) erhält man die Elemente der resultierenden 
Matrix M nach der bekannten Regel 
6 
Mya = My pi My, ia (10,11) 
Hat man n Schichten, so führen ganz entsprechende Überlegungen zu dem 
Ergebnis 


© (IT) = M Q (I) (10,12) 

mit den Abkiirzungen aad 
() = 0,18) 


M =m, m,_,... m,m |. (10,15) 


Die Kettenmatrix eines Schichtsystems gewinnt man also in 
einfacher Weise durch Multiplikation der Kettenmatrizen der 
einzelnen Schichten, wobei die Reihenfolge gerade umgekehrt 
zur geometrischen Anordnung ist. Für die Elemente der resultierenden 
Matrix findet man analog zu (10,11) 

6 6 


1 = 


Statt Gl. (10,16) anzuwenden, erscheint es jedoch a (rn — 1)-mal 
hintereinander die Beziehung (10,11) zu beniitzen. 
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Ebenes Problem: 

Liegt ein ebenes Problem vor, so gelten, wie man sofort einsieht, ganz ent- 
sprechende Regeln für die vierreihige Matrix I, und die zweireihige Matrix W, 
Statt Gl. (10,12) treten zwei Beziehungen auf. Die erste lautet 


©, (II) = M, Q, (1) A (10,17) 
analog ergibt sich die zweite Gleichung 
mit der Kettenmatrix q 


11. Das Reflexionsvermögen eines Schichtsystems 2 


Systems M sei als bekannt vorausgesetzt. Eine einfallende Kompressions- 
oder Scherwelle aus dem Halbraum I erzeugt an der Trennfläche x = —D 
 reflektierte Wellen, die im folgenden als Funktion der einfallenden Welle be- 
ae stimmt werden sollen. (Vgl. Abb. 1.) Dazu benötigt man außer der Ketten- 
matrix des Systems für den Halbraum I die Wellenmatrix und für den Halb- 
sie raum II die reziproke Wellenmatrix. 


fi Zunächst wird durch die Kettenmatrix ein Zusammenhang zwischen den 
vs ot Verschiebungs- und Spannungskomponenten am rechten und linken Rand des 
Systems hergestellt. Nach Gl. (10,12) gilt 
© (Il) = MQ (D. 1) 
se Außerdem kann man im Halbraum I für die Frontfläche des Sy stems x = —D 
die Wellenmatrix ansetzen: 
(— D) = Bi (— D) Bı C. (11,2) 
läßt sich im Halbraum II für die zr=0 die ‚reziproke 


Wellenmatrix aufstellen: 


Bu = Vil (0) (0) 


Nun gilt aber # ang 
— D) = 4 
Q (0) = O (ID). (11,5) 


Damit läßt sich durch Kombination der Beziehungen (11,1)—(11,3) $y als 
Funktion von ®; angeben: 


M, = m, Mn-ı-- - (10,20) 


Ein beliebiges Schichtsystem werde vorgegeben, und die Matrix des 
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mit der Matrix % 


3 = Bil (0) M Bz (— D) 


Der funktionelle Zusammenhang zwischen den Wellenampli- 
tudenim linken und rechten Halbraum wird also durch die Matrix § 
hergestellt,die entsteht, falls man die Kettenmatrix des Schicht- 
systems von rechts mit der Wellenmatrix an der Stelle x= —D 
bezogen auf das Medium I und von links mit der reziproken Wellen- 
matrix an der Stelle x = 0 bezogen auf das Medium II multipliziert. 
Es wird angenommen, daß diese Multiplikation durchgeführt sei, so daB § 
als bekannt vorausgesetzt werden kann. Die Gl. (11,6) steht im Mittelpunkt 
der weiteren Betrachtungen, und die im Rahmen dieser Arbeit aufgeworfenen 
Fragen werden durch die Diskussion der Beziehung (11,6) gelöst. In ausführ- 
licher Schreibweise lautet sie 


= Vor + + Bis + Bra Wi01 + Wyor + Bis 
0=Uon Bar Vor + + Bes + W201 + Bes Wyoı + Bas 
Wrou = Vo1 + Fae Vor + Wioı + W201 + Fas yor + Sas Wyoı 

0= Win = Vor + Vor Sas Wro1 + Sas Wrot Sas 
= Vor + Vor + Bsa Wioıt Sea + Bes Wor + 

Ugı + Sea Uor + W2ot + Bea Wror - Wyo1 See Wyot 
(11,8) 


Bei diesem Gleichungssystem ist die linke Seite der zweiten, vierten und 
sechsten Gleichung null, da im Halbraum II keine reflektierten Wellen auf- 
treten können. Für die weitere Rechnung ist es zweckmäßig, diese drei Glei- 
chungen getrennt zu betrachten. Dies führt auf zwei Gleichungssysteme. Das 
erste lautet 


= Vo + Uo + Bis Wio + Gra + Woo + Tre 7 
= Uo + Ue + Wee + Bae + Bas Woo + Fae Wro (11,9) 
yo = B51 Uo + Bsa Uo + Bsa Win + Bsa Veo + ss + Sse 


Für das zweite ergibt sich ER satel 


0 = Uo + Uo’ + Fas Boa + Ges Wyo + Ges Wyo 
9 = Fa + Sag UE + Bag + Wo + Bas Wio + Sas (11,10) 


0 = Fe Ue + Sea Uo + Seg Wo + ea Veo + Bes Myo + es Wyo- 


Hierbei sind bei den Wellenamplituden die Indizes I bzw. II weggelassen 
worden. Zur Unterscheidung wurde statt dessen eine auslaufende Welle im 
Halbraum II durch drei Striche gekennzeichnet. Zur Bestimmung des Re- ; 
flexionsvermögens des Schichtsystems dient Gleichungssystem (11,10). Nach 
Auflösung erhält man die Amplituden der reflektierten Wellen. Be 
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Bei einer einfallenden Kompressionswelle findet man 
Sar Gas Das 


der Abkürzung Wyo = - - U6 inal: 
| Sas Sas i, 
23 0402 25 O28 ‘th 


us | Ses Bea Wo Bas Bea See Wee (11, 15) 


ER 
vo 67 20 vo ’ 
Der gleichzeitige Einfall einer Kompressions- und einer Scherwelle läßt sich 


mit dem Gleichungssystem (11,8) nur behandeln, falls zwischen den Aus 
breitungsrichtungen der beiden Wellentypen das Brechungsgesetz erfüllt ist. 
Im allgemeinen wird dies jedoch nicht der Fall sein, und man gewinnt dann 


Ebenes Problem: 


Da beim ebenen Problem die Wellenmatrix, reziproke Wellenmatrix und 
die Kettenmatrix im gleichen Sinn zerfallen, muß auch die Matrix %, die den 
Zusammenhang zwischen den Wellenamplituden im linken und rechten Halb- 
raum herstellt, dem gleichen Gesetz gehorchen. Damit erhält man an Stelle 
der Gl. (11,6) zwei Beziehungen. Die erste lautet = © 


Bu = Oy P,ı 4 u (11, 18) 
mit der Matrix %, tikes. seh: 


B= M, By (—D. 
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(11,12) 


11,13) 


11,14) 


11,15) 


11,16) 


M Pu P.ı 
mit der Matrix %, 

Da nach Gl. (7,27) ®, nur eine Funktion von W, ist, bedeutet dies, daß beim R 
ebenen Problem die Scherwellenkomponente W, unabhängig von W, und U 


durch das Schichtsystem geht, während die Scherwellenkomponente. W,mit 3 


der Kompressionswelle U in Wechselwirkung tritt. Hierbei ist zu berück- 


sichtigen, daß bei der hier gewählten Darstellung die ausgezeichnete Fläche a E 


in der xy-Ebene liegt. Wie man sofort einsieht, gilt unabhängig vom Koor- — 
dinatensystem ganz allgemein: 
Die Scherwellenkomponente W,, die senkrecht auf der ausgezeichneten — 


Fläche steht, tritt in Wechselwirkung mit der Kompressionswelle, während die _ ar = 


Scherwellenkomponente W,, die parallel zur ausgezeichneten Fläche liegt, 
sich unabhängig ausbreitet. 
Dieser vereinfachte physikalische Sachverhalt beim ebenen Problem 


liefert auch für die Bestimmung des Reflexionsvermögens eines Schicht- 


systems einfachere Beziehungen. Für diese Aufgabe werden die Gln. (11,18) _ 
und (11,20) wieder in ausführlicher Schreibweise aufgestellt. Hierbei gewinnt 


man zwei Gleichungssysteme 
| 
U = Vo + Uo + Wee t+ Bu 


o 
I 


+ See Us + Beg W;o + Wr 
Wo = Bar Uo + Bae Vo + Wee + Sea 


(11,22) 


und 
Wi = Tas Who > Tse Were 


0 = Bes Wio + See Wro- 


Für die des Reflexionsvermögens sind bei diesen Systemen 
wieder nur die Gleichungen von Interesse, bei denen links null steht. Nach 
einfacher Auflösung folgt das Ergebnis. Bei einer einfallenden Kompressions- 
welle findet man 


“4 7’ yu ’ 


(11,25) 


Eine einfallende Scherwelle liefert | 
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12. Die Durchlässigkeit eines Schichtsystems 


Nachdem die reflektierten Wellen angegeben wurden, können unter Ver- 
wendung der Ableitungen des letzten Kapitels auch die durchgehenden 
Wellen berechnet werden. Dazu müssen die Lösungen (11,11) bis (11,17) 
in das Gleichungssystem (11,9) eingesetzt werden. 
so ergibt sich bei einer einfallenden Kompressionswelle 


a1 | 
Das | — Vaz Var das 


| 


> > | 
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Eine einfallende Scherwelle führt auf das Ergebnis 
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Wird dies durchgeführt 
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Ebenes Problem: 
Ver- Beim ebenen Problem erhält man das Ergebnis auf analoge Weise. Es 
aden werden also die Lösungen für die reflektierten Wellen (11,24) bis (11,26) 
1 ‘) in die Gleichungssysteme (11,22) und (11,23) eingesetzt. Bei einer einfallenden 
ihrt,  Kompressionswelle findet man 
Sa Öa Oa By 
Bei einer einfallenden Scherwelle folgen schlieBlich die ar 


> Vee | 


Wiyo= [055 Wyo (12,11) 


Damit sind die im Rahmen dieser Arbeit aufgeworfenen Fragen gelöst. 
Zum Schluß sei nochmals darauf hingewiesen, daß bei der Wellenausbreitung 
im Schichtsystem vorausgesetzt wurde, daß sämtliche Medien der Gruppe A 
der Tabellen 1 und 2 angehören, die durch das Auftreten von Schubspanungen 
bzw. von Kompressions- und Scherwellen gekennzeichnet sind. Nur diese 
Annahme ermöglicht die geschlossene Darstellung der letzten beiden Kapitel. 
(12,4) Im allgemeinen Fall, daß die Medien teils der Gruppe A und teils der Gruppe B 
angehören, müssen weitere Spezialfälle unterschieden werden, was in einer 
anschließenden Arbeit kurz behandelt werden soll. 


Herrn Prof. Dr. A. Papapetrou möchte ich an dieser Stelle meinen 
besonderen Dank für das fördernde Interesse an der Arbeit aussprechen. 


Berlin-Adlershof, Laboratorium für Kunststoffe der ome 
12,5) Akademie der Wissenschaften zu Berlin. 


Bei der Redaktion eingegangen am 1. Dezember 1955. 
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A 
b Sound Absorption in Liquids: _ 


# 
By S. Parthasarathy, M. Pancholy and S.S. Mathur 


Abstraet 


The thermal method developed in this Laboratory for the measurement of 
the ultrasonie absorption coefficient in liquids has been used elsewhere in a 
modified form. Re-examination of this new technique reveals that it is not 
capable of giving a value for the ultrasonic absorption coefficient of liquids. 


3 : It was shown!) in this Laboratory that if an ultrasonic beam having an 

Ben! intensity / at the source is fully absorbed in a medium, it produces an amount of 
heat H which is related to the sound intensity by the relation / = JH where J 
is the Joules constant for the mechanical equivalent of heat and has the value 
of 4,19 - 10? erg/see. On the basis of this principle a method?) was developed 
to determine the ultrasonic absorption coefficient of liquids by the measure- 
ment of the ultrasonic intensity at different points along the path of the beam 
by absorbing the whole of the beam at that point to produce heat in a liquid, 
the rise in temperature of the liquid being a measure of the heat produced by 
ultrasonics. The absorption coefficient can be determined from the relation 
H = H, e~*** where / in the usual attenuation formula is replaced by H as a 
consequence of the equivalence relation. Using this technique, absorption 
coefficients for a number of liquids were determined. The values for benzene 
and carbon disulphide as obtained by this method (and confirmed? later by 
repeating the experiment) are about 800 - 10-1” em—" sec? and 7500 - 10-1?em-! 
sec? respectively. The results thus obtained are in close agreement with those 


1) S. Parthasarathy, S. S. Chari and D. Srinivasan, Ann. Physik 12, 8 (1953); 
J. Acoust. Soc. Amer. 25, 335 (1953). 


2) S. Parthasarathy, S. S. Chari and P. P. Mahendroo, Z. Naturforsch. 8a, 
273 (1953); J. Physique Radium 14, 366 (1953); S. Parthasarathy, P. P. Mahendroo. 
Il Nuovo Cimento 10, 1196 (1953); S. Parthasarathy, P. P. Mahendroo and S. $ 
Mathur, J. Sci. Industr. Res. 12 (B) 457 (1953). 


2 Se, S.Parthasarathy, M. Pancholy and C. B. Tipnis, Nature 176, 611 (1955). 
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obtained by other methods of measurement, such as by radiation pressure, 
pulse technique, the new thermal etc‘). 

Recently Dutta and Samal°) have also carried out some experiments, on 
the basis of the equivalence relation, to determine the absorption coefficients 
for liquids. According to them, the absorption coefficients of benzene and 
carbon disulphide appear to be 75 - 10-17 em! sec? and 115 - 10-17 em! sec? 
respectively. These values are of a different order of magnitude compared 
to the usually accepted values. Their conclusions have been carefully examined 
by us. The object of this paper is to re-examine the technique and report the 
results of our investigations. 


2. Experiment 


The experimental arrangement employed is the same as that of Dutta 
and Samal. 

A beam of ultrasonic waves is produced in kerosene surrounded by a con- 
stant temperature bath. The experimental liquid is contained in a double 
walled metallic cylindrical vessel having asbestos packing. Two thin mica 
windows are fitted on the faces of the cylindrical vessel to allow the ultrasonic 
beam to pass in and out after absorption in the liquid. The vessel is put in 
the path of a strong ultrasonic beam produced in a tank containing kerosene 
and the steady difference in temperatures between the liquid in the vessel and 
the kerosene bath is recorded. 

If J be the intensity of the ultrasonic beam after passing through the 
first mica window and x, the length of the vessel, the ultrasonic energy ab- 


sorbed by the liquid every second is 


This energy is determined by producing the same difference in temperatures 
by heating the liquid electrically. Exactly the same procedure is followed for 
another vessel which differs from the first only in length, care being taken to 
keep the first mica window of this vessel at the same position relative to the 
erystal as that of the first. If x, be the length of this vessel, the ultrasonic 
energy absorbed per second in this case will be 


I {1 — exp (— 2,)} 


and can be found as before. The ratio of the ultrasonic energies absorbed in 
the two vessels (which was used by Dutta and Samal to estimate the value 
of x) can thus be determined. 

The experiment was performed with an ultrasonic beam having a frequency 
of about 1 Me./sec. Three vessels of the type described above and having 
lengths equal to 4,55 em, 3,6 em and 2,3 em respectively were used. 


*) L. Bergmann, “Der Ultraschall’ (1954 edition); J. J. Markham, R. T. Beyer 
and R. B. Lindsay, Rev. mod. Physics 28, 353 (1951); P. Biquard, Rev. Acoust 8, 
130 (1939); S. Parthasarathy, S. S. Chari and D. Srinivasan, J. Physique Radium 
14, 541 (1953); S. Parthasarathy, C. B. Tipnis and M. Pancholy, Z. Physik 140, 
156 (1955); G. W. Willard, J. Acoust. Soc. Amer. 12, 438 (1941); J. M. M. Pinkerton, 
Proc. Phys. Soc. [London] 62 (B) 129 (1949); S. Parthasarathy, C. B. Tipnis and 
M. Pancholy, Z. Physik 140, 504 (1955); J. R. Pellam and J. K. Galt, J. chem. Physics 

5) A.K. Dutta and K. Samal, Nature 174, 976 (1954). er 
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3. Results 2 


In graph (1) we heve shown the steady difference in temperatures that 
could be produced by a given electrical power for each of the three vessels 


3.6 


STEADY DIFFERENCE IN TEMPERATURES- °C 


ns 


fit 


used. The curves (1), (2) and (3) 
represent the calibration for the small, 
medium and the large vessels, having 
lengths equal to 2,3 em, 3,6 cm and 
4,55 em respectively. 

The first three columns of Table! 


. show for the various liquids, the 


steady difference in temperatures be- 
tween the liquid in the vessel and 
the kerosene bath, in each vessel, as 
maintained by an ultrasonic beam of 
constant intensity at the source. The 
electrical power required to produce 
the same difference in temperatures 
in any vessel is read off the correspon- 
ding calibration curve and is recorded 
in the last three columns of this 
Table. This energy for the small, 
medium and the large vessels is 
represented by W,, W, and W, 
respectively. 


4. Discussion 


04 Table II represents the ration of 
‘ ved the ultrasonic energies absorbed in 
= » any two of these vessels. The ratio 
e N vessels is also given for comparison. 
5 2 Electrical power to produce 
the same diff. in temp. 
ae (from the graph) in 
+ 
Liquid Small | Medium | Large u vessel | un 
vessel vessel | 2,3cm | 3,6cm | 4,55 cm 
3 i ry: 2 ey em| (watt) | (watt) | (watt) 
| 1 2 3 
1. Benzene 20 | 20 2,0 0,225 | 0,340 | 0,430 
2. Chloroform 2,9 2,9 2,9 0,320 0,490 | 0,620 
3. Xylene 2,5 2,5 2,5 0,280 | 0,420 | 0,540 


A glance at Table II shows that for all the three liquids used namely 
benzene, chloroform and xylene, which are known to have widely different 
values for the absorption coefficient (being as 10:3: 1) the ratio of the energies 
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Table II 
Ratio of Ratio of Ratio of 


Ultra- | Ultra- Ultra- 
Liquid Ellee. | Elec. Elec. 
AS Powers | paths | Power | paths | Power | paths 
za | | elt, | | 

1. Benzene 191 1,98 1,51 1,57 1,26 1,26 
9 Chloroform | 1,94 1,98 1,53 | 1,57 1,27 1,26 
3. Xylene wool Ve 1,98 1,50 | 1,57 | 1,29 1,26 


absorbed in any two of the vessels is the same. Exactly the same result 
was obtained by Dutta and Samal5). According to them, the ratio of the 
energies absorbed in two vessels having lengths equal to 4,6cm and 3,6 cm 
respectively was the same, namely 1,25, for both benzene and xylene. From 
these results they concluded that the absorption coefficients of these liquids are 
approximately the same. If, however, we compare this ratio with the ratio 
of the ultrasonic paths in the corresponding two vessels, we find that the two 
agree very closely (Table II). It appears therefore that the ratio of the energies 
absorbed in any two vessels is not a characteristic of the liquid but is only a 
function of the lengths of the vessels used, being equal to the ratio of the 
ultrasonic paths traversed in them. This conclusion admits of a simple ex- 
planation. 


It is known that «/»? for most of the liquids is of the order of 10-17 em sec?. 
It will be quite reasonable to assume that at a frequency of the order of a 
few megacycles per second and the length x = 3—4 cm the product « x is 
much less than 1. In the expansion for e~*** it is therefore justifiable to 
retain terms only upto the first power of x. The ratio of the ultrasonic energies 
absorbed in the two vessels thus reduces to the ratio of the lengths of the 
vessels, namely 2,/2,. This ratio is independent of the nature of the experi- 
mental liquid. Since from the order of accuracy attainable in such experi- 
ments it is not possible to account for the second order terms in x it appears 
that the above procedure cannot be used for the measurement of « specially 
for liquids having lower absorption coefficients. If however we increase 
the product « x by using a higher frequency for the ultrasonic beam, the ab- 
sorption will increase as v2 and it is likely that the beam may not pass through 
the vessel i.e. it will be absorbed completely before traversing the whole 
length of the vessel. We therefore conclude that the above procedure is not 
suited to give a value of « and that any interpretation of the above results as 


giving a value for « will be fallacious. 


5. Conclusion 


When a beam of ultrasonic waves of high intensity passes through a liquid 
contained in vessels of different lengths necessarily smaller than the length 
of the ultrasonic beam, then the steady difference of temperatures is the same 
for that liquid, but the energy (thermal or electrical) required for this steady 
state is different for each vessel depending on the path. The ratios of energies 
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in two such vessels is obviously the ratio of the lengths of the two beams of 
ultrasonic waves trapped in the vessels, whatever be the absorption coefficient 
of the liquid contained in the vessel. This experimental method does not 
enable one to determine the ultrasonic absorption coefficient of liquids as has 
been supposed by Dutta and Samal. 


New Delhi, National Physical Laboratory of India. 
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Wissenschaftliche Selbstbiographie 


3. Auflage. 34 Seiten mit einem Bildnis, 
1955. DIN A 5. DM 1.35 


Bildet: Lebensdarstellungen deutscher 
Naturforscher, Nr. 5. Herausgegeben von 
der Deutschen Akademie der Natur- 
forscher (Leopoldina). 
Physikal. Blatter: Planck hat uns eine 
Schrift nachgelassen, die neben dem Ent- 
wicklungsgang eines der größten Gelehr- 
ten manche wertvolle Einzelheit enthalt. 


JOHANN AMBROSIUS BARTH /VERLAG / LEIPZIG 


Handbuch der ältesten 


vorgeschichtlichen Metallurgie 
in Mitteleuropa 
Von Dr. HELMUT OTTO, Leuna, und Dr. Dr. WILHELM WITTER ¢ 


V, 222 Seiten mit 33 Tafeln. 5 Karten und zahlreichen Abbildungen 
in 31 Analysentabellen. 1952. 4°. DM 63.40, Leinen DM 65.40 


Chemische Technik: Ein außergewöhnliches Buch liegt hier vor. Dadurch, daß der 
Chemiker den Vorgeschichtler unterstützt, ließen sich bindende Schlüsse über Herkunft und 
Herstellungsart metallischer Gegenstinde der frühen Bronzezeit Mitteleuropas ziehen, 
die von größter Bedeutung für die Kenntnis dieser Zeit sind. In Form von fast einein- 
halbtausend Analysen von Funden aus ganz Europa und darüber hinaus ist die Forschungs- 
arbeit rund zweier Jahrzehnte von Wilhelm Witter, der leider das Erscheinen des Wer- 
kes nicht mehr erleben durfte, und seinem Mitarbeiter Helmut Otto niedergelegt. Auf ganz 
einfachen Methoden beruhend, ist hier ein geniales Werk von internationaler Bedeutung 
entstanden. E. Pechmann, Berlin 


Arbeitstagung Festkörperphysik Il 


vom 8. bis 11. April 1954 in Dresden 
Redaktion Dr. K. W. Béer, Berlin 


219 Seiten mit 205 Abbildungen im Text. 1955. Gr. 8°. DM 21.— 


Neue Hütte: Der dargebotene Stoff hat keineswegs nur für den theoretischen Physiker 
oder Physiko-Chemiker, sondern auch für den physikalisch-chemisch orientierten Metall- 
kundler besonderes Interesse. Man kann sagen, daß die Niederlegung der Tagungsberichte 
in Buchform sehr wertvoll ist, da sich in den Vorträgen eine große Zahl fruchtbarer, Ge- 
danken und Ansätze zur Theorienbildung zeigt. Die Beiträge sind von hohem wissen- 
schaftlichem Niveau. Thiimmler 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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Ost-Rassow 


ps Lehrbuch der chemischen Technologie 


26., neubearbeitete Auflage 


Herausgegeben von Professor Dr. Berthold Rassow + und Professer Dr. Wilhelm 
Karl Schwarze, Leipzig. XVI, 1269 Seiten mit 541 Abbildungen und 203 Tabellen 
im Text und 13 Tafeln. 1955. Gr. 8°. Ganzleinen DM 32.80 


Glas-, Email-Keramo-Technik: Es ist darauf hinzuweisen, daß das vorliegende 
Werk seit dem Jahre 1899, dem Erscheinungsjahr der ersten Auflage, die Aufgabe er- 
füllt, den Studierenden in das Gesamtgebiet der chemischen Technologie einzuführen und 
dem Praktiker und allen denen, die sich über chemische Industrien und ihre Erzeugnisse 
zu unterrichten haben, als Lehr- und Nachschlagebuch zu dienen. Der Wert dieses als 
Standardwerk der technischen Chemie zu bezeichnenden Buches liegt darin, daß es als 
allgemeines technisches Nachschlagewerk unerreicht ist, was insbesondere auch durch die 


hohe Auflagezahl unter Beweis gestellt wird. Karmaus 


Höhere Mathematik für den Praktiker 
Von Prof. Dr. G. JOOS, München, und Prof. Dr. TH. KALUZA fF, Göttingen 


8., verbessert2 Auflage 
XII, 299 Seiten mit 97 Abb. im Text. 1956. Gr. 8°. DM 21.60, Leinen DM 23.10 


Die Technik: Ohne sich in Einzelheiten zu verlieren, wird in neun Abschnitten der 
umfangreiche Stoff behandelt. Bewußt wurde auf streng logische Begründung verzichtet, 
weil man den Zweck des Buches, dem Praktiker die Lösung einer gestellten Aufgabe zu 
ermöglichen, der Darstellung mathematischer Beweise, die in jedem Falle gelten, verzog. 
Daraus ergibt sich die Anordnung des Stoffes, Das Buch ermöglicht dem Leser, das Be- 
arbeitete an Aufgaben zu erproben; denn faßt jeder Abschnitt schließt mit einer Anzahl 
von Aufgaben, deren Lösungen am Schluß des Buches gebracht werden, bisweilen unter 
Anführung des Lösungsweges, Viele Beispiele aus der Praxis fördern das Verständnis, 
Das Buch ist auch als Nachschlagewerk geeignet. Sieben Auflagen innerhalb 16 Jahren 
für ein Mathematikbuch beweisen, daß die Verfasser mit diesem Werk das Richtige ge- 
bracht haben. W. Richter 


Operatorenrechnung 


und Laplacesche Transformation nebst Anwendungen in Physik und Tedinik 
Von Prof. Dr. phil., Dr.-Ing. e. h. K. W. WAGNERT 


2., verbesserte Auflage. XVIII, 471 Seiten mit 126 Abb. i. Text und einem 
Porträt von Laplace. 1950. Gr. 8°. Geb. DM 42.80 


VDI-Zeitschrift: Der Verfasser führt in die ganze Theorie ein, indem er sehr 
wertvolle praktische Beispiele aus verschiedenen Disziplinen der Naturwissenschaft und 
Technik durchrechnet. Man lernt, die Verfahren zu handhaben, aus dem Gefühle heraus, 
daß sie die einzige Möglichkeit bieten, die Probleme anzugreifen. Das Studium des vor- 
liegenden Werkes kann Naturwissenschaftlern, Technikern und Mathematikern wärmstens 
empfohlen werden. U, Wegner 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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